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ABSTRACT 


Several  fundamental  aspects  of  a  unitary  classical  field  theory 
of  extended  elementary  particles  are  studied.  The  theory  considered  is 
based  on  a  Lorentz  invariant  Lagrangian  density,  constructed  from  the 
electromagnetic  four-vector  potential,  A^,  and  a  real  scalar  field,  f. 

The  Lagrangian  is  not  gauge  invariant.  Particle  models  are  found  for 
which  a  general  class  of  localized  solutions  of  the  field  equations  have, 
at  most,  a  discrete  spectrum  of  allowed  electric  charge.  Exact,  finite, 
algebraic  solutions  are  derived  for  some  particular  examples  of  discrete 
charge  systems.  The  rest  energy  of  these  solutions  is  found  to  be  finite 
and  non-negative,  while  the  angular  momentum  is  zero.  The  problem  of  the 
division  of  the  total  angular  momentum  in  classical  field  theory,  into 
orbital  and  intrinsic  components,  is  considered.  For  the  type  of  Lagrang¬ 
ian  studied,  the  integrals  defining  the  orbital  and  spin  angular  momenta 
are  found  to  be  proportional  for  localized  static  solutions.  An  attempt 
is  made  to  determine  the  stability  of  solutions  in  the  models  studied, 
but,  due  to  the  complexity  of  these  systems,  we  are  unable  to  reach  any 
firm  conclusion. 
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NOTATION 


Here  we  outline  some  of  the  conventions  and  notations  that  will 
be  used  throughout  the  thesis. 

Define  the  coordinate  four-vector 

(x^)  =  (xk,  ict)  =  (x,  ict), 

where  ' c'  is  the  speed  of  light  in  vacuum.  Except  as  otherwise  noted, 
Latin  indices,  (i,j,k,  etc.),  take  on  the  values  1,2,3;  Greek  indices, 
(a,$,y,6,  etc.),  take  on  the  values  1,2, 3, 4.  The  summation  convention 
is  used,  so,  for  example: 


T 

yv,v 


T  +  T  =  T 

yk,k  y4,4  yv,v* 


Here, 


we  have  denoted  partial  derivatives  with  a  comma: 

9T 

rp  _  yv  _  . 

x.  -  a  —  o  x  • 

yv,v  9x^  v  yv 

We  will  have  occasion  to  consider  integrals  of  the  form 


P 

U 


* 


a 

where  a  is  a  space-like  surface,  with  a  normal  specif ed  by  the  time¬ 
like  vector,  da  .  If  a  is  a  surface  of  constant  time,  then 

v 


(dav)  =  (0,0,0-|d3x)  , 


where  d3x  =  dx  dy  dz  =  dx^  dx2  dx^. 


In  this  case,  P  becomes 

V 


ix 
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i 

1  ■ 


•  .  vl  1 
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Unless  otherwise  specified,  the  range  of  3-dimensional  integrations 
is  all  space. 

Define  the  electromagnetic  four-potential 

\  ~  i(^  =  i<J>)  » 

in  terms  of  which  one  can  define  the  electromagnetic  field  tensor 


F  =  A  -  A 
pv  v,  y  y,v 


Electric  and  magnetic  fields,  ft  and  ft,  respectively,  are  defined  in  the 
usual  way  by 


* s  -*♦  -I  fi  :  9  =  -  *xt  ■ 


The  electric  charge  density,  p  ,  and  the  current  density,  j^,  are  combined 
into  the  four-vector  current  density,  j  ,  given  by 

y 

j  =  (ft,  icp)  =  ~r~  F 
Jy  4tt  yv,v 

We  adopt  the  following  definition  for  the  canonical  energy- 
momentum  tensor 


diC 


T  =  6  4-  A  .. 

yv  yv  a,y  3A 


-  f 


,y  3f 
a,v  ,v 


where  jC  =  *C( A  ;  f  ;  A  ;f  )  is  a  Lagrangian  density  depending  only 

y  y  >  v 

on  the  vector  field  A  ,  the  scalar  field  f,  and  their  derivatives. 

M 

£  does  not  depend  explicitly  on  the  coordinates  x  ,  so  one  has  T  =  0. 

K  r  J  y  yv,v 

The  conserved  energy-momentum,  P  ,  is  given  by 


(P  )  =  (P.  ,  IE)  E  (?,  IE)  E  ( 

y  k 


T  da  )  , 
yv  v 


where  ft  is  the  three-momentum,  E  is  the  energy,  and  a  is  an  arbitrary 


x 


. 


■ 


•  •• 


P  r; 


-•  i  .  mG 


space-like  surface. 

For  definitions  of  the  symmetrized  energy-momentum  tensor  and  the 
angular  momentum,  the  reader  is  referred  to  Chapter  5.  The  dimension¬ 
less  variables  used  in  the  remainder  of  the  thesis  are  introduced  in 
Chapter  1. 
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CHAPTER  1 


INTRODUCTION 

In  this  thesis,  we  will  study  several  aspects  of  a  classical 
field  theory  of  extended  elementary  particles.  As  will  be  elaborated 
in  the  following,  the  theory  is  of  the  type  known  as  unitary.  In  the 
models  studied,  nonlinear  field  equations  are  derived  from  a  Lorentz 
invariant  Lagrangian  density.  The  Lagrangian  is  not  gauge  invariant, 
but  we  are  still  able  to  define  a  conserved  electric  charge.  We  attempt, 
with  some  limited  success,  to  choose  the  Lagrangian  such  that  particle¬ 
like  solutions  of  the  field  equations  have  only  a  discrete  spectrum  of 
allowed  charges. 

In  a  unitary  theory,  particles  have  no  properties  independent  of 
the  field;  particles  being,  simply,  special  solutions  of  the  field  equa¬ 
tions.  These  particle-like  solutions  are  ones  in  which  there  are  local¬ 
ized  regions  where  the  field  has  a  non-constant  value.  Particles  are 
not  assigned  coordinates  which  satisfy  equations  of  motion  and  are  not 
attributed  with  mechanical  properties  independent  of  the  field.  The 
entire  energy,  linear  momentum,  and  angular  momentum  of  a  particle  are 
found  by  integrating  the  corresponding  quantity  for  the  field  over  the 
region  defining  the  particle. 

In  the  sixty-five  years  since  Mie  first  introduced  a  unitary 

*j- 

field  theory,  many  authors,  including  Born  and  Infeld  (1934);  N.  Rosen 
(1939);  Finkelstein,  Fronsdal,  and  Kaus  (1956);  Schiff  (1962,1969); 

4* 

Mie's  theory  is  discussed  by  Pauli  (1958).  Note,  in  the  body  of  the 
thesis,  references  will  be  referred  to  by  the  author  and  year. 

Complete  references  are  given  at  the  end  of  the  thesis. 
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and  G.  Rosen  (1965),  have  tried  this  approach  to  the  formulation  of  a 
theory  of  extended  elementary  particles.  Part  of  the  appeal  of  a  unitary 
theory  is  aesthetic,  in  that  the  disjoint  notions  of  particle  and  field 
are  reduced  to  the  single  concept  of  the  field,  but  there  are  other  advant¬ 
ages  as  well. 

In  the  Maxwell  theory  of  the  electrodynamics  of  point  particles, 
the  particles  appear  as  singularities  of  the  field.  The  field  equations 
do  not  apply  at  the  position  of  the  particle,  leading  to  the  require¬ 
ment  that  the  field  equations  must  be  supplemented  by  equations  of  motion 
for  the  particle.  In  a  unitary  theory,  one  can  impose  the  condition  that 
the  fields  be  finite  everywhere  for  the  solutions  that  represent  physical 
particles.  The  field  equations  are  satisfied  at  all  points  in  space  and 
therefore,  it  may  be  possible  that  no  additional  equations  are  required 
to  complete  the  description  of  the  system. 

In  a  non-unitary  theory,  in  which  particles  are  assigned  coordin¬ 
ates  which  obey  equations  of  motion,  there  is  difficulty  in  describing 
multiple  particle  processes,  especially  those  involving  a  change  in  the 
number  of  particles.  In  the  region  where  extended  particles  overlap, 
the  particle  concept  loses  its  meaning.  A  unitary  theory  does  not  have 
this  difficulty  because  the  concept  of  a  particle  is  not  essential  to 
the  mathematical  manipulations. 

Unitary  field  models  developed  to  the  present,  have  been  primar¬ 
ily  concerned  with  demonstrating  the  existence  of  particle-like  solu¬ 
tions  which  have  appropriate  interactions  with  other  particles  at  large 
separations.  Classical  theories  appear  to  be  a  considerable  distance 
from  treating  the  problem  of  interactions  at  short  range,  (where  one 
would  expect  to  require  the  input  of  quantum  ideas  in  any  case) . 
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Perhaps  one  of  the  major  reasons  for  the  slow  development  of  the  unitary 
theories,  is  that  the  field  equations  in  such  an  approach  are  necessarily 
nonlinear,  if  an  interacting  system  is  to  be  described.  In  a  linear 
theory,  the  fields  due  to  two  particles  may  be  added  to  give  another 
solution  of  the  field  equations,  but,  in  the  resulting  system,  the  pres¬ 
ence  of  one  particle  does  not  affect  the  other.  Again,  one  must  add 
equations  of  motion  for  the  particles  in  order  to  get  a  complete  descrip¬ 
tion  of  the  system.  On  the  other  hand,  in  a  unitary  theory,  the  inter¬ 
actions  between  particles  are  'built  in'  by  way  of  the  nonlinearities. 

We  now  look  at  some  of  the  other  features  of  the  models  that  we 
have  studied.  The  fields  used  to  constuct  our  particle  models  are  the 
electromagnetic  four-vector  potential  A^,  and  a  scalar  field,  f.  The 
field  f  was  introduced  by  Nash  (1975)  as  a  quantity  derived  from  a  hypo¬ 
thetical  fluctuating  electromagnetic  field  which  is  assumed  to  permeate  all 
space.  A  randomly  fluctuating  background  field  has  also  been  used  by 
Boyer  (1975)  and  Schiff  (1969),  among  others,  in  the  context  of  classical 
field  theories.  The  possible  origin  of  a  background  field  and  the  prop¬ 
erties  it  might  have,  will  not  be  discussed  here,  except  in  relation  to 
the  boundary  conditions  we  impose  on  the  solutions  of  the  field  equations. 
In  a  Vacuum',  where  only  the  background  radiation  is  present,  the  value, 
f^,  of  the  field  f,  is  assumed  to  be  a  constant.  We  allow  for  both  the 
possibilities,  that  is  or  is  not  equal  to  zero. 

The  fields  f  and  A  satisfy  coupled  nonlinear  equations  which  are 

P 

derived  from  a  Lagrangian^,  <C  ,  through  a  principle  of  stationary  action. 


t 


The  quantity  L  = 


JL  d3x  is  the  Lagrangian,  so  ^  is  most  properly 


referred  to  as  the  Lagrangian  density.  However,  for  convenience,  we 
follow  the  common  practice  of  calling  ^  ,  simply,  the  Lagrangian. 
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We  assume  that  «£  depends,  explicitly,  only  on  the  fields  f,  A  ,  and 

y 

their  first-order  derivatives,  f  and  A  ,  respectively.  Define  the 

,  v  y  >  v 

action.  A,  by 


A  = 


-c 


v 


d3x  dt  , 


where  V  is  an  arbitrary  four-dimensional  volume.  The  variational  prin¬ 
ciple,  6A  =  0,  requiring  that  A  be  an  extremum  with  respect  to  varia¬ 
tions,  of  the  fields  f  and  A  ,  which  vanish  on  the  boundary  of  the  volume 

y 

V,  gives  the  Euler-Lagrange  field  equations 


8 

P 


(d<C 


8 

P 


U 

8A 

a 

8*_ 

8f 


(a=l,2,3,4) 


There  are  several  advantages  to  using  a  Lagrangian  based  theory. 
First,  by  requiring  that  the  Lagrangian  be  a  scalar  under  Lorentz  trans¬ 
formations,  the  resulting  field  equations  are  automatically  in  covariant 

form.  The  Lorentz  invariance  of  the  Lagrangian  also  leads,  via  Noether’s 

i 

theorem  ,  directly  to  the  conservation  laws  for  energy,  linear  momentum, 
and  the  total  angular  momentum.  Besides  providing  a  framework  for  the 
discussion  of  conservation  laws,  the  existence  of  a  Lagrangian  can  be 
exploited  in  some  calculations,  for  example,  in  variational  approxima¬ 
tions.  Later,  the  extremum  properties  of  the  (integrated)  Lagrangian 
will  be  used  to  simplify  calculations  of  the  energy  in  the  systems  that 
we  have  studied. 

Another  conservation  law  which,  on  the  basis  of  present  experi¬ 
mental  evidence,  should  be  incorporated  into  any  model  of  elementary 


t 


See,  for  example,  Soper  (1976)  ,  pp  101-108  . 
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particles,  is  that  of  the  conservation  of  electric  charge.  One  way  to 
ensure  that  charge  is  conserved  is  to  require  that  the  Lagrangian  be 

-j* 

invariant  under  a  gauge  transformation.  However,  in  the  models  we  have 
studied,  the  Lagrangian  is  not  gauge  invariant,  so  another  method  must 
be  used.  A  conserved  charge  is  introduced  through  the  field  equations, 
rather  than  through  a  symmetry  of  the  Lagrangian.  The  Lagrangians  we 
consider  are  of  the  form 


(1.1) 


1 

4 


F  2  + 
yv 


1 

2 


(f  )2 


,  v 


where  the  term  jC  may  include  nonlinear  terms  and  terms  which  couple  the 
field  f  to  the  electromagnetic  field.  Part  of  the  set  of  equations  fol¬ 
lowing  from  the  Lagrangian  (1.1)  is  of  the  form 


(1.2) 


F 

yv,v 


We  use  the  equation  (1.2)  to  define  the  electric  current  density  four- 

vector  j  =  (]f,  icp)  ,  where  J  is  the  three-vector  current  density  and 

p  is  the  charge  density.  From  (1.2),  using  the  antisymmetry  property, 

F  =  -F  ,  we  find  that  the  current,  j  ,  satisfies  the  usual  continuity 
yv  vy  y 

equation 


j  =  0  . 

y  »y 


This  implies  that  the  electric  charge,  q,  given  by 


q 


[  j  da 
J  y  y 


9 


is  conserved. 

It  is  interesting  that  one  of  the  major  criticisms  of  Mie’s 


t 


Barut  (1964),  pp  142-144. 
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original  work  in  the  area  of  unitary  field  theory,  was  that  the  theory 
was  not  gauge  invariant.  Born  and  Infeld  (1934)  mention  that  Mie’s 
theory 

M  ...  breaks  down  because  Mie’s  field  equations  have  the 
unacceptable  property,  that  their  solutions  depend  on  the 
absolute  value  of  the  potential.” 

Pauli  (1958,  p  192)  gives  the  following  comments  about  Mie’s  theory: 

"  Once  we  have  found  a  solution  for  the  electrostatic  potential 
f  of  a  material  particle  of  the  required  kind,  <j>  +  const,  will 
not  be  another  solution,  because  the  field  equations  of  Mie’s 
theory  contain  the  absolute  value  of  the  potential.  A  material 
particle  will  therefore  not  be  able  to  exist  in  a  constant 

external  potential  field.  This,  to  us,  seems  to  constitute  a 
very  weighty  argument  against  Mie's  theory.” 

However,  a  constant  external  potential  is  an  idealization  which  may  not 
be  relevant  to  the  physical  world,  so  this  does  not  seem  like  a  strong 
objection,  to  us.  It  is  important  to  explore  the  consequences  of  fore¬ 
going  the  requirement  of  gauge  invariance,  without  prejudice,  until  more 
information  about  this  approach  has  been  accumulated. 

We  now  briefly  outline  the  subjects  that  are  discussed  in  this 
thesis.  In  Chapter  2,  we  examine  the  question  of  whether  or  not  there 
exist  field  equations  in  which  all  the  particle-like  solutions  have  the 
same  charge.  To  illustrate  some  of  the  general  results,  specific  exam¬ 
ples  are  given  of  field  equations  for  which  a  restricted  class  of  solu¬ 
tions  all  have  the  same  charge.  However,  these  solutions  are  demonstr¬ 
ated  to  have  some  unphysical  properties  which  limit  their  usefulness. 

In  Chapter  3,  exact  algebraic  solutions  are  found  for  a  few  particular 
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field  theory  models.  The  energy  of  these  solutions  is  calculated  in 
Chapter  4,  and  is  found  to  be  non-negative  in  all  cases.  In  Chapter  5, 
some  general  comments  are  made  about  angular  momentum  in  classical  field 
theories,  with  particular  emphasis  on  the  problem  of  dividing  the  total 
angular  momentum  into  spin  and  orbital  parts.  These  general  comments 
are  discussed  in  relation  to  the  exact  solutions  mentioned  earlier.  The 
angular  momentum  of  all  the  solutions  is  found  to  be  zero.  Finally, 
Chapter  6  deals  briefly  with  the  problem  of  the  stability  of  solutions 
in  classical  field  theory.  Attempts  to  apply  stability  analysis  to  the 
simplest  of  the  models  studied  in  the  previous  chapters,  met  with  little 


success. 


Throughout  the  thesis,  dimensionless  variables  will  be  used  to 
simplify  the  notation.  Lengths  will  be  measured  in  units  of  a  fundamental 


length,  A.  In  the  case  where  the  asymptotic  value,  f^,  of  the  field  f, 

is  non-zero,  it  defines  a  corresponding  length  A^.  It  would  be  natural 

to  make  the  identification  A  =A  ^ .  However,  in  order  to  allow  for  the 

possibility  that  A  q=0,  we  leave  A  unspecifed.  Charge  will  be  measured  in 
e 

units  of  —  ,  where  e  is  the  magnitude  of  the  electron  charge  and  q~  is  a 

q0  U 


pure  number,  which  can  be  chosen  according  to  the  needs  of  the  particular 
problem  at  hand.  For  example,  given  any  solution  of  the  field  equations, 
having  a  non-zero  charge,  q,  the  number  q^  can  be  chosen  such  that  q  takes 
on  any  desired  value.  The  choice  of  q^,  however,  cannot  affect  the  ratio 
of  the  charges  of  two  different  solutions  of  the  field  equations.  Finally, 
velocities  are  measured  in  units  of  the  speed  of  light,  c. 

Denoting  dimensionless  quantities  with  a  hat  (~),  we  make  the  fol¬ 


U)  =  {(*,), 


lowing  definitions: 


■ 


8 


A  = 


<ln  A, 


y  e  n0  “y  * 


:  _  a2  2  *  :  A2  2  r 

f  =  pqg  f  ,  so,  f  =  ?q2  fQ 


x4  2  y* 

?  qo*C 


Other  quantities  of  interest  are  the  energy  and  the  angular  momentum, 

e2  e2 

which  are  measured  in  units  of  - — ?  and  — ^respectively. 

Xqo  cqo 

In  the  remainder  of  the  thesis,  only  the  dimensionless  quantities: 

A  A  A  A 

x  ,  A  ,  f  ,  f*  ,  etc.  will  be  used  but  for  notational  convenience, 

y  y  0 

we  drop  the  hats  throughout.  Equivalently,  we  may  say  that  e=qQ ,  and 

c=  X  =  1,  in  the  notation  used.  Note  that  the  fundamental  constant,  fi 

q2 

has  the  value  "h  =  JjD  -  137  q2  in  these  units,  (a  being  the  fine-structure 

a 

constant),  so  the  common  natural  units  where  11=1,  correspond  to  the 

h 

choice  q^  =  a  . 


- 


' 


■ 

. 


' 


CHAPTER  2 


FIXED  CHARGE  SYSTEMS 

2 . 1  Introduction  to  Fixed  Charge  Systems 

One  of  the  most  general  principles  in  physics  is  that  the  electric 
charge  of  any  particle  is  a  simple  integer  multiple  of  some  fundamental 
unit.  Also,  the  total  charge  of  any  system  appears  to  be  strictly  con¬ 
served.  The  generality  of  these  observations  is  evidenced  by  the  fact 
that  they  apply  to  all  particles,  independent  of  their  mass,  spin,  and 
the  nature  of  their  interactions  with  other  particles.  For  example, 

•J* 

a  variety  of  experiments  show  that  the  difference  between  the  electron 
and  proton  charges  is  consistent  with  zero  and  has  an  upper  limit  of 
10  19  or  less,  times  the  electron  charge.  That  this  is  so,  despite  the 
fact  that  the  electron  appears  to  be  point-like  and  does  not  participate 
in  the  strong  interactions,  whereas  the  proton  has  structure  and  does 
interact  strongly,  seems  quite  remarkable.  It  is  natural  to  try  to  in¬ 
corporate  the  conservation  of  a  discrete  electric  charge  into  our  part¬ 
icle  models. 

In  the  systems  studied  here,  the  conservation  of  charge  has  a  dy¬ 
namical  origin,  in  that  it  derives  from  the  field  equations  rather  than 
from  the  invariance  of  the  Lagrangian  under  some  symmetry  operation. 
Specifically,  as  discussed  in  Chapter  1,  our  Lagrangian  is  not  gauge  in¬ 
variant.  Therefore,  the  usual  mechanism  for  ensuring  charge  conservation 
is  not  present  here.  In  the  following,  we  investigate  the  possibility 
that,  in  addition  to  ensuring  charge  conservation,  the  field  equations 

1  Dylla  and  King  (1973)  give  a  summary  of  some  of  the  experimental  results. 
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- 
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\1  ,  '  ! 


may  also  give  rise  to  the  discrete  character  of  this  conserved  charge. 
Dirac  (1931,  1948)  has  proposed  a  theory,  based  on  the  existence  of  the 
magnetic  monopole,  which  would  account  for  the  quantization  of  electric 
charge.  However,  experimental  searches  have  not  yet  given  any  firm  evi¬ 
dence  for  the  existence  of  a  magnetic  charge.  Thus,  it  is  possible  that 
the  origin  of  a  discrete  electric  charge  lies  in  another  direction,  such 
as  the  one  that  we  explore  here. 

2.2  The  Class  of  Lagrangian  Studied 


Ideally,  one  would  like  to  guarantee  that  all  solutions  of  the 
field  equations  have  charges  which  are  integer  multiples  of  some  fund¬ 
amental  unit.  There  are  an  unlimited  number  of  possible  Lagrangians 
one  could  consider  as  candidates  for  generating  field  equations  with 
this  property.  However,  our  search  was  limited  to  a  rather  specific 
form  of  Lagrangian,  which  retains  some  similarity  with  classical  electro¬ 
magnetism.  Furthermore,  for  the  most  part,  only  spherically  symmetric, 
time- independent  solutions  were  considered.  If  all  the  solutions  of  a 
certain  class,  are  restricted  by  the  field  equations  to  having  the  charge 
-q,  0,  or  q,  where  q  is  some  fixed  positive  number,  we  will  call  this  a 
fixed  charge  system,  with  respect  to  that  class  of  solutions.  * q f  will 
be  called  the  value  of  the  fixed  charge.  Both  ±q  are  allowed  because 
particles  observed  in  nature  are  seen  in  oppositely  charged  pairs. 

Except  for  some  isolated  examples,  only  Lagrangians  of  the  follow¬ 
ing  form  are  considered: 


(2.1) 


7T  B(f,A  )  (f  )2  +  C(f ,A  ) 
2  y  ,v  y 


u  F  f  +  7-  F  2 
y  yv  ,v  4  yv 


\ 

+  G(A  ,f) 


Here  B,  C,  and  G  are  arbitrary  differentiable  scalar  functions  of  the 


(.  .  ' 

. 


' 


■ 
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real  scalar  field  f  and  the  four-vector  potential  A^ .  Note  that  the 

electromagnetic  term,  F  2  has  the  usual  sign  and  magnitude,  u  is  a 

yv  y 

constant  four-vector.  For  convenience,  we  assume  that  the  function,  C, 

appearing  in  the  Lagrangian,  is  defined  so  that  u  2=-l.  We  will  examine 

P 

localized,  time-independent  systems,  where,  as  indicated  in  Appendix  1, 
the  three-momentum  of  the  system  vanishes,  so  we  may  speak  of  a  particle 


rest  frame.  In  such  a  reference  frame,  it  is  assumed  that  u  reduces  to 

y 

the  form 


(2.2) 


(  uy  )  =  (  o,  0,  0,  iu)  , 


where  u  is  a  real  constant,  equal  to  ±1.  The  definition  of  u  as  a  four- 

y 

vector  then  uniquely  determines  u^  in  any  reference  frame  related  to  the 
rest  frame  by  a  Lorentz  transformation.  It  is  seen  that  u^  is  propor¬ 
tional  to  the  four -velocity  of  the  rest  frame  of  the  particle. 

The  Lagrangian  (2.1)  was  considered  by  Nash  (1975)  in  the  same 
context,  but  in  the  present  work  a  more  systematic  notation  is  developed 
and  some  additional  cases  are  studied.  As  did  Nash,  we  restrict  our 
study  to  the  case  where  the  vector  potential,  A  =  0.  The  field  equations 
for  this  case  have  been  given  by  Nash.  He  found  that  these  equations  are 
considerably  simplified  if  the  parameters  B  and  C,  appearing  in  the  La¬ 
grangian,  depend  only  on  the  scalar  field,  f,  and  are  related  by 


B(f )  =  b 


C(f) 


)  2 


y 


where  fb'  is  an  arbitrary  constant.  With  this  form  for  B  and  C,  and 
the  previously  mentioned  restrictions  on  A  and  u^,  the  Lagrangian  (2.1), 
in  the  time-independent  case,  becomes: 


(2.3) 


JL_ 

4tt 


(CVf)2  +  uC 


1 

2 


(V*)2  +  G(cJ)  ,f ) 


' 


. 


I 


12 


The  corresponding  field  equations  can  be  written 


(2.4) 


using  u2  =  1. 


The  appearance  of  the  combination  C  ^f,  in  the  Lagrangian  (2.3), 
suggests  introducing  a  new  field  h=h(f),  defined  by 


f 


(2.5) 


h(f)  =  C(ff)  df'  ;  so  Vh  =  C  Vf  . 


Assuming  that  the  relation  h=h(f)  can  be  inverted  to  give  f  as  a  function 
of  h,  the  Lagrangian  (2.3)  takes  the  form 


\ 


Note  that  the  coefficients  of  the  derivative  terms  in  the  Lagrangian  (2.6) 
are  all  constants.  Taking  the  field  variables  to  be  <j>  and  h,  the  field 
equations  following  from  this  Lagrangian  take  the  simple  form: 


(b+i)  -  U  f  -  b 

(b+1)  ’2h=  f+ulf  ' 


(2.7) 


If  one  replaces  h  by  its  definition  (2.5)  in  terms  of  f,  these 
equations  are  seen  to  be  the  same  as  those  derived  treating  <j>  and  f  as 
the  field  variables.  Since  the  problem  of  fixing  the  charge  depends 
only  on  the  field  equations,  it  is  not  affected  by  the  change  of  variable 
from  f  to  h.  Therefore,  in  principle  it  appears  that  without  loss  of 
generality,  one  can  restrict  the  search  for  fixed  charge  systems  to 
Lagrangians  where  the  coefficients  of  the  derivative  terms  are  constants, 
as  in  (2.3).  However,,  in  practice  there  are  some  differences  between  the 
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cases  where  B  and  C  are  functions  of  f  and  when  they  are  constants.  This 
difference  results  from  self-imposed  restrictions  on  the  type  of  function 
one  chooses  for  G(<f>,f),  and  on  the  type  of  solutions  one  seeks  to  the 
field  equations.  For  example,  on  the  grounds  of  simplicity  one  might 
restrict  one's  study  to  the  cases  where  G  is  a  polynomial  function  of  its 
arguments.  However,  if  G  is  a  polynomial  in  f  it  is  not  necessary  that 
it  be  one  in  h.  This  is  illustrated  by  the  case  C(f)=f,  for  which  the 
field  h(f)=  ~  f2,  from  the  definition  (2.5).  If  we  choose,  for  example, 
G(4> , f )  =  (f — u<f>) 14 ,  then  G(<j>,f(h))  =  (  (2h)^-u<J>) 4 .  Similarily,  for  tract- 
ability,  one  might  require  that  the  solutions  one  studies  have  an  asymp¬ 
totic,  (r-*°°) ,  expansion  of  some  specific  form.  For  example,  in  a  spher- 

3.  b 

ically  symmetric  case  one  might  assume  that  h  has  the  form  h  =  —  +  — 2  +... 

O  o 

Then,  if  C(f)=f,  as  before,  f  =  (2h) 2  =  (— )  2  +  ...  does  not  have  an 
expansion  of  the  same  form,  i.e.  with  only  integral  powers  of  f. 

In  spite  of  the  preceding  remarks,  we  will  subsequently  consider, 
for  the  most  part,  only  those  Lagrangians  in  which  the  coefficients  of 
the  derivative  terms  are  constants.  Specifically,  the  Lagrangian  used  is 


(2.8) 


4tt 


“  B  (Vf)2  +  uC  Vf-Vcf, 


1 

2 


(v<j>)2  +  G(4>,f) 


9 


where  u,  B,  and  C  are  constant  and  G  depends  on  <j>  only  through  the  scalar 
combinations,  <p2  and  u<}> ,  and  is  a  polynomial  function  of  its  arguments. 
The  field  equations  following  from  the  Lagrangian  (2.8)  are 


(2.9) 


-  V2cf>  +  uC  V2f  = 

d<|> 

uC  V2<J>  +  B  V2f  =  |y  . 


If  B  ^  0,  then  the  Lagrangian  (2.8)  and  the  equations  (2.9),  can  be  simp¬ 
lified  somewhat,  by  introducing  a  new  field  variable,  f,  given  by: 


' 
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(2.10) 


In  terms  of  this  variable,  f,  the  Lagrangian  (2.8)  becomes 


(2.11)  *C  =  B  (Vf)2  -  |  D  (7(.)2  +  G(*,f) 


where 


(2.12) 


and 


(2.13) 


G(<M)  =  G(cj>  ,f ) 


The  reason  for  introducing  the  variable,  f,  is  that  the  cross  term, 
V<f>*Vf  ,  does  not  appear  in  the  transformed  Lagrangian  (2.11).  Treating 
<j>  and  f  as  the  independent  field  variables,  the  field  equations  derived 
from  the  Lagrangian  (2.11)  are 


(2.14) 


assuming  B  ^  0  and  D  #  0.  One  can  easily  show  that  the  equations  (2.14) 
are  equivalent  to  the  system  (2.9).  Therefore,  if  we  find  that  (2.14) 
admits  only  solutions  in  which  the  charge  has  a  particular  value,  then 
the  charge  is  also  fixed  in  the  original  system  (2.9).  Because  the  form 
of  the  equations  (2.14)  is  simpler  than  that  of  (2.9),  we  prefer  to  use 
the  variable  f,  rather  than  f. 


Recall  that  the  constant  was  introduced  as  the  value  of  the 


field  f,  as  r*00.  From  the  definition  (2.10),  we  find  that  f  satisfies 


this  same  boundary  condition  at  infinity:  f  ->  f^ 


as  r  -*•  00 .  We  now  intro- 


I  9 


j  ■  -  -t 
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duce  another  field,  g,  defined  as  that  part  of  f,  which  vanishes  at  spat¬ 
ial  infinity: 


(2.15) 


We  find  it  convenient  to  express  the  Lagrangian  in  terms  of  the  field  g, 
rather  than  f.  Again,  this  transformation  does  not  affect  the  field  equ¬ 
ations  or  the  problem  of  fixing  the  charge.  Because  f^  is  a  fixed  const¬ 
ant,  Vg  =  ^f  ,  and  the  Lagrangian  becomes 


(2.16) 


JL 

4tt 


D  (V<f>)2  +  G(<j),g) 


5 


where 


(2.17) 


G(cj),g)  =  G(cf>,f)  -  G((j>,f) 


The  relationship  between  the  various  forms  of  the  Lagrangian:  (2.8),  (2.11), 
and  (2.16),  can  be  made  more  explicit  by  writing  out  the  equalities  in 
(2.17)  in  an  expanded  form.  Previously,  we  assumed  that  G  was  a  poly¬ 
nomial  in  the  fields  <f>  and  f.  Consequently,  we  can  write  G  as  a  finite 
summation  of  the  form 


(2.18) 


G(<M) 


l  G  f"  . 

m ,  n  mn  Y 


The  summation  ranges  over  all  non-negative  integers  m,  n  for  which  the 

constants  G  are  non-zero.  Because  the  relations  (2.10)  and  (2.15), 
mn 

defining  f  and  g,  respectively,  are  linear,  G  and  G  are  also  polynomials, 
and  are  of  the  same  degree  as  G.  We  can  write 


G(<J>,f) 


l  G 
m,n  mn 


I  G 
m,n  mn 


f” 
,m  n 

<f>  8 


(2.19) 
and 

(2.20) 


G(<f>,g)  = 


■ 


. 


, 
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Using  the  definitions  (2.10)  and  (2.15),  and  using  the  binomial 
theorem  to  expand  powers  of  f  and  f,  one  finds  from  (2.17),  that 


(2.21) 

and 

(2.22) 


mn  n  '  K  0J  mn+k  * 

k=0 


m 


5mn  ■  I nO  <f  >1'  wk 

k=0 

.  V  7  (n+k+D !  ( uCvk  1 

.  -L  nlk’.l!  '  B  '  C  V  m~ 
k=0  1=0 


uC  vk 


k  n+k+1 


where  (n^)  =  •  If  a  fixed  charge  system  is  described  in  terms  of 


n 


the  variables  <f>  and  g,  as  we  will  do  in  the  following,  then  (2.21)  and 
(2.22)  enable  one  to  easily  find  the  Lagrangian  (2.8)  or  (2.11),  describ¬ 
ing  the  same  system  in  terms  of  the  original  variables. 


2.3  Fixing  the  Charge  Through  Asymptotic  Expansions 


We  now  examine  the  possibility  that  the  charge  can  be  fixed  in  a 
system  derived  from  the  Lagrangian 


(2.16) 


£  =.J- 


J  B  (Vg)2  -  |  D  (V4>)2  +  G4,g) 


by  appropriate  choice  of  the  constants  B  and  D,  and  the  polynomial 


(2.20) 


G(4>>g)  =  l  g"  . 

m,n 


The  field  equations  derived  from  the  Lagrangian  (2.16)  are,  for  B,  D  ^  0, 


(2.23) 


I 

D  3<f> 


V  2g  = 


I  iG 

B  3g  * 


In  order  to  clarify  our  approach  to  fixing  the  charge  using  the  field 
equations  (2.23),  we  first  present  a  simple  example  of  a  system  in  which 
all  members  of  a  class  of  spherically  symmetric  solutions,  (if  they  exist). 


_ 


. 


have  the  same  charge. 


For  the  specific  example,  we  choose  B  =  D  =  1,  and 

(2.24)  G(4>,g)  =  X  (g-u<|>)2  +  2  g3  (g-u<|>)  , 

where  X  is  an  arbitrary  non-zero  real  constant  and  u=±l,  is  the  fourth 
component  of  a  constant  four -vector,  introduced  in  (2.2).  The  field 
equations  (2.23)  become 


(2.25) 

V2<p  =  -  ~  =  2  X  u  (g— U(f) )  +  2  u  g3  . 

(2.26) 

V2g  =  —  =  2  X  (g-U(j>)  +  8  g3  -  6  u  <j>  g2  . 

Assume  that  a  solution  to  the  equations  (2.25)  and  (2.26)  exists, 
of  the  form 


(2.27) 

<f>  =  ^  +  r 3  +  I”4  +  ”r 5  +  *  •  * 

a  ,  b  .  c  d  .  e 

8  =  7  +  72  +  73  +  74  +  75  +  •  •  •  » 

in  which  all  the  parameters  appearing,  are  constant.  Substituting  (2.27) 
into  (2.25)  and  equating  the  coefficients  of  like  powers  of  r  on  each 
side,  one  finds  the  conditions 


(2.28.1) 

0  =  a  -  uq 

(2.28.2) 

0  =  b  -  u3 

(2.28.3) 

0  =  2X  (c-uy)  +  2a3 

(2.28.4) 

2u3  =  2X(d-uy)  +  6a2b 

(2.28.5) 

6uy  =  2X(e-ue)  +  6(a2c+ab2) 

6uy  =  2X(e-ue)  +  6(a2c+ab2) 


Similarily,  substitution  of  (2.27)  into  (2.26),  after  some  simplifica¬ 
tion  gives  the  condition: 
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(2.29) 


6c  —  2A (e-ue)  +  6(a2c+ab2)  +  6q2(c-uy)  . 


(  The  first  four  conditions  derived  from  (2.26),  analogous  to  (2.28.1) 
through  (2.28.4),  do  not  give  any  additional  constraints  in  this  example, 
so  we  have  not  written  them  out.  )  Subtracting  equation  (2.28.5)  from 
(2.29)  we  find 


6 (c-uy)  =  6q2 (c-uy)  . 


(2.30) 


But,  from  (2.28.3)  and  (2.28.1),  we  see  that  c-uy  is  non-zero  for  charged 
solutions,  so  that  it  can  be  cancelled  from  (2.30)  to  give 


(2.31) 


Therefore,  in  this  example,  if  a  solution  of  the  form  (2.27)  exists, 
with  non-zero  charge,  q,  then  q  =  ±1.  We  say  that  the  field  equations 
(2.25)  and  (2.26)  are  a  fixed  charge  system. 

The  method  used  to  find  systems  of  this  type  was  simply  to  choose 
a  form  for  the  function  G(<J>,g),  containing  some  arbitrary  parameters. 

We  then  write  down  the  field  equations  and  substitute  power  series  of 
the  form  (2.27),  which  we  hope  represent  actual  solutions  to  the  field 
equations  for  large  r.  Next,  the  parameters  appearing  in  G  are  chosen 
so  that  combining  the  field  equations  leads  to  a  relation  which  limits 
the  charge  to  a  fixed  value,  as  in  (2.31). 


In  treating  a  more  general  case,  where  G  is  an  arbitrary  polynom¬ 


ial  function  of  <p  and  g,  it  is  convenient  to  introduce  a  systematic  not¬ 
ation  which  will  simplify  the  writing  of  the  field  equations  and  the 


' 
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expression  of  our  results.  Define  the  parameters  a  and  8  from  the 

mn  mn 


field  equations,  by 


v2»=  I  L  gn 


(2.32) 


m,n 


mn 


n?  r  ,m  n 

v  g  =  )  o  t  <p  g 

L  mn  T  & 


m,n 


Comparing  (2.32)  to  (2.23),  one  finds  from  (2.20)  that 


a  =  G 

mn  B  m  n+l 


(2.33) 


6 


mn 


(pH-1)  g 
D  m+1  n 


Assume  there  exist  solutions  of  (2.32)  of  the  form 


<P  = 


oo  b 

1  — 
-n 


(2.34) 


n=l  r 

oo  g 


n 


g 


=  l  — 

n=l  r11 


The  charge,  b^,  is  assumed  to  be  non-zero  throughout  this  chapter. 
Define  a  parameter  e  by 


g- 


(2.35) 


z  = 


For  m  =  0,1,2  ...  ,  let 


m 


(2.36) 

(2.37) 


m 


P  =  T  e  a  ,,  =  a  n  +  e  a  ..+  ...+£  a « 

m  ,  Ln  m-k  k  mO  m-1  1  0m 

k=0 

3P  i 

„  m  ,  _  ,  m-1 

S  =  - —  =  a  ,  ,  +  2e  a „  0  o  +  . . .  +  me  a. 
m  9e  m-1  1  m-2  2  Om 


Similar ily , 
(2.38) 


m 


k=0 
T  =  3Qm 

m  JT 


kk 


(2.39) 


• 
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Because  the  field  equations  are  derived  from  a  Lagrangian,  these  para¬ 
meters  are  not  all  independent.  The  relations  (2.33)  can  be  used  to  show 
that 


(2.40)  mBP  =  --DT  +  eBS 

m  m  m 

To  gain  some  familiarity  with  these  definitions,  we  evaluate  the 

various  parameters  for  the  example  (2.24).  From  the  field  equations 

(2. 25) and  (2.26)  one  can  easily  read  off  the  values  of  a  and  8 

mn  mn 

using  their  definition  (2.32).  One  finds 


^10  ^  2uA.  ;  8hq  2u  » 


'01 


(2.41) 


03 


~2uX  ,  o&q ^  2A  ,  01q3  ®  *  ot  12  - 6u 


The  formulae  (2.36)  through  (2.39)  then  give 


s  2(e-u)A  =  uQ^  ;  =  2A  =  uT^  ,  (using  u2=l)  ; 

(2.42)  P3  =  2e2  (4e--3u)  ;  S3  =  12e(2e-u)  ; 

Q3  =  2ue3  ;  T3  =  6ue2 

Using  (2.42)  one  can  verify  that  the  relations  (2.40)  are  satis¬ 
fied.  Perhaps  the  most  helpful  property  to  remember,  when  dealing  with 

the  parameters  P  ,  Q_>  S  ,  and  T  ,  is  that  they  are  derived  from  the 

m  m  m  m 

k  1 

coefficients  of  terms,  <j>  g  ,  in  the  field  equations,  such  that  the  sub¬ 
script  m  =  k+i.  Therefore,  from  the  field  equations  (2.25)  and  (2.26), 
for  example,  since  there  are  no  terms  of  the  form  <j>2,  <j>g,  or  g2  appearing, 
one  can  immediately  conclude  that  P^  =  =  S3  =  =  0. 


Employing  the  notation  which  has  been  introduced,  we  now  summarize 


. 


. 
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the  results  of  our  search  for  polynomial  forms  for  G(<j>,g)  which  give 
fixed  charge  systems.  Some  of  the  details  of  the  investigation  are  pre¬ 
sented  in  Appendix  2.  There,  the  power  series  expansions (2. 34)  are  sub¬ 
stituted  into  the  field  equations  (2.32).  By  equating  coefficients  of 
like  powers  of  r  on  each  side  of  the  equations,  conditions  are  obtained, 
analogous  to  (2.28)  and  (2.29).  These  conditions  are  given  by  (A2.3) 
and  (A2.4).  Next,  using  these  conditions,  it  was  shown  that  the  require¬ 
ment  that  the  energy  of  the  solutions  (2.34)  be  finite,  gives  only  the 
restriction  =  0,  on  the  form  of  the  function  G(<J>,g)  given  by  (2.20). 
Finally,  conditions  are  established  for  which  the  charge  is  fixed  by  the 
field  equations  to  a  specific  value  for  all  solutions  of  the  assumed  form 
(2.34). 

The  details  of  how  the  value  of  the  fixed  charge  relates  to  the 
parameters  in  the  Lagrangian,  depends  on  the  form  of  the  polynomial  G. 
However,  there  is  one  feature  that  is  common  to  all  the  fixed  charge  sys¬ 
tems  that  were  found,  namely,  that 


(A2.10) 


D  =  e2  B 


apparently  must  be  satisfied  if  the  magnitude  of  the  charge  is  to  have  a 
single  allowed  non-zero  value.  Recall  that  B  and  D  are  the  constant 
coefficients  of  the  derivative  terms  in  the  Lagrangian 


(2.16) 


_A_ 

4tt 


\  B  (h)2  -  j  D  ( V<f> )2  +  G(<J,,g) 


9 


8 1 

and  e  =  •—  is  the  ratio  of  the  leading  terms  in  the  expansions  of  g  and 
bi 

<J) .  The  form  of  G  determines  the  value  of  e.  Note  that  if  g  goes  to  zero 
faster  than  r-1  as  r  -*»,  so  that  g^  =  0,  then  for  charged  solutions,  e=0. 
This  in  turn,  from  (A2.10),  requires  D  =  0.  Since  D  is  the  coefficient 
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of  the  '  electromagnetic  '  tern,  (^)2,  in  the  Lagrangian,  its  vanishing 
would  be  a  significant  departure  from  our  requirement  that  the  present 
theory  be  a  modification  of  classical  electrodynamics.  One  might  expect 
that  the  case  where  D  ^  0,  so  e  ^  0,  would  be  of  greater  relevance  in 
describing  physical  particles.  However,  as  will  be  demonstrated  later, 
this  type  of  model  exhibits  unphysical  behaviour  for  the  interactions 
of  the  charged  particle  solutions,  regardless  of  the  value  of  e. 

We  now  present  the  specific  results  and  illustrate  with  examples. 
Case  IA  :  ^  0,  ^  0 

If  the  parameter  =  3q^>  is  non-zero,  then  the  following  cond¬ 
itions  apparently  must  be  satisfied  if  spherically  symmetric  solutions 
of  the  form  (2.34),  with  finite  energy,  are  to  exist  with  a  non-zero 
charge,  fixed  by  the  field  equations  to  a  single  magnitude,  q.  (These 
results  are  derived  in  Appendix  2) 


(2.43) 


a00  P00  °  ’  a10  ef310  ;  si  eTl  ;  P2  eQ2  ; 


T1P3+T2P2  =  S1Q3+S2Q2  * 


where  z  is  given  by 


(2.44) 


If  Q2  4  0,  then  the  value  of  the  fixed  charge,  q,  is  given  by 


(A2.ll) 


,2  - 


-2Q, 


^2 


(TlP4-SlQ4)+(T2P3-S2Q3)+(T3P2-S3Q2)+(Ila02-Sl602)T5‘ 


An  example  of  a  system  of  this  type,  satisfying  all  the  conditions 


(2.43),  is  the  Lagrangian: 


'  -411 
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X  =  -f 

4tt 


\  (Vg)2  -  |  W)2  +  G 


(2.45) 


,  i.e.  B~D=1  , 


with  G  -  (a(g-u<j>)+Bg2)2  +  ^q2  83  (g-U(^)  , 


where  u  =  ±1,  is  the  fourth  component  of  a  four-vector,  given  by  (2.2), 
and  a  and  3  are  arbitrary  non-zero  constants. 

The  field  equations  are 


V2<J)  2ua  (a(g-u<j>)+Bg2)  +  ^283 
V2g  =  2a2 (g-u<j>)  +  6a3g2  -  4ua3<j>g  +  4(B2+^~2')g3  -  g2  . 
One  finds  that 


S1  =  uT  =  2a2  ; 


(2.46) 


P2  ”  u(^2  =  2ot^  »  S2=8ua3  ;  T2=4a3  ;  aQ2=6a3  ;  3q2=s2u«3  ; 

P3  =  2U(  5^  +  2$2  }  ;  s3  -  12(  5^  +  62  )  ; 


Q 


3 


2 

5qT  5 


6u 

5?  * 


From  the  relations  (2.46)  one  can  verify  that  the  conditions  (2.43)  are 
satisfied,  as  is  (A2.10).  Evaluation  of  (A2.ll)  shows  that  the  parameter 
q,  appearing  in  the  Lagrangian  (2.45),  is  the  value  of  the  fixed  charge. 

As  a  supplementary  comment  to  this  example  and  the  others  which 
will  be  presented,  we  give  an  integral  expression  for  the  rest  energy  of 
solutions  of  the  field  equations.  From  Appendix  1,  choosing  the  parameter 
y  =  1  in  (Al.ll),  we  find  that  the  rest  energy  is 


E 


4tt 


4G  -  <  +  «ti> 


(g-u<t>)  (a(g-u<|>)+3g2)  d3x  . 


We  are  not  able  to  make  a  definite  statement  about  the  sign  of  the  energy 


. 


fp 

. 

. 
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in  this  example,  without  finding  exact  solutions. 

Case  IB  :  f  0,  Q9  =  0 

The  previous  results  applied  only  when  Q2  t  0.  When  Q2  =  0  it 
is  still  possible  to  fix  the  charge  if  the  following  conditions  are 
satisfied  with  #  0. 


(2.47) 


aoo  ”  600 


=  0 


5 


where  as  before,  e  is  given  by  (2.44).  The  value  of  the  fixed  charge,  q, 
is  given  by 


(A2.13) 


,2  = 


-6Q. 


T1(P5-eQ5)  -  3Q2 


No  fixed  charge  system  could  be  found  with  =  0. 

An  example  of  a  system  satisfying  the  conditions  (2.47)  is 


JC 

4tt 


(2.48) 


~  (Vg)2  -  ~  (V4>)2  +  G 


with  G  = 


1 


a(g-u4>)2  +  (g2-<f>2)  ((~V+ 8 )<P2  +  (~T  ~  3 ) g 2 ) 


where,  as  before,  u=±l,  is  the  fourth  component  of  a  constant  four-vector, 
and  a  and  8  are  arbitrary  non-zero  constants.  As  indicated  by  the  nota¬ 
tion,  the  constant  q  appearing  in  G,  will  turn  out  to  be  the  value  of  the 
fixed  charge  for  spherically  symmetric  systems,  as  given  by  (A2.13). 


The  field  equations  in  this  example  are 


V2<}>  =  ua(g-ucj))  +  2<J)  ( C^r  +  8)<J>2~8g2) 
V2g  =  a(g-u<j>)  +  2g(8<J>2+C^-- 8)g2) 


(2.49) 
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From  the  field  equations  we  find  that 


e  =  u 


Tj  =  uS^  =  ua 


(2.50) 


P3  =  uQ3  =  F  !  T3  =  "4uB  ;  S3  =  ■?"  4S  • 


Remembering  that  u2=l,  one  can  verify  that  the  parameters  (2.50)  satisfy 


the  requirements  (2.47). 

In  this  example,  one  finds  from  (Al.ll),  (choosing  y=l),  that  the 
rest  energy  of  localized  solutions  of  the  field  equations  (2.49)  is 


d3x  =  ~  (g-u<j))2  d3x 


Therefore,  by  restricting  a  to  positive  values,  we  can  ensure  that  the 
the  rest  energy  is  non-negative. 

Another  example  of  the  case  IB,  (i.e.  T^O,  C^O) ,  is  provided  by 
the  polynomial  G(<f>,g)  given  by  (2.24),  which  is  the  example  that  was  used 
to  introduce  the  section  on  fixed  charges. 

Case  II  :  T^  =  0 

The  characteristic  feature  of  the  preceding  examples  is  the  appear¬ 
ance  of  linear  terms  in  the  fields  <}>  and  g,  on  the  right  hand  side  of  the 
field  equations  (2.32).  When  T^  =  0,  these  terms  do  not  appear  in  fixed 
charge  systems.  In  these  systems,  the  following  conditions  are  satisfied. 


We  also  assume  that  ^  0.  The  derivation  of  (2.52)  is  similar  to  the 
analysis  of  Appendix  2  for  Case  I. 


The  parameter  e  is  determined  from  the  conditions  P2=Q2=0*  Even 
though  these  equations  are  quadratic  in  c,  (see  definitions  (2.36)  and 


- 

■ 


:  .  .  ' 


■ 


26 


(2.38)),  the  two  equations  cannot  have  two  distinct  real  roots  in  common, 
if  =  e^2*  ^as  as  necessary  to  fix  the  charge).  This  is  a  result  of  a 
relationship  between  and  Q2  arising  from  the  fact  that  the  field  equa¬ 
tions  are  derived  from  a  Lagrangian.  Therefore,  e  is  uniquely  determined 
by  the  Lagrangian  for  both  Case  I  and  Case  II. 

The  value,  q,  of  the  fixed  charge,  for  spherically  symmetric  solu¬ 
tions  of  a  system  satisfying  the  conditions  (2.52)  is  given  by 


(2.53) 


,2  = 


-2Q, 


WeV  -  2^3 


An  example  of  this  type  of  fixed  charge  system  is  given  by 


£  -l 


—  (Vg)2  -  \  (VO)2  +  G 


(2.54) 


with  G  =  ag(g-uO)2  +  g3(g-u<f>) 


where  u=±l,  and  a  and  q  are  arbitrary  non-zero  constants. 

The  field  equations  derived  from  the  Lagrangian  (2.54)  are 


4u  ,3 


V20  =  2uotg(g-u<jO  +  g 


V2g  =  a(3g2-4u«{>g+<j>2)  +  (4g3-3u0g2) 


From  these  equations,  one  finds 


(2.55) 


£  =  U 


4u 

P3  -  uQ3  =  5^ 


;  S2  =  uT2  =  2ua 


S3  5q^ 


12 


5q 


7  u 


The  parameters  (2.55)  satisfy  the  relations  (2.52),  necessary  to  fix  the 
charge,  and  the  value,  q,  of  the  resulting  fixed  charge,  is  correctly 

given  by  (2. 53) . 


' 
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With  Y  -  in  (Al.ll),  we  find  that  the  rest  energy  for  this 


example  is 


(2.56) 


Alternatively,  with  y  =  1  in  (Al.ll),  we  find 


(2.57) 


r 


4G  "  ((*)ff+8if)  d3x  =  g(g~u4>)?  d3x  . 


From  neither  (2.56),  nor  (2.57),  can  we  make  a  definite  statement  about 
the  sign  of  the  rest  energy,  without  looking  at  specific  solutions  of  the 
field  equations. 

2 . 4  Forces  Between  Particles 

All  of  the  fixed  charge  systems  that  we  have  found  satisfy  the 
constraint 

(A2.10)  D  =  e2  B  . 

This  relation  has  important  consequences  with  regard  to  the  interactions 
between  particles.  Before  discussing  the  nature  of  the  forces  between 
particles,  there  are  some  preliminary  remarks  we  would  like  to  make  about 
the  solutions  to  the  field  equations  derived  from  the  Lagrangian  (2.16). 

We  are  able  to  calculate  these  forces  only  for  certain  pairs  of  solutions. 
The  reason  for  this  restriction  is  discussed  in  the  following. 

Some  of  the  systems  described  by  a  Lagrangian  of  the  form  (2.16) 
have  the  undesirable  feature  of  depending  on  the  sign  of  the  charge,  in 
that  the  parameter  u  must  take  on  two  different  values  to  ensure  that  an 
’ antiparticle '  solution  exists  for  each  particle-like  solution.  If  (<J>=<f>0, 
g=gQ,  u=uQ)  is  a  solution,  then,  in  general,  ((J)*-^,  g=gQ)  will  be  a  solu- 
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tion  only  if  u=-uQ.  The  examples  given  for  Case  IA  and  Case  II,  have  this 
property.  We  note  that  the  solutions  ((f>0,g0,u0)  and  (“<l>0»g0»-u0)  have  the 
same  energy;  a  property  observed  for  particle,  antiparticle  pairs  occur¬ 
ring  in  nature.  This  follows  from  the  expression  (A1.12)  for  the  energy, 
showing  that  it  depends  only  on  G,  which  in  turn  depends  only  on  <j>2,  u<}>, 
and  g,  by  our  initial  assumptions. 

The  change  in  sign  of  u,  in  the  Lagrangian,  in  order  to  incorporate 
charges  of  both  sign  into  the  model,  causes  difficulty  in  discussing  the 
interactions  between  oppositely  charged  particles.  In  the  Case  IB  ex¬ 
ample  (2.48),  however,  the  problem  is  avoided,  because  we  can  describe 
charges  of  both  sign,  without  changing  a  parameter  in  the  Lagrangian.  If 
(<})=c{)q,  g=gg)  is  a  solution  to  the  system  (2.49),  then  (<|>=-<J>q,  g=-g^)  is 
also  a  solution  with  the  same  energy,  (2.51),  but  having  the  opposite 
charge . 

Now  we  proceed  to  calculate  the  force  between  two  particles,  which 

f 

are  described  by  the  same  Lagrangian.  Thus,  for  systems  similar  to  the 
examples  of  Case  IA  and  Case  II,  we  consider  only  the  interaction  between 
identical  particles,  but  for  examples  similar  to  that  for  Case  IB,  the 
analysis  will  also  apply  to  particles  of  opposite  sign. 

A  particle  as  understood  here,  is  an  extended,  although  localized 
region  with  finite  energy,  where  the  fields,  A  and  g,  have  finite,  non- 
zero  values.  We  suppose  the  size  of  this  object  to  be  characterized  by 
a  distance,  a.  In  the  asymptotic  region,  at  distances  r>>a,  from  the 
particle,  the  nonlinear  terms  in  the  field  equations  become  small  comp¬ 
ared  to  the  linear  ones.  In  this  region  one  can  linearly  superimpose 


A  similar  calculation,  but  for  a  different  field  theory,  is  given,  for 
example,  by  Rosen  and  Rosenstock  (1952). 


' 

' 


the  fields  from  more  than  one  particle.  However,  we  have  no  guarantee 
that  this  linear  superposition  is  also  a  solution  to  the  field  equations 
in  the  region  where  the  nonlinear  terms  cannot  be  neglected.  One  might 
even  expect  multiple  particle  solutions  not  to  exist  in  the  'fixed* 
charge  systems  we  have  studied.  However,  one  must  remember  that  the 
charge  was  fixed  only  for  spherically  symmetric,  static  solutions,  where¬ 
as  systems  having  more  than  one  particle  present  would  almost  certainly 
be  non-spherical  and  time  dependent.  Whatever  the  actual  situation,  in 
this  section  we  will  assume  that  solutions  do  exist,  representing  more 
than  a  single  particle. 

Suppose  that  there  are  two  particles,  P+  and  P  ,  instantaneously 

at  rest  at  the  points  z-d  and  z=-d,  respectively,  along  the  z  axis.  At 

a  point  r(p,^p,z),  specified  in  cylindrical  coordinates,  far  from  both 

4'  — 

particles,  the  fields  due  to  P  and  P  are  approximately 


(2.58) 


and  4> 


9 


r  r 

where  r±  =  |  r  +  d  f  |»a,  T  being  a  unit  vector  in  the  z  direction. 

Z  2 

» y 

The  total  fields  at  r  are  linear  superpositions  of  the  separate  contri¬ 
butions  of  the  two  particles: 


(2.59) 


<f>(r)  =  >+(r)  +  <J)  (r) 
g(r)  =  g+(r)  +  g  (r)  . 


The  force,  acting  on  a  volume,  V,  enclosing  the  particle  P  , 
say,  is  the  time  rate  of  change  of  the  momentum,  P  ,  in  the  volume. 


From  Appendix  1  we  have: 


■/; 


'  . 
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(A1.5) 


pk  =  -i 


Tk4d3x  * 


where  T  is  the  canonical  energy-momentum  tensor.  Therefore, 


dP 

F  =  — - 
k  "  dt 


dx, 


Tk4d3x  “ 


Tk4,4d3x  ■- 


T  .  d  3X. 


V 


The  last  step  follows  from  the  momentum  conservation  law  T,  =0. 
Converting  the  last  integral  to  a  surface  integral  by  the  divergence 
theorem  gives 


(2.60) 


Fk=- 


T.  .  dS. 

kj  J 


Here  dS^  is  the  outward  normal  to  the  surface,  S,  enclosing  the  volume, 

V.  Taking  the  surface,  S,  to  be  the  xy  plane,  joined  to  an  infinite 
hemisphere  enclosing  the  space  z>0,  we  find  the  only  contribution  to  the 
integral  comes  from  the  xy  plane,  because  T^dS^  vanishes  on  the  hemi¬ 
sphere.  On  the  xy  plane,  the  outward  normal  is  in  the  negative  z  direc¬ 


tion  so  (dS_.)  =  (0,  0,  -  dx  dy) ,  giving  for  the  z  component  of  the  force: 


(2.61) 


F  = 
z 


T33  dx  dy  = 


xy  plane 


2tt 


0  0 


T33^Z==0,p  P  dP 


F  and  F  vanish  from  symmetry  considerations.  To  evaluate  F  we  need 
x  y  z 

the  expression  for  the  canonical  energy-momentum  tensor.  This  is  given 

by 


(2.62) 


T .  .  -  6  . 
ij  iJ 


*  -  g  *£ 
,i  3<J>  ,i  3g  , 

>  J  »  J 


-  6  .  .«C  +  T~ 
ij 


Bg  g  -D4>  <|> 

>  L  » J  »  ^  » J 


where  ^  is  given  by  (2.16). 


Evaluating  T^  from  (2.62),  using  the  fields  <j>(r)  and  g(r)  given 


by  (2.59),  we  find  from  (2.60): 


i 


r 


■ 


P  j 


■ 
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(2.63) 


Fz  4 


B(g{-g”)2  -  D(b|-b~)2 


0 


(p2+d2) ^ 


CO 

I  (  b(sK)2  -  D(bt+bI>2 )  f  7 

o 


_L 

4tt 


2tt 


0  0 


G($>§)  I  Z-Q  P  dP  dfc 


btb 


2tt 


(D-E2B)  .  |  j  G(4,,g)|z=0  p  dpdf  , 


0  0 


+  .  + 


_j_  _r  _  —  + 

where,  in  the  last  line  we  have  written  g^  =  eb^  and  g  =ebj.  In  the 
case  of  classical  electrodynamics,  D  =  1,  B  =  0,  and  G  =  0,  so  we  get  the 
usual  Coulomb  force  law: 


F  = 
z 


bX 


(2d) 


However,  in  the  case  of  the.  fixed  charge  systems,  the  condition 


(A2.10) 


D  =  e2  B 


must  be  satisfied,  so  that  the  'inverse  square'  term  vanishes 

Therefore,  the  expression  (2.63)  for  the  force  reduces  to 

2  7T  0° 


(2.64) 


F 

z 


_L 

4  tt 


G(<J>,g)  lz_Q  P  dP  d{?  • 


0  0 


identically. 


Previously  we  had  required  the  r  0  term  in  G(<P(r)  ,g (r) )  to  vanish 
to  ensure  that  the  energy  of  the  particle-like  solutions  is  finite,  (see 
Appendix  2).  However,  investigation  shows  that  for  the  spherically  sym¬ 
metric  particles  we  have  been  considering,  the  r-4  term  vanishes  as  well. 
This  leads  to  the  result  that  the  force  (2.64)  falls  off  faster  than  the 


t 


For  fixed  charge  systems,  e  is  uniquely  determined  by  parameters  ap¬ 
pearing  in  the  Lagrangian.  Therefore,  e  will  be  the  same  for  parti¬ 
cles  P*  and  P“,  because  they  are  described  by  the  same  Lagrangian. 
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inverse  square  of  the  separation  of  the  two  particles. 

Because  the  charged  solutions  to  the  fixed  charge  systems  studied 
show  no  Coulomb  interactions,  it  was  felt  that  further  investigation  of 
such  systems  is  not  warranted  at  this  time.  No  attempt  was  made  to  find 
exact  or  numerical  solutions  to  field  equations  of  the  type  shown  in  the 
examples . 

2.5  Other  Aspects  of  the  Charge  Fixing  Problem 

In  the  preceding  sections  we  examined  the  possibility  of  finding 
a  Lagrangian  for  which  the  time-independent,  spherically  symmetric  solu¬ 
tions,  with  the  vector  potential,  it  =  0,  all  have  the  same  magnitude  of 
charge.  Even  for  this  restricted  class  of  solutions,  no  completely  ac¬ 
ceptable  model  was  found.  By  further  restricting  the  type  of  solutions 
for  which  one  demands  that  the  charge  be  fixed,  it  is  possible  to  over¬ 
come  the  problem  of  the  vanishing  of  the  Coulomb  force.  However,  in  do¬ 
ing  so,  values  may  appear  in  the  spectrum  of  allowed  charges,  that  are 
not  integer  multiples  of  a  fundamental  unit.  This  will  be  illustrated 


shortly. 


Most  of  the  remainder  of  the  thesis  will  be  based  on  the  Lagrang¬ 


ian 


where  K  and  M  are  differentiable  functions  of  g,  and  6  is  a  constant. 
We  assume  that  M(g)  contains  only  powers  of  g  >  4. 


In  this  section  we  consider  only  the  static  case,  with  X  =  0,  in 


which  case  (2.65)  reduces  to 


(Vg)2  -  Y  (V<|>)2  +  4>2g2  +  M 


(2.66) 


■ 

> 


. 
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The  field  equations  following  from  this  Lagrangian  are 


(2.67) 


V2(}>  =  -2g2(f) 


(2.68) 


V2g  =  2cj)2g  + 


dM 

dg 


Inserting  the  power  series 


oo  b 

♦  -  I  -= 


n=l  r 


n 


;  8n 

8  =  in 

n=l  r 

into  the  equation  (2.67)  gives  the  condition  g^  =  0.  For  the  solutions 

in  which  the  charge  is  fixed  at  the  value  1  by  (2.68).  However, 

in  general,  the  value  of  the  fixed  charge,  q=b^,  is  given  by  q2  =  %n(n-l), 

where  n  is  the  smallest  integer  such  that  g  ^0.  Note  that  the  value  of 

n 

the  charge  does  not  depend  on  the  particular  value  of  g  ,  but  is  deter- 

n 

t 

mined  by  the  general  form  of  g.  Nash  (1975)  has  searched  for  solutions 
to  some  systems  of  the  form  (2.67)  and  (2.68),  but  none  was  found,  except 
in  the  case  n=2 .  If  solutions  also  exist  for  n^2,  then  the  charges  will 
likely  not  be  integer  multiples  of  a  fundamental  unit,  q  being  irrational 

•k 

for  'most'  values  of  n. 

The  system  (2.66)  is  sufficiently  simple  that  the  charge  can  also 
be  examined  for  some  non-spherical  solutions.  Assume  that  series  expan¬ 


sions  exist  of  the  form: 


+  If  one  had  a  more  fundamental  theory  of  the  origin  of  the  scalar  field 
g,  then  it  is  possible  that  the  asymptotic  behaviour  of  g  is  completely 
determined.  For  example,  one  might  find  that  g  constant* r"2  as  r*». 
In  such  a  case  the  charge  would  be  fixed  to  one  value,  rather  than  to  a 
spectrum  of  values  corresponding  to  different  n. 

The  smallest  value  of  n  >  2  such  that  q  is  rational  is  n  =  9  (  q  =  6). 


' 
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(2.69) 


+  .  y  bJ(9>)  ■  ..  ?  Sn<6’*> 

L ,  n  *  6  n 

n=l  r  n=l  r 


We  assume  that  b^(0,%)  and  ( 0 » ^»>)  are  regular  functions  of  the  angles 
0  and  ^  . 

2 

Denote  by  L  ,  the  angular  dependent  part  of  the  Laplacian  opera- 

2 

tor,  v  ,  in  spherical  coordinates: 


(2.70) 

and 


L2  H 


sinB  30 


f  sln0  ] 

0  36  '  +  sine 


(2.71) 


V2*  =  ±  ^ 

y  r  3  r2  r2 


for  any  twice  differentiable  function  \p(r,d, £) . 

Using  the  expansions  (2.69)  in  the  field  equations  (2.67)  and 
(2.68),  and  equating  the  coefficients  of  r~3  on  each  side,  one  finds  the 
condition 

(2.72)  L2  bx  =  2g2b1  . 

In  the  spherically  symmetric  case,  b^  is  a  constant  so  that  L2b^  =  0. 
Therefore,  for  charged  solutions,  where  b^O,  we  find  from  (2.72)  that 
g^=0,  as  has  been  indicated  earlier.  However,  the  equations  do  not  ap¬ 
pear  to  prevent  b^  from  having  some  angular  dependence.  This  case  is 
more  difficult  to  treat  and  may  be  subject  to  difficulties  of  interpret¬ 
ation  in  terms  of  the  charged  particles  that  are  seen  in  nature.  There 
fore,  in  the  subsequent  analysis,  we  assume  that  b^  is  a  constant  and 
consequently,  g^=0.  With  these  assumptions,  equating  coefficients  of  like 
powers  of  r  on  each  side  of  the  field  equations  gives  the  conditions: 
(2.73.1)  (  2  -  L2  )  b2  =  0 

(  6  -  L2  )  b3  =  -2g2b1 


(2.73.2) 


• 
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and 

(2-74.1)  (  2  -  L2  )  g2  =  2g2b2 

(2.74.2)  (  6  -  L2  )  g3  =  2S3bj  +  482blb2  ' 

2 

Equation  (2.74.1)  is  an  eigenvalue  problem  for  the  operator  L  of  the  form: 

(2.75)  L2  g2  =  X  g2  ,  with  X  =  2(l-b2)  . 

There  is  a  regular,  non-zero  solution  for  (2.75)  only  if  the  eigenvalue, X 

f 

has  the  value 

(2.76)  X  =  1(1+1)  ,  with  1  =  0,  1,  2,  ...  . 

However,  to  get  a  real  solution  for  the  charge,  b^,  we  are  restricted, 
here,  to  the  values  1=0  and  1=1.  Corresponding  to  the  value  1=0,  the 
charge  is  given  by  (2.75)  and  (2.76)  as  b^  =  ±1,  whereas  1=1  corresponds 
to  the  neutral  case  b^  =  0.  More  generally,  if  n  is  the  smallest  value 
for  which  4  0,  then 

L2  g  =  (  n(n-l)  -  2b2  )g 
ii  1  J  n 

giving  the  condition 

(2.77)  b2  =  j  (  n(n-l)  -  1(1+1)  )  ;  1  =  0,  1,  2,  .  .  .  ,  n-1  . 

Note  that  when  1  =  n-2,  (n>l),  we  get  b2  =  n-1,  so  that  the  possiblility 
exists  of  finding  solutions  with  the  square  of  the  charge  equal  to  any 
non-negative  integer.  No  non-spherical  solutions  of  the  system  (2.67), 
(2.68)  have  been  found  (  except  for  translations  of  spherically  symmetric 
solutions)  so  it  is  not  known  whether  solutions  with  these  non-integral 

**  See,  for  example,  Merzbacher  (1970),  p  182. 
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values  of  the  charge  are  actually  realized. 

In  the  system  just  considered,  the  charge  may  take  on  some  irra¬ 
tional  values,  even  when  restricted  to  spherically  symmetric  solutions. 
We  now  give  an  example  where  the  charge  spectrum  allowed  for  spherically 
symmetric  solutions  consists  only  of  integer  values.  The  Lagrangian 
density  is 


(2.78) 


•C  -i 


\  u<K?g)2  -  | 


(V<f>)2  +  u<J>^g2  +  yg( 


where  u  is  the  constant  fourth  component  of  a  four-vector  given  by  (2.2) 
and  y  is  an  arbitrary  constant. 

The  field  equations  following  from  the  Lagrangian  (2.78)  are 


V2<j>  +  |u(7g)  2  =  — |u<j>2g2 
<j>V2g  +  V<J>«Vg  =  <J)^g  + 

g 

One  finds  that,  if  g  goes  as  —  as  r-*»,  with  g^O,  then  the  charge  b^ 

r 

is  fixed  by  the  field  equations,  by  b^2  =  n2  ,  n=l,2,3...  .  Unfortun¬ 
ately,  as  in  the  preceding  example,  one  finds  that  angular  dependent 
solutions  apparently  allow  a  larger  spectrum  of  charges.  In  the  present 
example,  one  finds  b2  =  n2  -  1(1+1)  ;  1=0, 1,2,  ...  n-1,  if  non-spherical 

solutions  exist.  In  particular,  when  l=n-l,  we  findb^  =n.  Therefore, 
the  charge  spectrum  for  the  Lagrangian  (2.78)  is  the  same  as  in  the  case 
(2.67)  providing  solutions  exist  corresponding  to  each  n  and  1. 

These  examples  show  that  it  is  not  adequate  to  consider  only  spher¬ 
ical  solutions  in  trying  to  fix  the  charge.  Although  the  Lagrangians  were 
chosen  so  as  to  generate  a  discrete  charge  spectrum  in  the  spherically 
symmetric  case,  it  was  found  that  if  non— spherical  solutions  exist,  they 
are  also  confined  to  a  discrete  set  of  charges.  It  would  be  interesting 


. 

■ 

. 
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to  investigate  the  relationship,  between  the  two  types  of  solution,  in 
more  detail.  Another  problem  for  future  study  is  to  determine  in  what 
manner  the  charge  allowed  for  time-dependent  solutions  is  related  to  the 
charge  spectrum  for  the  static  case.  It  is  of  some  intrinsic  interest 
to  find  a  system  of  field  equations  in  which  the  solutions,  representing 
single  particles  at  rest,  all  have  the  same  charge.  However,  to  fully 
exploit  the  idea  of  fixing  the  charge  dynamically,  through  the  field  equ¬ 
ations,  it  would  seem  necessary  to  extend  the  treatment  to  the  time- 
dependent  case.  We  feel  that  the  resulting  model  would  then  be  much 
more  likely  to  have  relevance  to  the  physical  world  of  time-dependent 
problems  involving  multiple  particle  processes. 


, 


' 
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CHAPTER  3 


EXACT  SOLUTIONS  TO  SELECTED  FIELD  EQUATIONS 

3 . 1  Introduction 

*j* 

The  existence  of  exact  solutions  to  many  of  the  linear  differ¬ 
ential  equations  of  relevance  to  classical  and  quantum  physics,  has  been 
very  valuable  to  the  development  of  these  disciplines.  The  simple  har¬ 
monic  oscillator  is  perhaps  one  of  the  best  known  examples  of  a  system 
for  which  the  exact  solutions  enjoy  widespread  application.  Such  exact 
solutions  often  serve  as  prototypes  to  illustrate  some  fundamental  phys¬ 
ical  or  mathematical  principle  common  to  a  wider  variety  of  systems.  In 
addition,  exact  solutions  often  find  use  as  the  starting  point  of  an  ap¬ 
proximation  scheme,  used  to  describe  a  more  complicated  system  for  which 
the  solution  is  not  precisely  known. 

The  situation  in  the  case  of  nonlinear  differential  equations, 
however,  may  be  somewhat  different.  As  a  class,  they  appear  to  show 

•jjf 

less  regularity  in  their  properties,  compared  to  linear  equations.  Gen¬ 
eral  theorems  governing  nonlinear  systems  are  less  well  developed.  One 
might  think,  therefore,  that  the  discovery  of  a  nonlinear  system  for 
which  exact  solutions  can  be  found,  would  be  of  little  value  because  it 
might  not  share  many  common  features  with  other  nonlinear  systems.  How¬ 
ever,  it  is  our  view  that  exact  solutions  to  nonlinear  equations  are  im¬ 
portant,  if  only  to  develop  an  appreciation  for  the  variety,  if  not  the 

^  By  exact,  we  mean  expressible  in  closed  form  in  terms  of  well  known 
functions. 

For  a  brief  discussion  of  some  of  the  fundamental  differences  between 
linear  and  nonlinear  differential  equations,  see  Davis  (1962)  . 
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regularity,  displayed  by  nonlinear  systems.  If  a  large  body  of  data  is 
collected  on  the  properties  of  individual  nonlinear  equations,  then  hith 
erto  unsuspected  general  features  may  emerge. 

In  this  chapter  we  present  an  example  of  a  nonlinear  system  of 
equations  for  which  particle-like  solutions  have  been  found  in  the  time- 
independent,  spherically  symmetric  cases  for  which  the  vector  potential, 
A  =  0.  Attempts  are  then  made  to  generalize  these  solutions  to  include 
cases  where  the  vector  potential  is  non-zero.  Calculation  of  the  energy 
and  angular  momentum  of  these  solutions  is  deferred  to  later  chapters. 

We  look  for  solutions  of  the  field  equations  derived  from  the 
Lagrangian  (2.65)  introduced  in  the  previous  chapter.  Only  static  solu¬ 
tions  will  be  sought.  In  the  static  case  the  Lagrangian  becomes 


(3.1) 


* 


\  (Vg)2  +  \  (vxit)  2  -  ±  (V$) 2 


1  ,±. 
2 


+  -jj-  +  K(g)V»it  +  g2 ((j)?— it2)  +  M(g) 


In  our  attempt  to  find  exact  solutions  to  the  fields  equations  derived 

f 

from  the  Lagrangian  (3.1),  we  consider  only  the  cases  6-0  and  6_1*  For 
convenience,  we  also  restrict  the  form  for  the  functions  K(g)  and  M(g) 
to  polynomials  in  g.  Further,  we  assume  that  K(g)  contains  only  powers 
of  g  greater  than  or  equal  to  2,  and  that  M(g)  contains  no  powers  of  g 
less  than  4.  These  conditions  are  useful  in  ensuring  that  the  energy  is 
finite  and  also  have  a  bearing  on  the  allowed  values  for  the  charge. 
The  field  equations  derived  from  the  Lagrangian  (3.1)  are: 


+  1 
'  When  6=1,  the  terms 


—  p 
4  yv 


2  '1  (Ap,m) 


appear  in  the  full 

time-dependent  Lagrangian  (2.65).  This  combination  occurs  in  Fermi's 
treatment  of  the  quantization  of  free  electromagnetic  fields.  See, 
for  example,  Nishijima  (1969) 


•- 
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(3.2.1)  =  2g2£  +  ^K 

(3.2.2)  yZfy  =  _2g2<j) 

(3.2.3)  v2g  =  2(<,242)g  +  ||  $4  +  N (g)  , 

where  N(g)  =  ™  .  From  our  previous  assumptions  about  the  form  for  M(g) , 
N (g)  can  be  written 


(3.3) 


N(g)  =  l  6mgm  ;  (  M(g)  =  l 


3. 


m- 1  m  N 
g  ) 


ml  3 


m>4 


m 


1 


If  the  vector  potential  goes  to  zero  faster  than  — ,  as  r-x>o,  then 
for  N  of  the  form  (3.3),  the  equations  (3.2.2)  and  (3.2.3)  fix  the  charge 
to  the  value  q=l,  for  spherically  symmetric  solutions  with  the  power  ser¬ 
ies  expansions  (2.34).  If  g  goes  as  r“n,  as  r-*»,  then 


(3.4) 


2*i(n-l) 


h 


is  the  value  of  the  fixed  charge.  As  we  will  see  in  a  specific  example, 
there  exist  solutions  for  which  the  vector  potential  goes  as  asymptot¬ 
ically,  in  which  case  the  charge  is  not  fixed. 

Another  feature  of  the  systems  we  study,  is  that  the  Lagrangian 
(3.1)  is  invariant  under  the  transformation  <t>  ~4) .  Therefore,  if  <t>  is 

a  solution  to  the  equations  (3.2),  then  ~4>  is  also  a  solution:  static 
solutions  always  occur  in  oppositely  charged  pairs.  Further,  the  rest 
energy  of  both  solutions  of  such  a  pair,  is  the  same,  as  is  seen  from  the 
expression  (A1.15)  in  Appendix  1.  While  the  existence  of  a  degenerate 
pair  of  solutions,  with  charges  of  both  sign,  is  a  physically  important 
property  of  the  Lagrangian  studied,  there  are  some  related  problems 
which  arise.  One  difficulty  is  that  the  time-dependent  generalization 
of  the  Lagrangian  (3.1),  given  by  (2.65),  does  not  automatically  admit 


A 


. 


■ 


' 

*  .  -  '  )  w 


41 


solutions  with  charges  of  both  signs.  This  problem  arises  because  of  the 

term  K(g)  A  .  Unless  we  assume  that  K(g)  changes  sign  with  the  charge, 

then  this  term  will  not  be  invariant  if  we  replace  A  with  -A  .  But  if 

y  y 

we  assume  that  K(g)  is  proportional  to  the  charge,  then  we  have,  in  ef¬ 
fect,  two  different  Lagrangians.  As  mentioned  in  Chapter  2,  this  would 
cause  some  difficulty  in  discussing  interactions  between  oppositely  char¬ 
ged  particles.  Another  situation  is  one  in  which  K  contains  only  odd 
powers  of  g  and  M  contains  even  powers.  Then  the  Lagrangian  (3.1)  is  in¬ 
variant  under  the  transformation  (A  ,g)  -*  (-A  ,-g).  In  this  special  case 

y  y 

charges  of  both  sign  are  allowed  for  a  single  Lagrangian.  If  K=0,  then 

the  Lagrangian  (3.1)  is  invariant  under  A  -*  -A  ,  so  that  if  (£,<b)  is  a 

y  y 

solution  to  the  associated  field  equations,  then  so  is  (-A, -<}>).  Unfor- 
tunately,  (-A,4>)  and  (A,~<{>)  are  also  solutions.  Therefore,  unless  A=0 , 
solutions  will  occur  in  quadruplets  rather  than  in  pairs.  All  four  solu¬ 
tions,  related  by  the  various  changes  of  sign,  have  the  same  rest  energy. 
In  spite  of  these  properties,  both  the  cases  K=0  and  K^O  will  be  consid¬ 
ered  . 

The  only  solutions  we  have  found,  have  the  field  g  of  the  form 

(3.5)  g(r)  =  g(r)  =  ^  f 
where  c  is  a  constant, 

(3.6)  R  =  (  a2  +  r2  )  2 

and  ’a'  is  another  constant,  which  is  proportional  to  c.  Nash  (1975) 
has  studied  the  equations  (3.2.2)  and  (3.2.3)  in  the  case  ^=0,  K— 0,  with 


•cf  used'  here  has  no  relation  to  the  speed  of  light. 
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N(g)  =  $^g3  +  $^g4  +  8,_gb  .  He  obtained  exact  and  numerical,  spherically 
symmetric  solutions.  In  all  cases  g(r)  was  nodeless,  as  is  the  expres¬ 
sion  (3.5).  Unless  otherwise  specified,  throughout  this  chapter,  we  as¬ 
sume  that  g  has  the  form  (3.5).  Note  that  g  goes  to  zero  as  r~2  as  r-*30, 
so  that  when  the  charge,  q,  is  fixed,  its  value,  from  (3.4)  is  given  by 
q2  =  1.  Also  with  the  form  (3.5)  for  g,  Teshima  (1976)  was  able  to  find 
exact  solutions  in  the  case  6=K=0,  &=0,  with  a  polynomial  expression 
for  N(g).  In  the  following,  we  present  Teshima' s  result  and  attempt  to 
generalize  it  to  the  case  of  non-spherical  solutions  for  which  the  vector 
potential,  it,  is  non-zero. 


3 . 2  The  Case  K=0,  6=0 

We  now  discuss  special  cases  of  the  equations  (3.2),  starting  with 
the  systems  in  which  K  and  6  are  both  zero.  In  this  case,  the  field  equ¬ 
ations,  (3.2)  become 

(3.7.1)  Vx  (Vxit)  =  2g  21 

(3.7.2)  V2<j>  =  -2g2(j) 

(3.7.3)  V2g  =  2  (<j)2-it2)g  +  N  (g)  , 


with  N(g)  given  by  (3.3). 

We  first  deal  with  the  equation  (3.7.2).  Substituting  the  expres¬ 
sion  (3.5)  for  g  into  (3.7.2)  and  using  the  Laplacian  operator  in  spher¬ 
ical  coordinates,  as  given  by  (2.71),  one  has 


(3.8) 


1  9  2  (r<fr) 

r  ~3r2 


<j>  =  0  , 


with  L2  given  by  (2.70)  and  R  defined  by  (3.6). 

This  equation  can  be  solved  by  the  method  of  separation  of  variables. 
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Assume  that  there  exists  a  solution  of  the  form 


(3.9) 


<Kr) 


I  I  V  ^(r)  Ylm<8.*) 

1=0  m=-l  J'  lm 


where  Ti1m  are  complex  constants  and  Y1m  are  the  well  known  spherical 


1m 


harmonics.  Because  the  Y  satisfy  the  relation 


Y  =  (-l)m  Y* 

1m  1  -m 


(  where  denotes  complex  conjugation),  if  we  choose 


(3.10) 


"in.  =  (_1) 


m  * 
*1-1 


■m 


then  the  requirement  that  (j>(r)  be  real  means  that  i|^(r)  is  also  real. 

Using  the  fact  that  the  spherical  harmonics,  Y.  ,  are  eigenfunctions  of 

Xm 

4* 

the  operator  L2,  corresponding  to  the  eigenvalue  1(1+1),  one  finds  that 


the  'radial*  function  \p^  (r)  satisfies  the  equation 


(3.11) 


r  drz  R4  r^ 


’I'l  88  0  > 


for  1-0,  1,  2,  ...  .  To  solve  this  equation,  introduce  the  new  independ¬ 
ent  variable  P,  and  the  new  dependent  variables 


_  a2  a2 

P  =  R2"  =  a^TT2" 


(3.12) 


X1  = 


 R 


21+1 


—  ifi^  ;  1=0, 1,2,  . . . 


In  terms  of  the  new  variables,  (3.11)  becomes 


d2Xl 

(3.13)  p(p-l)^T-  + 


(21+3)p-|(21+3) 


dxl 

dp 


(21+3) (21+1) a2 -2c2 

4a2 


X±  =  0 


See,  for  example,  Merzbacher  (1970),  p  182 


- 


. 


. 


- 
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(3.13)  is  a  special  case  of  the  hypergeometric  equation^: 


(3.14)  p 

where  we  make  the 


(P-D  + 

dp^ 


(y+v+l)p-A 


identifications 


d  y 

^  +  yv  x  =  o 

dp 


9 


(3.15) 


U  =  1  +  1  + 


2c‘ 


v  =  1  +  1  - 


The  equation  (3.14)  has  the  solution 

(3.16)  X  =  P  F(y,v;A;p)  +  q  p1_X  F(y-A+1, v-A+1 ;2-A ;p)  , 


where  F(y,v;A;p)  is  a  hypergeometric  function  and  P  and  Q  are  arbitrary 
constants.  The  term  in  Q  is  not  finite  as  p*K),  (r*») ,  so  we  set  Q=0. 

For  solutions  in  which  X-y-v  ~  0,  as  is  the  case  here,  where  A-y-v=-(l4^) , 
the  term  in  P,  in  general,  will  diverge  at  P=l,  (r=0) .  An  exceptional 
case  is  that  in  which  y  or  v  is  a  non-positive  integer,  in  which  case 
F(y,V;A;p)  reduces  to  a  polynomial  inP,  and  is,  consequently,  finite  in 
the  closed  interval  O^P^l,  (0£r) .  Therefore,  from  (3.15),  the  quantity 


(3.17) 


must  be  a  positive  integer  if  we  are  to  get  a  finite  solution  for  x^* 
Furthermore,  for  a  given  value  of  n,  1  is  restricted  to  the  values 

1=  0,1,2,  ...  ,n-l 


4* 

The  properties  of  the  hypergeometric  function  used  in  this  section, 
are  to  be  found  in  Erdelyi  et.al.  (1953) 

£ 

In  the  context  of  generalized  hypergeometric  functions,  this  is  de¬ 
noted  2F1  (y  ,v  ;A  ;p  ) 


. 


* 
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(3.7.2) 

(3.18) 

(3.19) 


Returning  to  the  original  variables,  one  finds  that  the 
has  the  solutions: 


(j>(r) 


*n(?)  =  Jn 

1=0  m=-l 


g(r)  =  g  (r)  = 

n  R2 


equation 


where 


(3.20) 

(3.21) 
and 
(3.6) 


^nl  "  21+1  Hl+n+1,1  n*+" 1 ;  !+■  2  J  r2  )  » 

R 


n 


R  5  (a^r2)*5 


The  constants  and  the  real  constant  ’a’,  may  depend  on  the  index  n, 
which  can  take  on  the  values  n=l,2,3,  ...  ,  in  the  preceding  expressions. 

For  some  explicit  calculations,  it  is  convenient  to  express  the 
hypergeometric  function  appearing  in  the  solution  for  <p ,  in  terms  of  the 

1*  g 

Gegenbauer  polynomials  ,  C.  The  relation  is 


1 


(3.22) 


F(l+n+l,l-n+l;l+j;pr^ 


)  = 


n+! 

v21+l' 


,1+1 

n-1-1 


■2^2) 


:2+a2 


f  n+1 v  _  (n+1) l 

where  (21+1)  (21+1) ! (n-1-1) !  * 


Next  consider  the  other  field  equations  (3.7.1)  and  (3.7.3). 
(3.7.1)  is  identically  satisfied  by  1=0.  In  this  case  (3.7.3)  reduces  to 

(3.23)  V2g  =  2(f) 2  g  +  N  (g)  . 


Since  we  have  assumed  the  spherically  symmetric  form  (3.19)  for  g,  (3.23) 


^  Erddlyi  et.al.  (1953) 


■i 


1 


■ 


is  satisfied  only  by  a  spherically  symmetric  <f>.  Therefore,  in  the  ex¬ 
pression  (3.18)  for  <f>,  only  the  1=0  term  can  be  included.  The  amplitude 
of  <p  is  determined  by  the  condition  that  the  charge,  q,  is  fixed  by  q2=l, 
as  mentioned  earlier.  Using  the  fact  that  F(p,v;A;p=0)  =  1,  we  find 


(3.24)  ♦(*)  -  *n<?)  =  f  F(n+l,l-n;f;f|)  =  £ 


fr2-a2) 


r2+a2 


and  g(r)  =  gR(r)  is  given  by  (3.19). 


With  <(>  and  g  given,  (3.23)  serves  to  determine 


(3.3) 


N(g)  =  l  B„g 


m^3 


m 


One  finds  that 


(3.25) 


e  (n> 

m 


Vn)  - 


16  (n2-4) 
3  (4n2-l) 


n=l ,2,3, 


•  *  e 


1 


nc 

n 


2  0  m-2  ,  . .  k  _  ,, 

fa  \ m-2  r  (-4)  ,2n+k.  # 

V  ;  i  L k+1  K2k+V  9 
n  k=0 


m=4,5. 


2n  ;  n>  2 


All  other  8  are  zero, 
m 

To  summarize,  the  system  of  equations  (3.7)  is  exactly  satisfied 
by  the  fields  A,  <f> ,  and  g,  if  N(g)  is  given  by  (3.3)  and  (3.25),  with 


A  =  0 

(3.26)  ij>  =  ^  F(n+l,l-n;~;^z-) 

c 

n 

8  R2  * 

C  O  \  i 

o  9  %  4n  -l  H. 

for  n= 1 , 2 , 3 ,  ...  ,  q=±l,  R  =  (az+r^)  ,  and  c  =  a  — - — 

n 

The  solutions  (3.26)  are  those  found  by  Teshima  (1976)  that  were 

previously  mentioned.  An  interesting  feature  of  these  solutions  is  that 

d>  has  (n-1)  radial  nodes  in  the  open  interval  0Sr<°°.  Therefore,  for 
n 

large  n,  the  solutions  have  a  correspondingly  large  number  of  nodes. 


■ 


• 
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This  is  in  contrast  to  the  more  usual  situation,  in  which  the  known  ex¬ 
act  solutions  in  a  nonlinear  field  theory  are  nodeless  or  have  only  a 
small  number  of  nodes. 

To  get  a  better  feeling  for  the  nature  of  these  solutions,  we 
write  out  the  explicit  form  of  <p  and  N  for  the  smallest  allowed  values 
of  n. 

n=l 


n=2 


n=3 


N(g) 

16  3 
=  ~T8 

<f>(r) 

=  £ 

R 

N(g) 

8a4  4 

=  -/8 

;  c 

ii 

o 

N3 

II 

l**a 

4>(r) 

=  f  d- 

2a2  , 

N(g) 

°6o 

II 
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3 

a4  4 
■JTS*  + 

64  a 

3  c1 

<f>(r) 

“f  (l- 

16  a2 

3  R2 

,  16  a4 
3  r4 

)  • 

5  32  a®  c  r35)h 

-  T"T5'gb  J  :  c=c?=hrJ  « 


W8' 


3  c 


A  general  feature  of  these  solutions  is  that,  for  n-3,  N(g)  includes 
3  2n 

terms  from  g  to  g  .  One  can  also  see  the  increase,  with  n,  of  the 
number  of  radial  nodes  in  <|>(r).  The  charge  of  all  the  above  solutions 
is  q=±l. 

Next  we  investigate  the  possibility  of  finding  static  solutions 
in  which  the  vector  potential,  A,  is  non-zero.  When  A  =  0,  the  angular 
momentum  of  any  of  the  time-independent  localized  systems  we  have  studied, 
is  zero.  However,  when  A  is  non-zero  we  may  be  able  to  find  a  system 
representing  a  particle  with  non-zero  intrinsic  angular  momentum.  This 
would  be  an  important  addition  to  our  particle  models. 

The  equation  (3.7.1)  for  the  vector  potential  is  satisfied  by 
any  solution  of  the  system 


■ 
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(3.27)  V2?  =  ~2g2it 

(3.28)  V*£  =  0  , 

although  (3.7.1)  may  have  other  solutions  as  well.  The  reason  for  look¬ 
ing  at  the  pair  of  equations  (3.27)  and  (3.28),  is  that  (3.27)  can  be 
solved  immediately  to  give  the  Cartesian  components,  A^,  of  the  vector 
potential.  This  follows  because  the  structure  of  (3.27)  is  the  same  as 
equation  (3.7.2),  which  we  have  already  solved  for  the  scalar  potential, 
<J>.  Comparing  (3.27)  to  (3.18)  gives 

n-1  1 

(3.29)  A  (r)  -  l  l  (k.)lm  Ylm(e>t)  ^nl(r)  ;  n-1, 2,3 . 

1=0  m=-l 

where  (k.)-  are  the  Cartesian  components  of  the  constant  vectors  lc-  , 
i  1m  lm 

which  satisfy 


(3.30) 


(-1)"  K  m 


y 


and  which  may  depend  on  n. 

By  proper  choice  of  the  constants  the  vector  potential  given 

by  (3.29)  can  be  made  to  satisfy  the  supplementary  condition  (3. 28), and 
thereby  providing  a  solution  to  the  equation  (3.7.1).  The  problem  of 
imposing  the  condition  V*A=0  is  treated  in  Appendix  3.  From  (A3. 9)  we 
have  the  result  that  (3.7.1)  is  satisfied  by 

A  =  A(r,0)  1^  , 

a+n+l .  l-n+1 ; 

x  sin0  F( l+y,y -y;  2 ;  sin20 )  , 

/ 

where  is  a  unit  vector  in  the  direction  of  increasing^  ,  and  are 


n-1  r 

A(r  ,0 )  =  l  k,- 
1=1  R 


' 


i 


I  .B 


we  can 
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arbitrary  real  constants.  Using  the  definition  (3.20)  for  <b  , 

nl 

simplify  this  expression  to 


n-1 

(3.31)  A(r,0)  -  l 

1=1 


K- 


4>ni(r)  sin0  f(1+|»'2  “'J»2;sin20)  ;  1  =  A  L 


The  question  now  arises  whether  the  solution  of  (3.7.1)  for  A, 
given  by  (3.31),  and  the  solution  of  (3.7.2)  for  <J>,  given  by  (3.18), 
also  satisfy  the  remaining  field  equation 


(3.7.3)  V2g  =  2(<j>2-l2)g  +  N(g)  . 


Since  we  have  assumed  that  g=g (r)=^2+^2'»  a  necessary  condition  for 
(3.7.3)  to  be  satisfied  is  that  the  angular  dependence  in  (<j)2-A2)  vanish. 
Other  than  the  case  ^=0,  we  have  been  unable  to  find  any  solution  of  the 
given  forms  (3.31)  and  (3.18),  which  have  this  property.  There  appears 
to  be  at  least  two  reasons  for  this  failure.  One  problem  arises  from 
the  relative  sign  between  the  term  in  t}>2  and  that  in  it2.  The  highest 
value  of  1,  l=n-l,  in  the  expression  (3.18)  for  <j>  gives  rise  to  terms 
in  <f)2  of  the  form  (cos0)  ,  whereas  the  solution  (3.31)  for  A  gives 
terms  of  the  form  sin2 6 (cos0 ) (l-cos20) (cos0)^n  ^ ,  in  it2=A2.  Upon 
calculation  of  4>2-it2  one  finds  that  the  terms  in  (cos0)  reinforce, 


rather  than  cancel,  one  another. 


f 


Another  reason  why  no  solution  could  be  found  is  that  the  condition  V*A-0 
apparently  requires  that  the  ~,  ( 1=0) ,  term  in  expansions  for  A,  should 
vanish,  whereas  for  charged  solutions  the  ~  term  in  <f>  does  not  vanish. 


Because  of  the  recent  interest  in  ’Euclidean'  field  theory  (  see,  for 
example,  Jackiw  (1977)  )  we  also  looked  for  solutions  to  systems  in 
which  the  combination  <f>2+l2,  rather  than  <J>2-A2,  appears.  However,  no 
exact  solution  could  be  found  when  K=0.  A  solution  to  a  system  of 
this  type,  with  K^O,  will  be  given  in  the  next  section. 


' 


•' . 
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This  difference  in  the  form  of  <j>  and  A  appears  to  prevent  cancellation  of 
the  angular  dependence  in  (J)2~A2. 

As  noted  in  Appendix  3,  the  solution  (3.31)  for  the  vector  poten¬ 
tial,  is  not  necessarily  the  most  general  solution  of  (3.7.1).  Other 

possibilities  include  solutions  where  V*A  ^  0  and  in  which  g  is  not  of 

c 

the  simple  form  g^^-Tj.  In  particular,  one  could  look  for  solutions  in 
which  g  included  some  angular  dependence.  Some  variational  calculations 
were  attempted  in  order  to  investigate  the  possibility  of  a  non-spherical 
g.  However,  the  results  were  generally  inconclusive  and  will  not  be  pre¬ 
sented  here.  One  of  the  reasons  for  this  failure  is  that  we  had  no  means 
at  our  disposal  to  judge  the  accuracy  of  the  variational  'solutions'. 

The  Lagrangian  is  extremized  by  the  exact  solutions  to  the  field  equa¬ 
tions,  but  for  an  unrestricted  class  of  variational  trial  functions,  the 
maximum-minimum  character  of  this  extremum  is  unknown. 

3.3  The  Case  6=1 ,  K=Q 

The  restrictions  on  finding  exact  solutions  imposed  by  the  condi¬ 
tion  7-A=0,  introduced  in  the  preceding  section,  can  be  avoided  by  intro¬ 
ducing  a  term  in  (V*A)2  into  the  Lagrangian.  When  5=1  and  K=0  in  the 
Lagrangian  (3.1),  the  field  equations  become 

(3.32.1)  V2X  =  -2g2it 

(3.32.2)  V24>  =  “2g2<f> 

(3.32.3)  V2g  =  2(<})2-it2)g  +  N(g) 

For  large  r,  we  expect  solutions  of  the  equations  (3.32)  for  <f> 
and  it,  to  approach  solutions  of  the  system 


■ 
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l2t  =  0 
V2(j)  -  o  . 


These  are  the  Maxwell  equations,  in  the  static  case,  if  one  imposes  the 
gauge  condition  V»A=0.  However,  there  may  exist  solutions  of  (3.32) 
which  do  not  satisfy  V*A~0.  Such  a  solution  would  not  tend  toward  a  solu¬ 
tion  of  the  Maxwell  equations,  in  the  limit  r-*°°,  that  is,  in  the  region 
where  their  usefulness  has  been  well  established.  However,  because  of 
the  paucity  of  nonlinear  systems  for  which  exact  solutions  are  known, 
it  was  thought  worthwhile  to  investigate  solutions  of  (3.32)  for  which 
V*A^0,  even  though  it  is  expected  that  the  results  may  have  no  phys¬ 
ical  relevance.  These  comments  also  apply  to  the  other  systems  studied 
for  which  6=1  (Section  3.5). 

Equations  (3.32.1)  and  (3.32.2)  have  general  solutions  given  prev¬ 
iously  by  (3.18),  (3.19),  and  (3.29)  as 


n~l 

1 

l 

1=0 

I 

m=-l 

ki 

lm 

Yi* 

♦nl 

n-1 

1 

I 

1=0 

I 

m=-l 

hi 

lm 

Yim 

♦nl 

c 

n 

R2 

y 

where  $  ,  c  ,  and  R  are  defined  by  (3.20),  (3.21)  and  (3.6),  respect 

nl  n 

ively.  If  we  choose  the  constants  and  k  such  that 


(3.33) 


1=0 


=  i  F(n+l,l-n 


3  a2. 
;2;R?) 


then  (3.32.3)  is  satisfied  for  the  same 

section,  namely,  N(g)  =  I  3  g  ,  with  3 

m>3 


choice  of  N(g)  as  in  the  previous 

given  by  (3.25)  . 
m 


K  .  M 


- 
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The  condition  (3.33)  is  satisfied  if 


(3.34) 


(  1 , Y  ,m,m'  not  all  0) 


(TW2 


-  <*oo>2  - 


This  is  easily  verified  by  substituting  the  series  (3.18)  and  (3.29)  into 

_i< 

(3.33)  and  using  the  value  of  the  spherical  harmonic  =  (4tt)  2  in  the 
one  term  that  does  not  cancel  on  the  left-hand  side. 

The  solution  we  have  obtained  to  the  system  (3.32)  suffers  from 
the  fact  that  a  continuum  of  values  exists  for  the  charge,  q.  From  the 
expression  (3.18)  for  <j>,  one  obtains 


(3.35) 


q  (4tt)^  n°° 


Substituting  (3.35)  into  (3.34)  gives 


(3.36) 


q2  =  1  +~k  (*00>2  ’ 


so  that  the  magnitude  of  the  charge  can  take  on  any  value  greater  than 
or  equal  to  1.  The  charge  is  1  only  when  £Q0  =  0,  in  which  case  X=0,  as 
follows  from  the  relations  (3.34). 


3 • 4  The  Case  6=0,  K^Q 

The  next  case  considered  is  6-0,  but  with  no  immediate  restric¬ 
tions  on  K(g) .  The  field  equations  are 

(3.37.1)  Vx(Vxl)  =  2g 2t.  +  VK 

(3.37.2)  V2(f>  =  -2g24> 

V2g  =  2(p~P)g  +  N(g)  +  ~  V-l 


(3.37.3) 


■ 


. 


■ 
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When  K=0  we  found  solutions  to  (3.37.1)  and  (3.37.2),  given  by  (3.18) 
and  (3.31),  which  also  satisfied  (3.37.3)  only  in  the  case  it=0.  It  was 
originally  hoped  that  the  introduction  of  the  non-zero  K  would  lead  to 
solutions  with  a  non-zero  vector  potential.  A,  and  possibly  a  non-zero 
magnetic  field,  it  =  Vxit.  However,  the  only  solutions  to  the  system 
(3.37)  that  we  have  been  able  to  find  with  K^O  are  of  the  form  it=A(r)  "t  , 
where  t-  is  the  radial  unit  vector.  The  curl  of  any  spherically  symmet¬ 
ric  radial  vector  field  is  zero  so  such  a  solution  has  no  magnetic  field. 

With  i£=A(r)  and  ,  equation  (3.37.1)  gives  immediately , 

r  K*- 

the  result 


(3.38) 


1  dg  dK  =  r  dK 
2gz  dr  dg  c  dg 


m 


For  simplicity,  assume  that  K  =  xg  ,  for  some  constant  t  and  where  m  is 
an  integer  greater  than  or  equal  to  2.  From  (3.38)  we  then  find 


(3.39) 

The  second  field  equation, 
by  choosing 


,  xm  m- 1 
A  =  —  r  g 

(3.37.2),  can  be  satisfied  in  the  usual  way 


*  -  1  *o0  =  t  Fd+n.l-n;^ 


(3.40) 


g  =  8n  =  W 


rUn  -l\h  _ 
with  c  =  (—  0  •  1  a 
n  v  2  J 

From  section  3.2,  recall  that  the  functions  g^  and  <j>n0  satisfy 
the  equation 


> 


(3.41) 


v2Sn  =  2®n^n0  +  Vs) 


f  ^  |  '  i  $  $  v 

- 


* 
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n  jjj 

where  NQ(g)  =  2.  ^m8  >  with  $  determined  by  (3.25).  Substituting  (3.41) 

m-3 

into  (3.37.3),  one  finds  that 

(3.42)  (1-q2)  v2g  =  N (g)  -  q2N0(g)  -  2A2g  +  || j.t  . 

All  the  functions,  except  the  unknown  N(g),  appearing  in  (3.42)  can  be 
calculated  in  terms  of  g,  with  the  result 

N(g)  =  q2NQ(g)  +  “(1-q2)  (cg2-4a2g3) 

(3.43) 

-  (  (3-2m)c  +  2ma2  g  )  g2m  2  . 

Previously,  we  have  assumed  that  N(g)  contains  only  powers  of  g 
greater  than  or  equal  to  3.  When  m-3,  this  requirement  on  the  form  of 
N(g)  is  equivalent  to  fixing  the  charge  to  the  value  q  “1,  because  other¬ 
wise  a  term  in  g2  would  appear  in  the  expression  (3.43)  for  N(g).  With 
q2~l ,  we  find  from  (3.43)  that 

(3.44)  N(g)  =  N0(g)  -  f  (3-2m)c+2ma2g  )  g2”'2  ;  m-3, 4,5,  ... 

The  case  m=2  is  somewhat  different.  One  finds  from  (3.43)  in  this  case: 

N (g)  =  q2NQ(g)  +  (l-q2+2T2)  (cg2-4a2g3) 

Therefore,  to  eliminate  the  g2term  from  N(g),  we  need 

(3.45)  q2  =  l+2x2. 

With  the  charge  given  by  (3.45),  the  expression  for  N(g)  reduces  to 

(3.46)  N(g)  =  (1+2t2)  N0(g)  ;  m=2  . 

For  a  given  value  of  the  parameter  t  appearing  in  K—rg the  charge  will 
be  fixed  at  the  value  given  by  (3.45),  when  m=2,  so  q2=l  only  for  x=0. 
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in  this  case,  but  for  m^3,  q2=l  for  all  x. 

The  particular  form  K~xgra  is  not  necessary  to  find  solutions.  By 

an  appropriate  choice  of  K  and  N  one  can  find  a  wide  variety  of  solutions 

c 

for  which  =  A(r)  t  ,  <j>  =  q<b  A,  and  g  =  — y  . 

r  nu  R^- 

As  noted  earlier,  the  current  interest  in  ’Euclidean’  field  theory 

o  ”>2 

has  led  us  to  examine  systems  in  which  the  combination  <j>  +A  ,  rather  than 
appears  in  the  Lagrangian.  Specifically,  choose  m=2,  so  K=xg2. 

The  Lagrangian  studied  was 


4tt 


j  (Vg)2  -  ”  W)2-  |  (VxA)2  +  RV-A  +  (A2+<|>2)g2  +  M(g) 


This  leads  to  the  field  equations 


Vx(Vx^)  =  2g2X  -  2xgVg 


(3.47) 


V 2  4>  =  -2g2<j> 

v2g  =  2(t^«M2+<J>2)  +  N(g)  . 


For  appropriate  N(g),  solutions  were  found  of  the  form 

c  f4n2-l^  2 
n  _  ( — « — J  a 

8  =  “  v-  2  ’  - 
K  R2 

2 

<t>  =  F(2+n,2-n;4;f2-)  cosB  ,  where  k  is  a  constant, 

R  ^  R 

A  =  Ar(r)  lr  +  A^(r,B)  1^  , 

2xr  1 

where  A  =  »  with  x  -  ±  — , 

r  R  2^ 

ler  5  a 2 

and  A^,  =  ±  <J>  cotB  =  ±  F(2+n, 2-n;— ;^y)  sin0 


,  cos6  ,  t  2x  ■>  .  ,  sin6  + 
In  the  limit  r-*»,  we  find  <j>  +  k  and  A  -  r  lr  ±  k  r7 


' 
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The  solution  is  seen  to  be  neutral,  and  describes  a  particle  with  equal 
electric  and  magnetic  dipole  moments  of  magnitude  k.  (  The  radial  term 
in  A  does  not  contribute  to  the  magnetic  field  and  therefore  does  not 
affect  the  dipole  character  of  the  solution.) 

No  other  solutions  to  the  system  (3.47)  were  found. 


3 • 5  The  Case  6=1,  K^Q 

The  least  satisfactory  results  were  obtained  in  the  case  where 
both  the  parameter,  5,  and  the  function, K(g) ,  are  non-zero.  With  6=1 
the  field  equations  become: 

(3.48.1)  V2it  =  -2g 21  ~  VK 

(3.48.2)  V2<j)  =  — 2g2<J> 

(3.48.3)  V2g  =  2(<j>2-l2)g  +  V-X  +  N(g)  . 


Again,  throughout  this  section  we  assume 


(3.49) 


c  /4n2-l'|  2 

g  =  -rT  *  where  c  =  cn  =  [—  j  a 

n-1  1 

*  =  I  I  nlm  Ylm  ^nl  * 
1=0  m=-l 


and 


with  ^  given  by  (3.20). 

When  K=0,  we  found  in  (3.29)  that  (3.98.1)  has  the  solution 


(3.50) 


n-1 

i  =  i 


1 

l 


1=0  m=-l 


Here  we  have  used  a  superscript  *h’  to  denote  that  (3.50)  is  a  solution 
of  the  homogeneous  equation  ,  (K=0)  .  To  find  the  general  (regulai ) 


■ 


*s 
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solution  of  the  inhomogeneous  equation  (3.A8.1)  we  only  need  to  find  any 
finite  solution,  t?  9  which  satisfies  (3.48.1)  and  add  it  to  the  homogen- 


eous  solution 

(3.50).  We  look  for  a  particular  solution  of  the  form 

(3.51) 

N 

'w' 

SI 

% 

Substituting  (3.51)  into  (3.48.1),  one  finds  that  AP(r)  satisfies  the 


equation 

(3.52) 

1  4%^  +  2(^-4)  AP  =  . 

r  drz  '-R4  r2;  dr 

i  P 

The  method  of  variation  of  parameters  can  be  used  to  find  A  be¬ 
cause  the  general  solution  is  known  to  the  corresponding  homogeneous  eq¬ 
uation 


(3.53) 

1  i£(rAP)  +  1  ■ )  Ap  =  0 

r  d7^  +  V  7^3  A  0  ’ 

The  equation  (3.53)  is  a  special  case,  (1=1),  of  the  equation  (3.11). 

From  (3.16)  we  find  that  the  general  solution  of  equation  (3.53)  is 

q  2  2 

(3.54)  AP  =  S  “3  F(n+2,2-n;-2  ;^-)  +  T  r F (n+^ , *S-n ; ~h ;|j) 

In  the  method  of  variation  of  parameters  we  look  for  a  solution  in  which 
S  and  T  are  functions  of  r.  Substituting  (3.54)  into  (3.52),  one  finds 
that  the  latter  equation  is  satisfied  if 


(3.55) 

dS  dK  ^2  . 

d7=d7T  and 

(3.56) 

dT  _  dK  ^1 
dr  “dr  W 

where  ij»  ,  ip2,  and  W  are  given  by 


Ross  (1964),  p  118. 


■  5  s 


. 

. 


. 
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bJ  F(n+2’2  n;2’R2^  "  ^nl 

(3*57)  ^2  =  r  F(n+is>^“n;~is;^2) 

d\p 

w  =  *idF  ~  *2dT  * 

W  has  the  simple  form  ,  -^r,  where  id  is  a  constant  determined  from  the 
functions  4^  and  . 

No  attempt  was  made  to  integrate  the  equations  (3.55)  and  (3.56) 
for  all  values  of  n=2,3,  ...  ,  but  the  case  n=2  was  sufficiently  simple 
that  the  integrations  could  be  done  for  K  of  the  form  K=Tgm,  with  m^2. 
When  n=2 ,  one  finds  from  (3.57),  that 


leading  to 

(3.58) 

(3.59) 


2mTcm  f  (r7+9a2r5-9a**r3-a6r) 


S  (r)  =- 


R 


2mf5 


dr 


T(r)  = 


2mxc 


m 


2m+5 


dr 


R 


These  integrals  can  be  done  in  terms  of  elementary  functions  for  the 
cases  where  m  is  a  positive  integer  or  half  a  positive  integer,  but  the 
general  expressions  will  not  be  used  here. 

The  particular  solution  AP,  from  (3.54)  and  (3.57),  is  given  by 


(3.60) 


AP  =  S(r)  +  T(r)  4^  . 


4. 

The  form  of  W,  up  to  a  multiplicative  constant,  is  determined  by  the 
differential  equation  (3.52)  for  4^  and  4>2*  See,  for  example, 

Ross  (1964) ,  p  376. 


■ 


■ 
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One  can  show  that  T(r)  and  go  to  zero  as  r*»,  for  K  of  the  form  K=xgm, 

with  m^2.  However,  from  (3.57),  as  r-*»,  so  that  T  must  vanish  as 

1  n 

fast  as  ~2"»  in  this  limit,  if  A1  is  to  vanish  asymptotically.  From  (3.59) 
we  find 


(3.61) 


lim  _  s  v  2nrr  c 
r-Ko  T(r>  =  — T 


m 


oo 


0 


R 


2mf5 


dr 


The  choice  of  the  lower  limit  of  integration  just  corresponds  to  a  part¬ 
icular  integration  constant.  This  choice  ensures  that  T(r)  =  0, 
which  is  necessary  for  AP  to  be  finite  as  r->0.  From  (3.61)  we  find 


(3.62) 


lim  TCm  rCy)  r(m+l) 

£  T(r>  ■  7^  ,5. 


3  a 


r  (m+2) 


lim 


where  V  is  the  well  known  special  function.  (3.62)  shows  that  *7“  T  does 

^*->00 

not  vanish  for  any  positive  m,  so  that  the  resulting  A*5  will  not  be  finite. 


One  finds,  for  example,  when  m=2,  that  A  ,  as  given  by  (3.60),  is 


AP(r)  -  %  r 


These  considerations  show  that  no  acceptable  solution  can  be  found  in 
which  K  has  the  simple  form,  K=xgm.  However,  if 


(3.63) 


K  =  xgm  +  pgP  ;  with  p>m*2. 


the  extra  parameter,  n,  can  always  be  chosen  such  that  ^  )co  T(r)  -  0,  and 
an  AP(r)  can  be  found  which  is  finite  everywhere.  For  example,  if  p=m+l, 
then 


(3.64) 


a2  (2m+5) 
c  (2m+2)  T 


is  the  appropriate  value  for  p.  The  general  expression  (3.60)  for  A 
in  such  a  case,  would,  however,  be  cjuite  complicated  and  the  explicit 


solution  was  not  found. 


■ 


, 


The  preceding  discussion  has  concerned  only  the  field  equation 
(3.48.1).  Yet  to  be  treated  is  the  equation  (3.48.3).  Only  in  the  case 
n=2  was  any  reasonable  progress  made  in  analyzing  this  equation.  Substi¬ 
tuting  <J>  and  g  from  (3.49),  and  using  A  =  Ah  +  AP  t  ,  with  Ah  of  the  fora 
(3.50),  it  was  found  from  (3.48.3)  that  the  angular  dependent  terms  could 
be  cancelled  for  a  very  specific  form  of  K(g).  This  K  is  given  by  (3.63) 
and  (3.64)  with  m=|,  that  is, 

5 

(3-65)  K(g)  =  g2  (  1-^-  g  )  • 

If  K  was  not  of  the  special  fora  (3.65),  no  solutions  could  be  found,  al¬ 
though  the  possibility  that  some  other  expression  for  K  could  lead  to  sol¬ 
utions,  was  not  entirely  removed. 

Because  of  the  appearance  of  the  fractional  powers  in  the  express- 
sion  (3.65)  for  K  and  the  fact  that  the  analysis  was  restricted  to  the 
case  n=2 ,  the  study  of  the  case  where  6  and  K  are  both  non-zero  was  felt 
to  be  of  less  general  interest  than  that  of  the  preceding  cases.  In  the 
subsequent  analysis  of  the  energy  and  the  angular  momentum  of  the  solu¬ 
tions,  the  case  6=1  and  K^O,  will  be  excluded. 


■■ 
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CHAPTER  4 


ENERGY  OF  EXACT  SOLUTIONS 


For  the  time  -  independent  systems  studied  in  Chapter  3  we  can  cal¬ 
culate  the  energy  corresponding  to  the  exact  solutions  that  have  been 
found.  The  energy,  so  calculated,  is  the  rest  energy  or  rest  mass  of  the 
particle  described  by  the  solution.  As  indicated  in  Appendix  1,  section 
A1.3,  the  rest  energy  is  given  by 


E  = 


where  is  given  by  (3.1).  Alternate  expressions  for  E  can  be  obtained 


by  using  integral  relations,  satisfied  by  solutions  of  the  field  equa¬ 
tions  obtained  by  considering  scale  and  amplitude  variations  of  the  stat¬ 
ic  Lagrangian.  For  some  systems,  such  as  the  example  of  case  IB,  (see 
equation  (2.51)  ),  the  sign  of  the  energy  can  be  determined  from  these 
integral  expressions  for  E,  without  knowing  the  form  of  the  exact  solu¬ 
tions.  However,  in  the  examples  considered  in  Chapter  3,  this  was  not 
possible,  so  the  integrals  have  to  be  carried  out  using  the  explicit 
form  of  the  solutions. 


Our  calculations  of  the  energy  will  be  based  on  the  expression 


(  (  Vg)2  +  2(P-^)g2  )  d3x  ; 


(A1.15) 


which  was  derived  in  Appendix  1.  Note  that  we  do  not  need  the  explicit 
form  for  K(g)  or  M(g)  to  evaluate  the  integral  (A1.15). 

The  first  case  to  be  considered  is  6  =  K  =  0.  The  solution  for 
which  the  energy  is  to  be  evaluated,  in  this  case  is  ,  from  (3.26)  and 


(3.22): 
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A  =  0 


(4.1) 


*  =  q,f,nO(r)  =  Cn-1 


rr2-a21 


,2+a2 


g  =  g 


n 


n  ”  R2 


where  q2=l,  R=  (a^fr2)  2,  and  c  = 

n 


f4n2-l^ 


a  ;  n= 1 , 2 , 3 ,  ... 


Denoting  the  energy  of  the  system  by  E^Cn),  (A1.15)  gives 


(4.2) 


E0{n)  12tt 


(K>2  - 


2<{)2.g2 
yn0  n 


d3x 


Transforming  to  spherical  coordinates  and  doing  the  angular  integration 
gives 

rco  r  dg  2 

(4.3)  EQ(n)  =  ±  f~")  -  2d>2.  R2  r2dr 


(4.4) 


(4.5) 


l 

3 


0 


2  2 
3  n 


* 

f2r2  1 

fcl 

fr2-a2^ 

2  ' 

[w  ^ 

[Cn-1 

r2+a2 

> 

* 

r2  dr  ,  using  (4.1) 


0 


|c  (  21  -  "2 1 9  )»  where 

3  n  I  n^  L 


\  5 


rH  ,  _  _ 7T 

R8  r  32a2 1  a 


,  and 


0 


f  r2 

r  i 

r* 

fr2-a2^ 

t 

R6 

L  - 

n-1 

r2+a2 

4 

dr 


0 


1 

a-u2)^ 

r  i  ) 

r1  ('u') 

- 1  a  | 

U  u  ; 

n-1 

k  / 

du 


-1 

where  u  = 


r2-a2 
r2+a2  * 


integral  in  (4.5)  is  just  the  normalization  integral  for  the 
Gegenbauer  polynomials  ,  having  the  value  (Erdelyi  et.  al.  (1953),  p  177) 


TT 


I2  16a^ I  a 


(4.6) 


' 


. 

■ 
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Substituting  (4.4)  and  (4.6)  for  1^  and  I^,  respectively,  into  the 
expression  (4.3)  for  the  energy,  one  has 

c2 

E0(n)  =  (1~^>  > 


or  using  the  definition  of  c  , 

n 


(4.7) 


IT 


Vn)  =  4^ 


(4n2-l) (n2-l) 


n-1 


,  n  1,2,3,  ... 


The  rest  energy  is  seen  to  be  positive  for  n=2,3,  ...  ,  but  is  0 
for  n=l.  Because  the  parameter,  a,  can  depend  on  n,  the  spectrum  of 
rest  energies  is  not  restricted  to  the  explicit  dependence  on  n  given  in 

(4.7) ,  but  can  take  on  any  form,  (except  for  n=l,  for  which  is  always 

0) .  'a'  characterizes  the  size  of  the  particle,  in  the  sense  that  the 

value  of  g  at  r=a  is  one  half  the  maximum  value,  which  occurs  at  r=0. 
Therefore,  to  say  that  'a'  depends  on  n  is  equivalent  to  saying  that  the 
particles  corresponding  to  different  n  may  be  different  sizes. 

The  next  case  to  consider  is  6=1,  K=0,  for  which  we  found  solutions 
given  by  (3.18),  (3.29),  and  (3.19)  as 

1  *taTl»(e*W  ♦nl<r) 

1=0  m=-l 

(4.8) 

1=0  m=-l 
c 

8  =  gn  =  W 


The  parameters  and  are  chosen  such  that 


(3.30) 


<J)2  —  ^2  _  (j. 


2 

n0 


Therefore,  the  energy,  as  given  by  (A1.15),  is 


. 
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(4.9) 


E(n) 


1 

12  7T 


1 

1 2  7T 


(^g)  2  +  2(l2-<J>2)g2 


d  he 


(\)  2  - 


d  he 


=  E0(n)  »  from  (4.2), 


where  Eg(n)  is  given  by  (4.7). 

Therefore,  the  rest  energy  spectrum  for  the  solutions  (4.8)  is  the 
same  as  that  of  the  spherically  symmetric  solutions  (4.1).  To  each  dis¬ 
crete  value  of  the  energy,  however,  there  is  an  infinite  number  of  solu¬ 
tions,  corresponding  to  various  values  of  the  constants  and 
Even  if  one  were  to  restrict  solutions  to  a  specific  value  of  the  charge, 
q  =  (4tt )  2  n  ,  there  are  still  an  infinite  number  of  solutions  for  each 
of  the  discrete  energies. 

The  final  case,  for  which  we  can  calculate  the  energy,  is  6-0, 
K=xgm.  The  solutions  here  are  given  by  (3.39)  and  (3.40)  as 


t  mx  ,  Nm-1  f 
A  =  —  r  (g  )  1 

c  n  r 

n 


(4.10) 


<t>  =  q  <t> 


n0 


n 


8  =  8n  =  R? 


Again,  starting  from  the  expression 


(A1.15) 


E  = 


12tt 


(Vg)2  +  2(A2-<f>2)g2 


d3x 


we  find,  substituting  in  the  solutions  (4.10) 


E  = 


(4.11) 


12  7T 


„  2x2  2  ,  2m2x2  2  2m  ,  3 

(Vgn)2-2q2^0g2  +  -^-r  gn  dx 

n  m 
00 


=  (l-S)  +  |™2x2(c2)m-1 


n 


n 


0 


R 


4m 


dr 


•  •  i 


*  M 

. 


where  the  previous  result  (4.7)  has  been  used  in  evaluating  the  first  two 
integrals.  The  remaining  integral  is 


65 


(4.12) 


d_  =  l!  1  (4m-7) ! 

^4m  2  ^  i4m-5  (2m-4) ! (2m-l) ! 


Substituting  the  expression  (4.12)  into  (4.11)  gives 

2  2m- 2 

f f.  1  p  _  ir(4n2-l)  (  q2.  771  cn  m2(4m-7)! 

(4-13)  E  -  (1‘^  +  7774^5  T2M-4)  !  (2m-l7T 


Using  the  definition  (4.7)  for  EQ(n)  we  find 


(4.14) 


E  =  EQ(n) 

■  Vn) 


(1+2t2)  ,  when  m=2,  and 


1  3x2n2m2  (4m-7)! 

f  4n2-f 

m-2  ' 

(n2-l)  (2m-4) ! (2m-l) I 

32i?1 

4 

for  m  £  3. 


The  difference  in  the  form  for  m=2,  and  m^3,  follows  from  the  difference 
in  the  charge  in  the  two  cases.  From  (3.45),  the  charge,  q,  is  given  by 
q2=l+2x2  when  m=2,  whereas  q2=l  for  m^3. 

Again  in  this  example  we  find  that  the  rest  energy  is  positive. 
The  last  case  considered  has  the  property  that,  if  the  parameter  x,  ap¬ 
pearing  in  K(g),  is  unrestricted,  there  is  a  continuum  of  allowed  rest 
+ 

energies  .  This  undesirable  feature  can  be  eliminated  by  giving  x  some 
fixed  numerical  value.  For  any  given  x  ,  there  is  a  discrete  spectrum 
of  energies  corresponding  to  different  values  for  n. 

Note  that  if  x  is  different  from  zero,  then  the  nodeless  'ground 
state',  n=l,  no  longer  has  zero  energy,  if  m-3.  In  particular,  from 
(4.13): 


E(n=l,m-3)  = 


3tt  x2m2(4m-7) .' 
16  (2m-4) ! (2m-l)  ! 


32a' 


m-2 


t 


Also,  when  m=2,  there  is  a  continuum  of  charge  values. 


We  see  that  the  non-zero  vector  potential.  A,  resulting  from  the  term 
K  V*A,  in  the  Lagrangian,  can  contribute  to  the  energy,  even  though  it 
does  not  produce  an  electric  or  magnetic  field.  This  is  a  particular 
example  of  how  the  potential  in  a  non-gauge- invariant  theory  can  take  on 
physical  significance. 

In  the  examples  considered  in  this  chapter,  we  have  seen  that, 
even  in  a  classical  theory,  it  is  possible  to  have  the  energy  restricted 
to  a  spectrum  of  discrete  values.  Also,  we  have  found  that  the  rest  en¬ 
ergies  of  all  the  particle-like  solutions  obtained,  are  positive.  In  the 
next  chapter,  we  consider  another  fundamental  property  of  elementary  part¬ 
icles,  namely,  the  angular  momentum. 


CHAPTER  5 


ANGULAR  MOMENTUM  IN  CLASSICAL  FIELD  THEORY 


5. 1  Total  Angular  Momentum 


•j- 

In  the  treatment  of  classical  field  theory,  many  authors  discuss 
the  conservation  of  angular  momentum.  However,  there  seems  to  be  only  a 
minimal  discussion  of  the  circumstances  under  which  it  is  useful  to  div¬ 
ide  the  total  angular  momentum  into  orbital  and  intrinsic,  (spin),  con¬ 
tributions.  The  main  topic  of  this  chapter  is  this  division  of  the  total 
angular  momentum. 

* 

It  follows  from  Noether's  theorem  that  the  total  angular  momentum 
of  a  system  is  conserved  if  the  Lagrangian  is  invariant  under  infinites¬ 
imal,  homogeneous  Lorentz  transformations  and  does  not  depend  explicitly 
on  the  coordinates,  x  .  The  conserved  total  angular  momentum  has  spatial 

u 

components,  J^,  given  by 


(5.1) 


* 


where  a  is  a  space-like  hypersurface  with  a  normal  specified  by  da^. 

The  angular  momentum  density,  J  to  be  defined  shortly,  satisfies  the 
continuity  equation 


(5.2) 


J 


yvA ,  A 


=  0  . 


The  relation  (5.2)  ensures  that  the  definition  (5.1)  is  independent  of 


These  include  Aharoni  (1965),  Barut  (1964),  Rzewuski  (1967)  and  Soper 
(1976). 

See,  for  example,  Soper  (1976),  p  105. 
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the  choice  of  the  surface  a. 
one  has 

(5.3)  Jk  = 


Taking  a  to  be  a  surface  of  constant  time, 


1 

2  £klm 


The  definition  of  the  angular  momentum  density,  J  , ,  depends  on 

yvA 

the  Lorentz  transformation  properties  of  the  fields  involved  in  the  sys¬ 


tem  of  study.  Define  an  infinitesimal  homogeneous  Lorentz  transformation 


by 


x  x'  =  x  +  6x  ,  with 

y  y  y  y 

(5.4)  6x  =  e  x  ; 

y  yv  v 

c  =  -e 
yv  vy 


The  four-vector  potential  A  transforms  according  to 

y 


(5.5) 

(5.6) 


A  (x  )  ->  A"(x")  =  A  (x  )  +  ^  e 
y  A  y  A  y  A  l 


op 


(S  )  =  -i  (  6  6 

op  yv  ay  pv 


(S  )  A  (x  ) 
ap  yv  v  A 

"66  )  . 
av  py 


,  where 


For  the  systems  we  have  studied,  in  which  the  Lagrangian, ,  is  con¬ 
structed  from  a  scalar  field,  g,  and  the  vector  field  A^,  the  angular 
momentum  density  is  then  defined  as 


(5.7) 

(5.8) 

(5.9) 


Here  T 

yv 


J  =  x  T 
yvA 

S 


:  -  x  T  +  S 

vA  v  yA 


I  ,  =  i  (S  )  A 
yvA  yv  ap  a 


=  A 


v^yA  '  yvA 

*c 

3.C 


y  3  a. 


-  A 


v  ,A 


v  3A 


y  »A 


is  the  canonical  energy-momentum  tensor. 


,  where 


,  using  (5.6) 
given  by 


T 

yv 


<5  ~  A 

yv0^  o  ,y 


8A 


“  g 


9tC 


a  ,v 


>y  3g 


(5.10) 
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For  problems  involving  other  tensor  or  spinor  fields  than  A  ,  the  trans- 

M 

formation  law  analogous  to  (5.5)  is  used  to  define  appropriate  general¬ 
izations  of  the  matrices  S  ,  which  introduces  additional  terms  into  the 

0p 

definition  (5.7)  of  the  angular  momentum  density. 

While  the  definition  (5.7)  is  the  one  which  arises  directly  from 
the  invariance  of  the  Lagrangian  under  homogeneous  Lorentz  transforma¬ 
tions,  there  is  another  common  definition  of  angular  momentum,  which  we 
0 

denote  by  J  This  second  definition  is  based  on  the  symmetric  energy- 

momentum  tensor  0  ,  introduced  by  Belinfante  (1940)  and  Rosenfeld  (1940), 

y  v 

which  exists  for  any  Lorentz  invariant  Lagrangian.  0  is  related  to  the 

yv 

canonical  energy-momentum  tensor,  T^,  by 


(5.11) 

0  =0  E  T  +  g  ,  . 

yv  vy  yv  yvA,A 

(5.12) 

g  E  -%  (  S  +  S 

yvA  v  yvA  Avy 

S  ,  being  the  tensor  introduced  in  (5.8). 
yvA 

(5.13) 

S  =  -  S  , 

yvA  vyA 

leads  to  the  result 

(5.14) 

gyvA  SyAv 

,  where 


which,  in  turn,  ensures  that 


(5.15) 


0  =  T  +  g  =  0  ,  (using  T  =  0)  . 

y  v , v  yv, v  'yvA,Av  yv,v 


Therefore,  the  energy-momentum 
(5.16)  P®  = 


0  da 
yv  v 


is  conserved.  Furthermore,  since  0  and  T  differ  by  a  divergence, 

r  v  ' 


'«  • 


■ 


■ 
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(5.17) 


P6  =  P  = 

y  y 


T  da 
yv  v 


for  fields  which  fall  to  zero  sufficiently  fast  at  spatial  infinity. 

Now  define  the  angular  momentum  density,  J°  ,  by 

yvA 


(5.18) 


J0  =  x  0  -  x  0  . 

yvA  y  vA  v  yA 


0 

and  the  corresponding  total  angular  momentum,  J^,  by 


(5.19) 


9  -  I  J.  .  da. 


^k  2eklm 


1mA  A 


'klm 


Xl9mX  S 


From  the  result  (5.15)  and  the  symmetry  of  0  one  can  easily  show  that 


(5.20) 


JyvA,A  0 


0 

which  ensures  that  the  definition  (5.19)  for  is  independent  of  the 
surface  a.  Choosing  a  to  be  a  surface  of  constant  time,  we  have 


(5.21) 


J  =  -i  e 
k  klm 


x.0  .  d3x 
1  m4 


0 

We  now  compare  the  total  angular  momentum,  ,  defined  by  (5.21), 
to  the  first  definition,  given  by  (5.3)  and  (5.7).  Using  the  defin¬ 
ition  (5.11)  of  the  symmetrized  tensor  0^v,  one  finds  from  (5.21): 


(5.22) 


T0  _  .  _ 

J.  “i£.  i 

k  klm  J 


x- (  T  ,+g  .  )  d3x 

1  m4  m4A , A 


and  from  (5.3)  and  (5.7): 


(5.23) 


Jk  "  1£klm 


Xl<  Tm4H42Sl»4  5  d3x 


Combine  (5.22)  and  (5.23)  to  give 


(5.24) 


Jk  “  Jk  "1£klm 


(  x.g  )  d3x 

1  m4A , A  lm4 


i  . 


1 
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Using  gm^  0,  (which  follows  from  (5.14)),  (5.24)  can  be  written 


0 

J,  -  =  -i  e.  _ 

k  k  klm 


(xlgm4n\n  "  8m41  **  Slm4 


d3 


x 


The  last  two  terms  in  this  expression  cancel  when  one  substitutes  the 
definition  of  g  from  (5.12).  This  leaves 


(5.25) 


J?  -  J  =  -i  e  [  (x  g  ,  )  d3x 
k  k  klm  J  l&ra4n  ,n 


This  last  integral  can  be  converted  to  a  surface  integral  at  spatial  in¬ 
finity.  Therefore,  if  the  fields  fall  off  fast  enough  as  r-*»,  then  the 

6 

right-hand  side  of  (5.25)  will  vanish  and  we  find  =  J^.  However, 
there  may  be  cases  in  which  the  fields  do  not  go  to  zero  fast  enough  so 
one  should  always  be  aware  that  there  is  the  possibility  of  a  difference 
in  the  two  definitions  of  the  total  angular  momentum.  In  the  discussion 
of  the  division  of  the  total  angular  momentum  into  spin  and  orbital  parts, 
we  find  it  convenient  to  refer  to  the  total  angular  momentum  rather 
than  When  the  two  definitions  may  not  be  equivalent,  however,  both 

K 

forms  are  used. 


5 • 2  Orbital  and  Spin  Angular  Momentum 


The  total  angular  momentum  density,  J  is  given,  from  (5.7), 
as 

(5*7)  JpvA  “  xyTvA  “  xvTyA  +  SyvA 

Looking  at  the  terms  that  appear  in  (5.7),  we  define 

(5-26)  SvA  '  VvA  '  VuA  and 

(5.27)  Lfc  =  h  eklm  |  LlmA  =  eklm  j  x1TnL\  %  • 

a  a 


■ 


. 


yv 


* 
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Also  define 


(5.28) 


eklm 


a 


a 


where  (5.9)  has  been  used  to  obtain  the  second  equality  in  (5.28). 

Recall  that  the  momentum  of  the  system  with  canonical  energy-mom- 


Therefore,  the  quantity  L 


defined  by  (5.26),  has  a  form  similar  to 


jky 


the  density  of  the  angular  momentum,  rxp,  of  classical  particle  dynamics. 
This  similarity  has  led  some  authors  to  call  L^,  as  defined  by  (5.27), 
the  orbital  angular  momentum  of  the  system.  The  term,  S^,  given  by  (5.28), 
is  then  identified  as  the  intrinsic  angular  momentum,  or  spin,  of  the  sys¬ 


tem.  The  total  angular  momentum,  J^,  is 


J,  =  L,  +  S. 
k  k  k 


Soper  (1976),  p  114,  mentions  that  the  volume  integral  of  S  ^  is  the 
operator  for  the  spin  in  quantum  field  theory,  lending  support  to  the 
hypothesis  that  is  the  spin. 

However,  there  are  several  problems  associated  with  these  identi¬ 
fications  for  the  orbital  and  spin  angular  momenta.  In  general,  the  def¬ 
initions  (5.27)  and  (5.28),  of  Lk  and  Sfc,  are  not  independent  of  the 
spacelike  surface  o.  So,  for  example,  if  o  is  a  surface  of  constant  time, 
t,  for  one  observer,  a  Lorentz- transformed  observer  cannot  use  his  surface 
of  constant  time,  t",  and  expect  to  calculate  values  for  Lk  and  Sk  which 


are  consistent  with  those  calculated  by  the  first  observer.  One  has  the 
problem  of  deciding  which  hypersurface,  o,  to  use  in  the  evaluation  o! 


4 


■ 


. 

. 
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the  integrals  for  and  S^, 


The  condition  for  and  to  be  independent  of  q  is 


t 


(5.29) 

(5.30) 


L  =  0  ,  or 

1mA ,  A 

S.  .  .  =  0  . 

1mA ,  A 


(5.29)  and  (5.30)  are  also  just  the  conditions  necessary  for  the  orbital 
and  spin  angular  momenta,  and  S^,  to  be  separately  conserved.  It  is 
easy  to  show  that  (5.29)  and  (5.30)  are  satisfied  if  the  canonical  energy- 
momentum  tensor  is  symmetric,  for  in  that  case  we  have,  using  the  defini¬ 
tion  (5.21)  of  L  .  and  the  result  T  =0  : 

yvA  yv,v 


(5.31) 


L  =  (  x  T  -  x  T  )  =T  -T  =0  . 

yvA,A  y  vA  v  yA  ,A  vy  yv 


Therefore,  when  the  canonical  tensor  is  symmetric  the  orbital  and  spin 
angular  momenta  (5.27)  and  (5.28)  are  well  defined  and  conserved.  This 
result  has  led  Hilgevoord  and  Wouthuysen  (1963)  and  Hilgevoord  and  de  Kerf 
(1965),  to  propose  new  definitions  for  the  spin  and  orbital  angular  mom¬ 
enta,  so  as  to  be  conserved  for  free  field  systems,  even  in  the  cases 

* 

where  the  canonical  tensor  of  the  usual  formulation  is  not  symmetric. 
However,  as  of  the  1965  paper,  the  method  had  been  applied  only  to  those 
systems  for  which  the  field  equations  could  be  expressed  as  Klein-Gordon 
equations  subject  to  some  subsidiary  conditions.  (The  Dirac  equation  is 
of  this  type.) 


f 


* 


(5.29)  and  (5.30)  are  not  independent  because  J  ^  =  +  anc* 

JyvX,A  =  °' 

More  correctly,  the  form  of  the  definitions  (5.27)  and  (5.28)  is  re¬ 
tained,  but  a  Lagrangian  giving  rise  to  a  non-symmetric  energy-momen¬ 
tum  tensor  is  replaced  by  a  new  Lagrangian,  plus  subsidiary  conditions 
on  the  field  equations,  for  which  the  canonical  energy-momentum  tensor 
is  symmetric. 


* 


y> 


\ 


- 
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We  now  examine  some  other  aspects  of  the  spin  and  orbital  angular 
momentum.  Suppose  that  one  has  chosen  the  surface,  a,  to  be  a  surface  of 
constant  time  for  some  observer,  so  that  and  are  unambiguously  de¬ 
fined.  Then  from  (5.27)  and  (5.28)  we  have 


(5.32) 


Lk  1  eklm 


x.,  T  .  d3x  ,  and 
1  m4 


(5.33) 


Sk  1  £klm 


A,  d3X  . 


1  8A 


m,4 


Now  look  at  the  behaviour  of  L,  under  translations.  The  orbital 

k 

angular  momentum,  calculated  in  a  reference  frame  translated  by  an  amount 

a  ,  from  a  frame  in  which  L.  is  given  by  (5.32),  is 

y 


(5.34) 


Lk  1  eklm 


x'  T' . (xx)  d3x^  ,  where 
1  m4  u 


(5.35) 


x'  =  x  -  a 


y 


All  known  fields  are  scalars  under  translations,  so 


(5.36) 


T'  (x')  =  T  (x  ) 
a8  y  y 


Changing  the  integration  variable  in  (5.34)  to  the  unprimed  coordinates, 
using  d3x^  =  d3x,  and  substituting  from  (5.35)  and  (5.36),  one  has 


(5.37) 


=  -i 

£klm 

=  -i 

klm 

II 

-i  e, 

XI  WV  d3x 

(xl+al>  Tm4(Xy>  d3x 


klm  al 


T  ,  d3x  . 
m4 


Xn  the  rest  frame  of  a  particle,  from  Appendix  1,  we  know  that  the  integral 
appearing  in  (5.37)  vanishes,  since  it  is  just  the  three-momentum,  Pm>  to 
within  a  factor.  Such  a  rest  frame  exists  when  the  system  is  described 


* 
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by  time-independent  localized  fields.  From  (5.37)  we  see  that  for  such 
solutions  =  L^.  In  such  a  case,  is  translationally  invariant,  or 
equivalently,  we  may  say  that  the  orbital  angular  momentum  is  independent 
of  the  point  about  which  it  is  calculated.  This  result  is  consistent 
with  our  expectation  that  the  orbital  angular  momentum  should  vanish  in 
all  reference  frames  in  which  the  three-momentum  is  zero.  However,  we 
have  not  been  able  to  show  that  the  integral  (5.32)  for  does  vanish 
for  all  time-independent  localized  solutions.  If  L,  is  non-zero  in  any 

K 

zero  momentum  frame,  it  would  describe  what  would  be  called  an  intrinsic, 
rather  than  an  orbital  property.  If  one  could  show  that  1^=0,  in  every 
such  zero  momentum  frame,  then  the  interpretation  of  as  orbital  angu¬ 
lar  momentum  would  be  much  more  satisfactory  than  it  presently  appears 
to  us. 

Another  point  illustrating  that  the  division,  of  the  total  angular 
momentum  into  spin  and  orbital  parts,  is  not  clear  cut,  is  that  for  a 
variety  of  Lagrangians  one  can  show  that  the  integrals  for  and  are 
proportional  for  all  localized  time- independent  solutions.  As  an  example, 
consider  the  class  of  Lagrangian  constructed  from  the  following  scalar 
quantities: 


(5.38) 


y 

=  h  A 

A 

x  =  h  z 

a. 

P  ,o 

V 

=  h  A 

A 

K  =  g  a 

o  ,p 

o,p 

,0  < 

0) 

=  A 

C  =  A  A 

0,0 

0  0 

Let  the  Lagrangian  density  have  the  form 


(5.39) 


X  =  Up  +  Vv  +  (w,x  >5  »C  >8)  » 


4 


* 
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where  U  and  V  are  constants  and  c C  is  an  arbitrary  differentiable  function 

its  arguments.  The  general  form  (5.39)  includes  all  the  Lagrangians 

studied  in  Chapter  3  for  which  exact  solutions  were  found.  The  ordinary 

electrodynamics  Lagrangian,  X  =  F  2,  is  the  case  U=-V~r-  and X  =0. 

16tt  yv  4tt 

In  Appendix  4  we  derive  the  following  results  for  the  angular  momenta  of 
time- independent ,  sufficiently  localized  solutions  of  the  field  equations 
derived  from  Lagrangians  of  the  form  (5.39): 


(A4.17) 


t  =  U 


iilxit  d3x 


1  =  (2V-U) 


txt  d3x 


J  =  £  +  t  =  2V 


txt  d3x 


y 


based  on  the  definitions  (5.32)  and  (5.33),  and  ]£=-V<f) 

As  indicated  previously,  the  interpretation  of  t,  as  the  orbital 
angular  momentum  is  satisfactory  only  if  it  vanishes  for  static,  localized 
solutions.  We  see  from  (A4.17)  that,  with  the  exception  of  the  case  U=2V, 
the  vanishing  of  L  implies  that  the  spin  angular  momentum  is  also  zero. 

So,  our  conclusion  is  that  either  there  is  no  time-independent  system  of 
the  form  (5.39)  which  can  describe  a  localized  particle  with  non-zero 
spin,  or  the  identification  of  L,  given  by  (5.31),  as  the  orbital  angular 
momentum,  is  not  justified. 

The  results  (A4.17)  also  illustrate  a  problem  mentioned  by  Belin- 
fante  (1940)  in  connection  with  the  division  of  the  total  angular  momentum 
into  spin  and  orbital  parts:  two  different  Lagrangians,  generating  the 
same  field  equations,  can  give  a  different  splitting  of  the  total  angular 
momentum.  To  take  a  specific  example,  consider  the  following  two  Lagra¬ 


ngians: 


4 

‘  • 


. 
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(5.40) 

X  =  -~ 

I  4tt 

(5.41) 

X  XI  ! 

where  X  =  X  ( 

g  „g  „  ;  g 
>a  ,a 

T  F  2  +  \  (A  )2  + 
{  4  uv  2  a, cr 


X  - L  f  I  (A  )2  +  _f' 

TT  4tt  2  ^a,p; 


>a  a  a  a 

The  field  equations  derived  from  both  Lagrangians  are 


A 

a,pp 


9 

P 


(*£L) 

9gJ 


=  *£ 
3g 


The  Lagrangian,  iCT ,  given  by  (5.40),  corresponds  to  \]=~y=~~  in 

1  4tt 

whereas  corresponds  to  U=0;  V=-~  .  We  find  from  (A4.17) 

angular  momenta  for  the  Lagrangians  cC ^  and  are 


(5.39), 
that  the 


<rIT=  o 

t11  =  -2? 

J11  =  -2* 

Vcj).  Although  the  total  angular  momenta 
J"  and  J'"1'  are  equal,  (5.42)  shows  that  unless  all  the  angular  momenta 
are  zero,  the  division  into  spin  and  orbital  pieces  depends  on  the  form 
of  the  Lagrangian. 

We  regard  the  problem  of  defining  the  spin  and  orbital  angular 
momenta  in  classical  field  theory  as  remaining  unsolved.  We  have  not 
been  able  to  show  that  a  common  definition  of  the  orbital  angular  momen¬ 
tum  vanishes  in  the  rest  frame  of  a  particle,  as  it  should.  The  defin¬ 
ition  suffers  from  a  lack  of  criteria  to  determine  what  hypersurface  to 
use  in  evaluating  the  relevant  integrals.  There  is  also  an  ambiguity 
arising  from  a  dependence  on  the  form  of  the  Lagrangian.  A  proportion- 


(5.42) 


51  =  t 
il  =  -3? 
31  =  -2? 


where  S  =  j  d3x  ;  15  =  - 
T  TT 


■ 
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ality  between  the  orbital  and  spin  angular  momenta  for  time  -independ¬ 
ent  localized  solutions  in  a  variety  of  cases,  and  the  translational 
invariance  of  the  orbital  angular  momentum,  in  the  same  cases,  suggests 
a  possible  ’intrinsic'  component  to  what  as  been  defined  as  the  orbital 
angular  momentum.  Until  these  difficulties  are  resolved,  to  us,  it  seems 
that  the  best  approach  is  to  consider  the  total  angular  momentum,  J,  in 
a  zero-momentum  frame,  (rest  frame),  as  the  intrinsic  angular  momentum 
of  the  particle.  One  then  still  has  to  resolve  the  problem  of  the  mean¬ 
ing  of  the  separate  conservation  of  and  S^.  when  the  canonical  energy- 
momentum  tensor  is  symmetric. 

5.3  Angular  Momentum  for  Exact  Solutions 

We  now  calculate  the  angular  momentum  of  the  exact  solutions  found 
in  Chapter  3.  The  Lagrangian  for  the  systems  of  Chapter  3  is 


(2.65)  <C=-h  \  |(8  ,.)2  +  Tf....2  +  |(A.  ,.)2  +  K(g)Ai  ,,  -  g2A2  +  M(g) 


4tt 


•  y 


4  yv  2  y,y' 


This  is  of  the  form  (5.39),  with  U=  -  V=  .  In  terms  of  the  variables 


and  C  introduced  in  (5.38),  jC  is 


(5.43) 


X 

4tt 


v  -  y  +  -^)2  +  K(g)o)  +  X  -  g2C  +  M(g) 


Another  Lagrangian  leading  to  the  same  field  equations  is 


(5.44) 


X  =-f 

II  4tt 


v 


+  (6~l)y  +  Km  +  X-  S2C  +  M(g) 


which  corresponds  to  U  = 


1-6 


and  V  =  -  .  When  6=0,  X  -X 


but 


4tt  -  '  4tt .  ’  I  II 

when  6=1,  we  will  have  to  distinguish  between  X^  and  cC^  when  calculat¬ 
ing  angular  momenta. 

First  consider  the  case  6=0.  From  (5.43)  in  (A4.5)  and  (A4.7), 


■ 
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we  find  the  'orbital'  angular  momentum 


(5.44) 


L.  =  e  x..  (A  <p  -K<p  )  d3x 

k  4tt  klm  J  lv  n,m  ,n  r ,m 


and  the  spin'  angular  momentum 


(5.45) 


$  =  ~  f  M  d3x  ;  £  =  -V<J) 


The  total  angular  momentum,  3,  is  given  by 


(5.46)  3  =  £  +  t  . 

When  6=0  and  K=0,  the  only  solutions  we  were  able  to  find  are 
given  by  (3.26).  The  vector  potential,  1=0,  for  these  solutions.  From 
(5.44),  (5.45),  and  (5.46)  we  see  that  if  K=0  and  1=0,  then 


£  =  $  =  3=0  . 


When  6=  0  and  K^O,  and  specifically  for  K=xgm,  we  found  solutions 
of  the  form  l=A(r)  T  ,  with  A  given  by  (3.39),  and  <J>=<Kr),  given  by  (3.40). 
Since  £--£<}>=-  for  these  solutions,  £  is  parallel  to  1.  We  find  from 

(5.45),  using  this  result,  that  £=0.  Also,  since  l=A(r)  and  <f>=<j>(r) , 
one  has 


(5.47) 


A 

n,m 


<t> 


,n 


A 

r 


6  + 
mn 

Xn  d<j> 
r  dr 


x  x 
n  m 

r 


d_  A 
dr[r 


and 


Substituting  (5.47)  into  (5.44)  and  using  eklm=_ekml»  one  finds  that 
L=0.  Combining  this  with  our  previous  result  for  S  gives  J  =  L  =  S  =  0. 
Therefore, for  all  the  exact  solutions  found  with  6=0,  the  angular  momenta 
are  all  zero.  However,  in  the  particular  case  where  K=Tg  ,  (3.39)  gives 


1  -  ■»-  :  *  m  tifjrs|  H  -1*  >  .  *1 
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A(r)=-~j-,  so  A  goes  to  zero  as  r"1,  as  r-*».  In  this  case  there  is  the 
possibility  that  the  definitions  of  the  angular  momentum  based  on  the 
canonical  and  symmetrized  tensors,  may  be  different.  From  (5.25),  these 
definitions  are  related  by 


(5.25) 


J,  =  -i 

k  klm 


(xlW,n  d3x 


One  finds  from  (5.9)  that 


(5.48) 


yvA 


J 

4tt 


A  (F  -K6  .  )  -  A  (F  .  -K6  ) 

y  Av  vA'  v  yA  yA 


Using  (5.48)  in  the  definition  (5.12)  of  g 


yvA  * 


we  get 


A  F  +  K(A  6  -A  6  J 
y  Av  A  yv  v  yA 

Substituting  (5.49)  into  (5.25)  and  using  the  results  (5.47),  one  finds 

6 

that  the  right-hand  side  of  (5.25)  vanishes,  leaving  .  There¬ 

fore,  we  may  say  unambiguously  that  the  angular  momentum  is  zero  for  all 
the  solutions  with  6=0. 

Next  we  consider  those  solutions  with  6=1  and  K=0,  given  by  (3.18) 
and  (3.29).  In  this  case,  using  the  Lagrangian,  @C  ,  given  by  (5.43), 
we  find  from  (A4.5),  that  the  static  angular  momenta  are 


(5.49) 


g. 


yvA 


4tt 


(5.50) 


_1 

4ir  eklm 


and  from  (A4.7) 

(5.51)  S 


However,  using  the  Lagrangian, 
and  (A4.7),  that 


J 


)  d  3x 
m 


EXA  ddx 

given  by  (5.44),  we  find  from  (A4.5) 


t11  =  ?n  =  0  ,  so  j11  -  0  . 
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With  considerably  more  effort,  we  can  show  that  and  as  given 
by  (5.50)  and  (5.51),  also  vanish.  To  see  this  start  by  manipulating 
the  expression  (5.50)  for  L.1  into  the  form 


(5.52)  L*=-^eklm 


(x  A.<j>  .)  ~(x. A. d>  )  .+A-<f) 

{  1  J  ,  J  ,m  1  3y,m',j  ly,m  J 


d3x 


or  in  vector  form,  using  (5.51)  in  the  last  term 


(5.53) 


t,1  »--L 

4tt 


lim 

r-x» 


Vx((^-2)r)  +  (A.rxJ)  .  d3x  -  ^  , 

.  3  - 


±  j  M  rxS  dS  J  ?xg  t-n  dSj  -  t1  , 

s  s 


where  n  is  the  outward  normal  to  the  surface  S.  We  can  take  S  to  be  a 
sphere  of  radius  r,  so  that  n  =  1^.  As  a  result,  rxn  =  0  and  the  first 
integral  in  (5.53)  vanishes.  In  the  limit  r-*»,  the  electric  field, 
for  the  solution  (3.18)  is  radial  so  that  rx^-K).  Therefore,  the  second 
integral  in  (5.53)  vanishes  as  well,  leaving 


(5.54) 


This  implies  that  the  total  angular  momentum  =  t]1  +  ^  =  0.  We  can 
go  one  step  further,  to  show  that  t1  and  21  vanish  separately.  From 
(A4.6)  in  Appendix  4: 


(A4. 6) 
(5.55) 


S,1  =  ~r~  e,  t  A1(l>  d3x  ,  so 
k  4it  klm  J  1  ,m 


S1  =  J_£ 
k  8^  klm  . 


(x  (<t>A  .-<!>  .A  ))  .  -  (<J>A  ) 

1  m, j  ,J  m  ,j  m  ,1 J 


d3x 


1 


+ e  *  (<j>  ..A-^A  ..)  d3x  , 

8tt  klm  J  1  ,JJ  m  m,j j 


where  the  second  expression  follows  as  an  identity  from  the  first,  as  can 
be  seen  by  expanding  the  derivatives  and  simplifying.  The  integrand  of 


.. 
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the  second  integral  in  (5.55)  can  be  shown  to  vanish  directly  from  the 
field  equations  (3.32.1)  and  (3.32.2).  The  first  integral  in  (5.55)  can 
be  transformed  to  a  surface  integral  at  spatial  infinity.  Using  the  ex¬ 
plicit  form  of  the  solutions  (3.18)  and  (3.21),  in  the  limit  r-*»,  one  can 
show  that  this  surface  integral  vanishes  also,  so  that  =  0.  Therefore, 
from  (5.54)  we  have 

(5.56)  t1  =  t1  =  31  =  0  . 


Before  stating  that  the  angular  momentum  of  the  solutions  is  zero 

0 

we  now  check  that  the  definition  of  the  total  angular  momentum,  J^,  based 
on  the  symmetrized  energy-momentum  tensor,  gives  the  same  zero  result. 

From  (5.25) 


(5.25) 


-1  £ 


klm 


(xiW 


n 


d3x 


with  g  ,  defined  by  (5.12)  and  (5.8).  Converting  the  integral  (5.25)  to 
yvA 

a  surface  integral  at  spatial  infinity  and  using  the  form  of  the  solutions 
(3.18)  and  (3.29),  in  the  limit  r-*30,  one  finds  that 


(5.57) 


The  detailed  calculation  of  the  result  (5.57)  is  rather  lengthy  and  will 
not  be  presented  here.  The  result  (5.57)  holds  for  both  the  Lagrangians 
and  cC-j-j. 

Xn  summary,  we  conclude  that  all  the  angular  momenta,  ^i,  and 
£j*0  zero  for  the  exact  solutions  we  have  found.  The  difficulty  of  inter— 
pretation  which  would  arise  from  a  non-zero  orbital  angular  momentum  is 
therefore  not  present.  Also,  we  have  seen  that  the  two  definitions  of 
the  total  angular  momentum,  based  on  the  canonical  and  symmetrized  energy- 


' 


■ 


. 
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momentum  tensors,  respectively,  agree  for  these  solutions.  Furthermore, 
the  two  Lagrangians  studied  which  led  to  the  same  field  equations,  also 
gave  the  result  for  the  angular  momentum. 


CHAPTER  6 


STABILITY  ANALYSIS 

6.1  Two  Criteria  For  Stability 

The  question  of  whether  or  not  solutions  of  the  field  equations 
are  stable  against  small  perturbations  is  one  which  frequently  arises  in 
the  study  of  nonlinear  classical  field  theories.  It  is  not  clear  what 
stability  requirements  need  to  be  imposed  at  the  classical  level.  One 
possibility  is  that  there  need  be  no  restrictions,  based  on  stability, 
on  the  acceptability  of  a  solution  to  a  classical  field  theory.  One 
would  then  rely  on  quantization  of  the  classical  solution  to  introduce 
some  element  of  stability.  In  this  regard  one  is  reminded  of  Bohr's  or¬ 
iginal  quantum  treatment  of  the  hydrogen  atom,  in  which  a  classical  sys¬ 
tem  with  the  same  structure  would  collapse  due  to  the  emission  of  radi¬ 
ation. 

While,  experimentally,  there  are  no  known  mesons  which  are  absol¬ 
utely  stable,  they  do  exist  for  a  finite  time  before  decaying.  This  sug¬ 
gests  that  an  appropriate  classical  model  would  be  one  in  which  the  solu¬ 
tions  are  metastable,  that  is,  the  time  scale  for  the  decay  of  a  particle¬ 
like  solution  is  long  compared  to  some  other  characteristic  time  of  inter¬ 
est.  An  example  of  a  system  for  which  metastable  solutions  exist,  has 
been  given  by  Rosen  (1965). 

A  third  alternative  in  the  treatment  of  stability  of  solutions  of 
classical  field  theories  is  to  require  that  the  particles  be  absolutely 
stable  against  small  perturbations.  In  this  chapter  we  study  two  criter¬ 
ia  that  have  been  used  in  the  literature  for  deciding  stability.  The 
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first  of  these  tests  is  based  on  the  premise  that  the  energy  of  a  stable 
static  solution  should  be  a  minimum  with  respect  to  static  variations  of 
the  fields,  about  that  solution.  The  second  test  requires  for  stability, 
of  a  static  solution,  that  any  slightly  different  solution  oscillate  a- 
bout,  rather  than  diverge  from,  the  given  static  solution  as  time  increas¬ 
es.  After  presenting  a  simple  example  to  illustrate  the  two  types  of  an¬ 
alysis  and  the  relation  between  them,  we  attempt  a  similar  study  of  one 
of  the  systems  for  which  exact  solutions  were  found  in  Chapter  3. 

We  first  look  at  an  example  of  the  application  of  the  minimum  en¬ 
ergy  requirement  with  respect  to  static  variations  about  a  static  solu¬ 
tion.  Using  this  criterion.  Derrick  (1964)  showed  the  instability  of 
solutions  for  a  class  of  Lagrangian  constructed  from  a  real  scalar  field, 
ip .  Specifically,  consider  the  following  Lagrangian: 


(6.1) 


0I>  )2 

»v 


+  r0|>) 


y 


where  F  is  an  arbitrary  function  of  the  scalar  field,  ip.  The  field  equa¬ 
tion  is 


(6.2) 


V2ip 


8  ib2  _  dr 

~  dip 


The  rest  energy  for  a  static  real  solution  ip  =  ipQ(r)  of  the  equation  (6.2) 
is 


(6.3) 


£  d3x 


(  WQ)2  +  r(*0)  )  d3x  . 


Now  calculate  the  value  of  the  integral  (6.3)  when  ipQ  is  replaced 

A 

by  a  perturbed  function,  >p,  given  by 


ip(r) 


m 


4-  e 


1  -*• 

r<r) 


+  0(e2) 


y 


(6.4) 


> 


' 


where  e  is  a  parameter  sufficiently  small  that  the  second  order  terms  can 
be  neglected.  The  function,  i|/(r)  is  arbitrary  except  that  it  is  suffic¬ 
iently  well  behaved  that  it  can  be  regarded  as  a  perturbation  to  the  term 
Expanding  E  (ij))  in  a  Taylor  series  in  e  we  find  from  (6.3)  that 

(6.5)  E  (\j/)  =  e(^q+£i^)  =  E  (i^q)  +  £  <5E  +  e2  62E  +  0(e3)  , 


where 


(6.6) 

(6.7) 


6E  =  ~ 
de 


£=0 


e=0 


and 


Using  the  fact  that  \p^  satisfies  the  field  equation  (6.2),  one 
finds  that  <$E=0.  This  is  simply  a  consequence  on  the  fact  that  the  field 
equation  was  derived  from  a  Lagrangian  through  a  variational  principle. 
Using  (6.3)  for  the  energy,  one  finds  from  the  definition  (6.7),  that 


(6.8) 


(Vi/)2  +  Oj;1) 


d2r 

dipZ 


d3x 


Now,  following  Derrick,  if,  as  a  particular  case,  we  define  the  perturbed 

A 

function  ip,  by  a  scale  transformation  of  the  original  solution  iK: 


(6.9) 


^ (r)  -  ^Q(ar)  ;  a>0  , 


1  -*  -*■  -> 

(a=l+e  corresponds  to  ip  (r)  =  r»Vi p^  =  ) 


then  one  finds  from  (6.8),  or  more  easily  by  making  a  change  of  variable 
in  (6.3),  that 


(6.10) 


M0)2  d3x  . 


Because  the  second  variation,  (6.10),  of  E  is  negative,  (for  \pQ? 0) ,  the 
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solution,  \pg,  does  not  provide  a  minimum  of  the  energy  functional  e(^), 
with  respect  to  variations  of  the  form  (6.4).  The  solution  is  therefore 
said  to  be  unstable. 

Using  the  same  example,  (6.1),  we  now  follow  the  treatment  of 
Rosen  (1965),  in  examining  a  second  criterion  for  testing  stability. 

Given  a  time- independent  solution,  ipg(r),  of  the  field  equation  (6.2),  we 
assume  that  there  exists  a  time-dependent  real  solution,  ip(r,t),  which 
differs  slightly  from  Now  we  imagine  that  due  to  some  transient  ex¬ 

ternal  influence,  the  system  described  by  the  field  i p^  is  disturbed  so  as 
to  be  replaced  by  a  system  in  which  the  field  has  the  value  ij>.  The  sys¬ 
tem  is  stable  if  all  such  solutions,  \p,  remain  near  the  original  solution, 
as  time  increases.  To  study  this  possibility  assume  the  following 
specific  form  for  ip : 

(6.11)  ^ (r , t)  =  i|>0(r)  +  £  \p J(r )  cosa>t+  0(e2)  , 

where,  again,  £  is  a  parameter  which  is  sufficiently  small  that  the  terms 

1 

of  second  order  can  be  neglected.  The  function,  ^(r),  giving  the  spatial 
dependence  of  the  perturbation,  is  determined  by  the  requirement  that  ip 
be  a  solution  of  the  field  equation  (6.2).  We  will  demand  only  that  (6.2) 
be  satisfied  only  to  first  order  in  £. 

1 

We  assume  that  o>  is  real  or  purely  imaginary,  so  that  ^(r)  is  a 
real  function.  If  a  solution  of  the  form  (6.11)  exists  for  oi  purely  imag¬ 
inary,  then  as  the  time,  t,  increases,  ip  will  eventually  differ  from  ipQ  by 
a  significant  amount.  In  this  case  we  say  that  the  solution  <p^  is  unsta¬ 
ble.  On  the  other  hand,  if  w  is  real,  the  solution  ip  performs  small  osc¬ 
illations  about  In  this  case  the  solution  is  stable  against  perturb 

ations  of  the  particular  form  (6.11),  but  we  cannot  conclude  that  the 
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system  is  actually  stable  because  there  may  exist  solutions  into  which 
ipg  could  evolve,  which  are  not  of  the  form  (6.11). 

Substituting  the  proposed  solution  (6.11)  into  the  field  equation 
(6.2),  one  finds  that  it  is  satisfied  to  first  order  in  z  if  \p^  satisfies 
the  equation: 


(6.12) 


a2r 

aip2 


\/j0  =  -  u>2  ip 


i 

o  * 


It  is  important  to  note  that  the  existence  of  a  solution,  ip^,  to  (6.12) 
does  not  guarantee  that  ip,  given  by  (6.11)  is  a  solution  of  the  nonlinear 
system  (6.2)  for  any  non-zero  value  of  the  expansion  parameter  e.  This 
is  because  we  have  considered  only  first  order  terms  in  e.  Determination 
of  the  conditions  under  which  a  solution,  ip,  of  the  form  (6.11)  exists 
to  the  full  nonlinear  system  (6.2)  is  a  difficult  mathematical  problem 
that  we  have  not  considered.  Assuming  that  such  solutions  do  exist,  the 

question  of  stability  is  thus  reduced  to  determining  the  sign  of  the 

?  1 
eigenvalue  w  ,  in  the  eigenvalue  problem,  (6.12),  for  ip^. 

Rosen  indicates  that  the  equation  (6.12)  can  be  derived  from  a 
variational  principle: 


6w2=  0  ^  ,  with 


(6.13) 


0) 


2  -  2 
^  -  0) 


l.o  ,  d2r 
(V*0)2  + 


dip ‘ 


1' 2 


0/0 

ip-ip0  0  J 


d3 


x 


(t(;J)2  d3x 


t  9  d  : 

By  we  mean  -^7  w 


for  an  arbitrary  variation  with 


0T 


given  by  (6.13). 
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Note  that  the  form  of  to2  given  by  (6.13)  is  just  that  obtained  by  mult¬ 
iplying  (6.12)  by  \J)q  and  integrating.  The  condition  <5u)2  =  0  means  that 

2  •  1 
to  is  an  extremum  when  evaluated  at  ip^.  It  may  further  be  shown  from 

the  properties  of  Sturm-Liouville  equations,^  that  this  extremum  is  a 

minimum.  Therefore,  assuming  that  a  solution  ip^  of  (6.12)  exists,  then 

if  (6.13)  is  evaluated  for  some  trial  function  \p^  instead  of  ip*,  and  it 

is  found  that  to2  is  negative,  then  to2  will  also  be  negative  for  the  exact 

solution,  ipQ.  Taking  the  trial  function  to  be  defined  by  a  radial 

scale  variation  of  the  solution  ip^,  as  in  (6.9),  one  finds  that  w2  given 

by  (6.13)  is  negative.  This  follows  from  our  previous  result  (6.10)  be- 

1  2 

cause  the  numerator  in  (6.13),  evaluated  at  ,  is  just  6ZE,  given  by 
(6.8)  and  (6.10).  Therefore,  again,  we  conclude  that  the  solution  \p^  is 
unstable. 

Using  both  the  (static)  criterion  for  instability,  that  62E<0  for 
arbitrary  static  variations  about  ipQ,  and  the  (dynamic)  criterion,  that 
m2<0  for  time-dependent  perturbations  of  satisfying  the  field  equa¬ 
tions,  we  found  that  static  localized  solutions,  ip^,  to  the  equation  (6.2), 
are  unstable.  In  this  example  there  is  a  very  close  connection  between 
the  two  tests  of  stability.  Comparing  (6.8)  to  (6.13)  we  see  that  if 
the  static  variation  ip 1  is  equal  to  the  function  appearing  in  the 
time-dependent  solution  (6.11),  then  the  following  relation  is  obtained: 

„  i,  2  52e^ 

(6.14)  w2(>0)  =  — - - 

|(^q)2  d3x 


Sturm-Liouville  differential  equations  for  several  independent  vari- 
bles  are  discussed  by  Hellwig  (1967). 
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Since  ip^  is  a  real  function,  (6.14)  shows  that  <52E  and  oj2  have  the  same 
sign,  when  evaluated  in  terms  of  It  is  important  to  note  that  this 

identity  of  sign  has  been  shown  only  for  ip^  satisfying  the  condition  (6.12) 
that  ip  given  by  (6.11)  be  a  solution  to  the  field  equation  (6.2)  to  first 
order  in  £.  However,  the  function 

(6.15)  f 1  =  r  , 

corresponding  to  the  static  scale  variation  (6.9),  does  not  satisfy  (6.12), 
in  general.  Thus,  even  though  it  happens  that  62e(iJ^)  is  negative  for  ^ 
given  by  (6.15),  we  can  only  relate  this  result  to  a  variational  approx¬ 
imation,  m2(^),  to  the  exact  value  of  w2  =  co2  C^q)  • 

The  nature  of  the  stability  of  a  particle  is  a  statement  regarding 
the  evolution  of  a  system  with  time.  Therefore,  we  regard  the  ’dynamical' 
test  of  stability  based  on  the  sign  of  m2  as  more  easily  understood  than 
the  ’static’  test  based  on  62E.  From  this  point  of  view,  the  relation 
(6.14)  may  be  regarded  as  a  justification  for  the  use  of  62E  in  deter¬ 
mining  stability.  However,  there  are  situations  in  which  62E  cannot  be 

o 

directly  related  to  m  ,  as  indicated  in  the  previous  paragraph,  and  as 
will  be  seen  when  we  apply  these  tests  to  another  example.  In  these 
cases,  unless  there  is  some  further  argument  in  favor  of  deciding  stab¬ 
ility  from  S2E,  we  feel  that  the  sign  of  co2  is  the  more  meaningful  test. 

6.2  Stability  Analysis  For  A  System  From  Chapter  3 

We  now  attempt  to  examine  the  stability  of  one  of  the  simplest 
systems  of  chapter  3,  namely,  the  case  6=K=0,  in  the  Lagrangian  (3.1). 

The  full,  time-dependent  Lagrangian  is: 


' 


1 


. 


(6.16) 


JC 


JL 

4tt 


liFpv2+I(8,M)2-*2AM+M^) 


The  corresponding  field  equations  are,  written  out  in  detail: 


->■  ->  -> 
Vx(Vxa)  +  V 


'  -  2~* 


(6.17) 


9 


9  9 -*■ 

+  -  2g2A  =  0 


9t- 


V2<f>  +  yt:(V*A)  +  2g24>  =  0 
y2g  -  +  2(A2-<J>2)g  -  N (g)  =  0  , 


where  N(g)  = 


dM 

dg* 


Static  solutions  to  (6.17)  for  <j>(r)  and  g(r),  with  A=0,  are  given 
by  (3.26).  First  look  at  the  energy  and  its  second  static  variation,  6  E 
for  solutions  of  this  type.  From  Appendix  1,  the  energy  of  solutions  to 
the  field  equations  (6.17)  in  the  static  case,  is 

e£  d3x 


(6.18)  E 


g»<f>»A 

k  > 

a 

4tt 


|(V*A)2  -  |(V<102  +  |(Vg)2  +  (4’2-A2)g2+  M 


d3x 


To  first  order  in  a  small  parameter  e,  define 


<j>(r)  =  <J>(r)  +  e  <P'(r) 


(6.19) 


g(r)  =  g(r)  +  e  g"(r) 


A(r)  =  A(r)  +  e  A  (r)  , 


where  <|>(r) ,  g(r),  and  A(r)  are  a  static  solution  (for  example,  those 
given  by  (3.26)  )  to  the  field  equations  (6.17).  Defining  62E  by  (6.7), 
that  is, 

(6.7) 


1  d2E 
<5  E  =  — 


2  de‘ 


£  =  0 


we  find  the  following  expression  for  62E: 


(6.20)  62e(£J,£)  =~ 


(VxZ')2-  (fy')2+  ($g')2+  (~+2(4,2-P))(g')2 


+  2g2  ((^')  2-(A')  2)  +  8gg'(<j><|> '”*•£') 


d  3 


Because  of  the  complexity  of  (6.20),  we  have  not  been  able  to  determine 
the  sign  of  <$2E  for  any  general  form  of  the  perturbations  g',<p',  and  A'. 
However,  in  the  particular  case  of  the  exact  solutions  (3.26),  we  were 
able  to  do  an  explicit  calculation  of  62E  for  a  static  variation  of  the 
form  + 


(6.21) 


<Kr)  “►  <j>(r)  =  s<f>(sr)  ;  s>0 

-¥  A  -> 
g(r)  g (r)  =  sg(sr) 


This  variation  is  chosen  rather  than  a  simple  scale  variation  because  it 
preserves  the  value  of  the  charge  (i.e.  the  coefficient  in  an  expansion 
of  <f>),  whereas  the  scale  variation  does  not.  Writing  s=l+e,  where  e<<l, 
comparison  with  (6.19)  shows  that  this  corresponds  to 


4> '  (r)  ~  <p  (r)  +  r*^<f> 
g"(r)  =  g(r)  +  r*^g 


One  finds  from  (6.7)  and  (6.18)  that 


(6.22) 


1  d2E 


E  "  2  de^ 


1  d2E 


e=0 


2  ds2 


s=l 


n+1  0 


En(n)  ;  n=l,2,3,  . . . 


where  E^(n)  is  the  (non-negative)  energy  given  by  (4.7).  In  the  case  of 
the  'ground  state',  n=l,  the  energy  EQ(n=l)=0,  so  that  62E=0,  but  for 
n-2,  Eq(u)>0  so  62E<0.  Therefore,  with  the  exception  of  the  solution 
for  n=l,  if  the  sign  of  62E  is  a  valid  test  for  stability,  then  we  con¬ 
clude  that  all  the  exact  solutions  are  unstable. 


t 


This  is  sometimes  called  a  dilation  transformation. 
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Next  consider  time-dependent  perturbations  of  the  solutions  g(r), 
<J>(r)  ,  and  A(r).  To  first  order  in  £,  let 

J(r,t)  =  <j> (r )  +  e  (  <J>L(r)  elU>t  +  4>X (r) *  e“1W  C  ) 

& 

(6.23)  g(r,t)  =  g(r)  +  e  (  g^r)  e1(0t  +  g*(r)*  e  1W  C  ) 

i(?,t)  E  td)  +  e  (  tl(t)  eia)t  +  tl(i)*  ), 

where  g  ,  <J>^  and  t}  are  arbitrary,  non-singular  complex  valued  functions 
which  drop  off  at  least  as  fast  as  g,  <}>,  and  A,  respectively,  as  r-*».  If 
a)  is  not  purely  imaginary,  then  substitution  of  the  expression  (6.23) 
into  the  field  equations  (6.17),  retaining  only  terms  up  to  first  order 
in  e,  and  requiring  that  the  coefficients  of  the  linearly  independent 

•k 

terms,  eia3t  and  e  separately  vanish,  gives 

(6.24.1)  Vxtfxt1)  +  iu>  (u2+2g2)ll  -  4ggL£  =  0 

(6.24.2)  V2^1  +  ha  V-t1  +  2g24> 1  +  4<f>gg1  =  0 

(6.24.3)  V2g"*  +  (w2+2(l.2-(j)2) -)g^  +  4g(A,it*-<j><|>*')  =  0  , 

as  well  as  the  complex  conjugate  of  these  equations.  If  w  is  purely 

* 

»  i(j0  ^ 

imaginary,  then  e'  =  e  and  we  get  only  one  set  of  equations,  name¬ 

ly,  the  real  part  of  the  equations  (6.24).  The  imaginary  parts  of  X1 , 

4> ■*■ ,  and  g^  are  arbitrary  in  this  case. 

If  a  sum  of  terms,  corresponding  to  different  frequencies,  w,  had 
been  included  in  the  expressions  (6.23),  then  the  equations  for  the  amp¬ 
litude  of  the  perturbation  of  one  frequency,  decouple  from  those  for  an¬ 
other.  Therefore,  even  when  more  than  one  frequency  component  is  present, 
we  need  only  consider  equations  of  the  form  (6.24).  Attempts  made  to 


■ 
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find  exact  solutions  to  (6.24)  using  the  exact  solutions,  (3.26),  for  <j> 
and  g,  with  &=0,  were  not  successful. 

Now  look  more  closely  at  the  equations  (6.24).  Multiplying 
(6.24.1)  by  tl9  adding  <J>^  times  (6.24.2),  and  subtracting  the  result 
from  g^  times  (6.24.3)  gives: 

(V*itV  -  (V)2  +  (Vg1)2  +  (g+2(*2-it2))  (gV 


(6.25) 

^fgV.*1)  = 


+  2g^f(<J>  f  -(^)2)  +  8gg 


d3 


X 


(  (J1)2+(g1)2  j  d3x 


The  linear  term  in  to  has  dropped  out  because  it  can  be  converted  to  a 
surface  integral  at  spatial  infinity.  We  have  assumed  that  the  perturb¬ 
ations  die  off  fast  enough  as  r-*»  to  make  this  integral  vanish. 

Since  the  equations  (6.24)  apply  only  when  to  is  not  purely  imag¬ 
inary,  (6.25)  also  has  this  restriction.  When  to  is  purely  imaginary,  the 
appropriate  equations  to  be  satisfied  are  the  real  parts  of  (6.24).  The 
integral  expression  for  this  negative  to2  is  exactly  of  the  form  (6.25), 
but  with  g^,  <f>\  and  replaced  by  their  real  parts,  Re(g^),  Re(4>^),  and 
Re  (X1)  ,  respectively.  Now  compare  this  result  to  the  expression  (6.20) 
for  <$  E.  If  we  let  the  real  functions  g  ,  <f>  and  A  ,  appearing  in  the 


1  ,1 


-fl 


integral  for  <$ZE,  equal  the  real  parts  of  g  ,  4>  ,  and  A  ,  respectively, 
then  we  see  that  when  to  is  purely  imaginary,  the  following  equality  holds: 


(6.26) 


2  8 it  62E 

to^  =  -j - - - 

((A')2+(g")2)d3x 


Since  A "  and  g'  are  real  functions,  (6.26)  shows  that  to2  and  62E  have  the 
same  sign, (negative) ,  when  to  is  purely  imaginary.  When  to  is  not  purely 
imaginary,  a  similar  equality  holds: 
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8  TT  $2E 

f(l1)2+(g1)2)  <*3X 

for  all  gA,  <p and  Ax  satisfying  the  equations  (6.24).  However,  as  will 
be  shown  shortly,  in  this  case  g\  <{>^,  and  ^  are  not  real  functions  so 
that  the  corresponding  62E  is  not  a  real  static  variation  of  the  type 
previously  supposed.  Also,  since  and  g^  are  not  real,  we  cannot  guar¬ 
antee  from  (6.27),  that  <$2E  and  uj2  have  the  same  sign. 

Assuming  that  g\  <f>\  and  are  real  in  (6.24)  leads  to  the  cond¬ 
itions 

-H  ->•  i 

y  (2vA  -V<|>  )  =  0 

(6.28)  y  V-A1  =  0 

1  n 

y  v  g  =  o  , 

where  y  =  Re(w)  and  v  =  Im(co) .  If  y^O,  then  we  find  from  (6.28)  that 
either  4> ^  =  0  and  =  0,  or  v=0  and  =  0.  The  latter  alterna¬ 

tive  is  not  consistent  with  our  requirement  that  the  perturbations  vanish 
at  spatial  infinity.  Therefore,  there  is  no  acceptable  non-zero  solution 
to  (6.24)  if  g\  4> ^ ,  and  A*  are  real  and  to  is  not  purely  imaginary.  So, 
under  these  conditions  we  are  not  able  to  relate  to2  to  62E  for  a  real 
static  variation. 

Similar  to  the  result  for  the  scalar  field  example,  (6.14),  the 

relations  (6.26)  and  (6.27)  between  62E  and  to2  hold  only  when  the  perturb- 
11  ->1 

ations  g  ,  <p  ,  and  A  satisfy  the  appropriate  differential  equations, 
(6.24).  However,  we  have  not  been  able  to  solve  (6.24).  To  do  practi¬ 
cal  calculations  to  determine  the  sign  of  62E  or  u)2 ,  it  would  be  conven¬ 
ient  to  be  able  to  use  a  variational  trial  function,  such  as  a  scale 
varied  solution  as  was  done  for  the  scalar  field,  ip ,  example.  This 
method  could  be  used  in  the  latter  case  because  the  variational  principle 


(6.27) 


or 


J 
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6to2  =  0,  given  in  (6.13),  is  an  extremum  having  a  minimum  character.  We 
will  show  that  the  expression  (6.25)  for  w2  is  an  extremum  with  respect 


to  a  certain  type  of  variation.  However,  the  appearance  of  derivative 
terms  in  the  numerator  with  both  signs  suggests  that  it  will  be  difficult 
to  determine  the  maximum-minimum  character  of  this  extremum.  A  further 
problem  is  that  the  relation  (6.25)  is  only  valid  if  w  is  not  purely 
imaginary  and  as  indicated  previously,  in  this  case  the  functions  g\  cf>\ 
and  A  cannot  all  be  real. 


To  show  that  w2,  given  by  (6.25),  has  some  extremum  properties. 


consider  variations  of  the  form 


(6.29) 


infinity.  Substitute  (6.29)  into  the  expression  (6.25).  For  ease  of  not¬ 
ation,  let 


(6.30) 


Then  one  finds: 


(6.31) 


1  Ef£2'£3'° 


w 


L 


a 
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Integrating  by  parts  in  the  first  term  in  (6.31)  and  dropping  the  surface 
term  gives : 


(6.32)  =-| 

=  0 


using 


+  (oj2+2  (A2~(j)2)_  )  gl  +  4g(i£.^’-<j)(j)*) 

(6.24.3)  . 


d3x 


Similarily, 


(6.33) 


 3ti)2 


3s, 


el  e2  e3  ® 


-2^<j)  •  V(f)  +  4g2c|>1<j/2^  +  Sgg^cj/2^ 


d3x 


2 

A 


v24i1+2g24,1+4<(,gg1 


*(2)  d3X  . 


Substituting  in  (6.32)  from  the  equation  (6.24.2)  gives: 


(6.34) 


V2  “-r4  { *(2)^1  d3x 


Finally, 

(6.35) 


63<o2 


3(1)2 


3s. 


el=e:2=:£3=^ 


’  2(Vxit1)-(Vx£(2))  -  4g2i1-X(2)  -  8gg114(2)  ]  d3X 


f-  | 2jiX(2)  d3 


X 


=  —  f  Vx(VxX^)-(o)2+2g2)l1-4gg1X  -X^2^  d3x  . 

J  V  V 


Substituting  from  (6.24.3)  leaves 


(6.36) 


2  io) 


V*14(2)  d3X  =  ^  J  4,17-X(2)  d3: 


Therefore,  expanding  u)2  (e ^  e3>)  a  Taylor  series  in  s^,  an^  £3* 


■ 
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we  have  from  (6.32),  (6.34)  and  (6.36): 


(6.37) 


0)" 


(£l,£2,e3^ 


=  to2 


ere2=e3=0 


2io) 


’e2*<2)^21-e3*1^.J(2)  ' 


d6x 


+  terms  of  higher  order. 


Therefore,  if  =0,  then  w2  given  by  (6.25)  is  an  extremum  with  respect 
to  the  restricted  class  of  variations  about  the  exact  solutions  g*,  cj>\ 
and  A*,  such  that  ^•^(2>=o,  because  in  this  case,  the  first  order  terms 
in  £^,  anc*  vanash  in  (6.37). 

This  need  to  impose  a  gauge  condition  such  as  ^*^=0,  in  order  to 
make  w  an  extremum,  arises  from  the  linear  terms  in  u)  in  (6.24),  which 
in  turn,  come  from  the  first  order  time  derivatives  in  the  original  field 
equations.  The  expression  (6.25)  is  not  a  'Lagrangian*  for  the  system 
(6.24).  That  is,  requiring  Sw2=0  for  arbitrary  variations  does  not  gen¬ 
erate  the  equations  (6.24).  This  is  to  be  compared  to  the  example  in¬ 
volving  the  scalar  field  ij>,  for  which  only  quadratic  terms  in  a>  appear 
and  the  expression  (6.13)  for  w2  is  a  ’Lagrangian*  for  the  equation  (6.12). 

The  explicit  appearance  of  linear  terms  in  a)  can  be  eliminated 
in  the  equations  (6.24)  if  the  fields  <j>  and  "K  satisfy  the  Lorentz  cond¬ 
ition  A  =0.  However,  as  our  theory  is  not  gauge  invariant,  we  should 
m  ,  y 

consider  other  gauges  as  well.  If  one  could  show  that  a  solution  is  un¬ 
stable  with  respect  to  perturbations  satisfying  some  additional  gauge 
restriction,  then  we  could  conclude  that  the  original  system  is  unstable. 
However,  if  one  finds  stability  with  respect  to  perturbations  in  a  spec¬ 
ial  gauge,  this  would  not  show  stability  of  the  system  with  respect  to 
unrestricted  perturbations. 
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Although  our  analysis  of  the  stability  of  solutions  of  the  system 
(6.17)  has  not  progressed  very  far,  there  are  several  things  that  have 
been  learned.  General  relations,  such  as  (6.27),  exist  between  <52E  and 
w2  if  one  computes  62E  for  complex,  rather  than  only  real,  static  vari¬ 
ations.  In  the  particular  case  that  m  is  purely  imaginary,  assuming  sol¬ 
utions  exist,  o>2  and  <52E  have  the  same  (negative)  sign.  We  have  seen 
that,  in  a  simple  case,  62E  evaluated  for  a  scale  varied  solution,  for 
example,  can  give  information  about  stability  through  its  relation  to 
resulting  from  a  minimum  variational  principle.  However,  because  of  the 
complexity  of  the  system  involving  g  and  ,  we  have  been  unable  to  est¬ 
ablish  comparable  results. 


. 

*■ 


CHAPTER  7 


SUMMARY 

We  have  seen  that  there  exist  nonlinear  field  equations  for  the 
electromagnetic  potential,  A  ,  and  a  real  scalar  field,  f,  such  that  the 

H 

charge  is  fixed  for  all  solutions  of  a  certain  class.  Specifically,  there 
are  systems  for  which  all  spherically  symmetric  solutions  with  power  ser- 
ies  expansions  of  the  form  (2.33),  have  a  charge  which  is  restricted  to 
the  values  -q,  0,  or  q,  for  some  number,  q.  However,  charged  particle 
solutions  of  this  type  are  found  to  interact  with  one  another  according 
to  a  force  which  drops  off  more  rapidly  than  the  inverse  square  of  their 
separation.  Other  systems,  with  a  (possibly)  larger  spectrum  of  discrete 
charges,  have  the  proper  Coulomb  force  law  for  particles  momentarily  at 
rest.  Furthermore,  exact  solutions  were  found  to  the  field  equations  of 
this  type,  which  are  finite  everywhere  and  have  a  finite,  non-negative 
rest  energy. 

These  results  were  derived  in  what  is  essentially  a  phenomenolog¬ 
ical  approach  to  finding  a  nonlinear  modification  of  the  Maxwell  equa¬ 
tions  suitable  for  describing  extended  particles.  An  investigation  of 
a  detailed  model  of  the  background  radiation  field,  assumed  to  give  rise 
to  the  scalar  field,  f,  might  delimit  the  possible  form  for  the  nonlinear 
contributions.  A  more  fundamental  theory  might  also  lead  to  definite 
statements  about  the  asymptotic  boundary  condition  appropriate  for  the 
field  f  (  or  g  ) ,  with  important  consequences  for  the  determination  of 
the  allowed  charge  values.  We  regard  the  present  work  as  an  exploratory 
study,  to  give  an  indication  of  the  potential  of  the  concept  of  fixing 
the  charge  through  the  field  equations.  In  this  spirit,  there  are  still 
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a  number  of  questions  to  be  examined.  In  addition  to  the  possibility  of 
looking  at  other  forms  for  the  Lagrangian,  even  for  the  systems  of  the 
type  already  considered,  there  remains  the  important  question  of  whether 
or  not  the  charge  can  be  fixed  for  time-dependent  solutions.  One  might 
also  examine  the  relation  between  the  charge  fixing  problem  for  angular 
dependent  solutions  and  the  simpler  case  of  spherical  symmetry. 

While  some  of  the  systems  we  have  studied  have  many  desirable 
features,  such  as  a  discrete  charge  spectrum  and  finite,  non-negative 
energy,  there  are  improvements  or  generalizations  to  be  sought  in  other 
areas.  The  angular  momentum  for  all  the  exact  solutions  that  we  found, 
turned  out  to  be  zero.  It  would  be  interesting  to  find  a  classical  field 
theory  model  in  which  there  are  particle-like  solutions  with  non-zero 
total  angular  momentum.  In  this  regard,  one  might  also  look  for  solutions 
which  are  periodic  in  time,  rather  than  static. 

Also  concerning  angular  momentum,  we  have  seen  that  there  are 
problems  associated  with  the  division  of  the  total  angular  momentum  into 
spin  and  orbital  components,  in  the  case  where  the  canonical  energy-mom¬ 
entum  tensor  is  not  symmetric.  These  include  the  fact  that  the  separa¬ 
tion  depends  on  the  choice  of  hyper surf ace  used  in  the  definitions,  and 

■j* 

on  the  specific  form  of  equivalent  Lagrangians.  The  observed  propor¬ 
tionality  of  the  spin  and  orbital  angular  momenta  in  some  cases  of  inter¬ 
est,  the  translational  invariance  of  the  orbital  angular  momentum  for 
localized  solutions  in  a  zero-momentum  frame,  and  the  failure  to  exhibit 
a  general  proof  that  the  orbital  angular  momentum  vanishes  for  such  solu¬ 
tions,  all  cast  doubt  on  the  definitions  used  for  the  quantities  involved. 


t 


Equivalent  Lagrangians  are  those  generating  the  same  field  equations. 
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If  a  zero  momentum,  localized  solution  was  found,  it  would  provide  mot¬ 
ivation  for  an  attempt  to  find  a  new  definition  for  orbital  angular  mom¬ 
entum  in  classical  field  theory.  In  any  case,  further  study  in  this  area 
seems  indicated. 

Finally,  with  regard  to  the  stability  of  solutions  in  classical 
field  theory,  we  were  unable  to  make  definite  statements  about  the  part¬ 
icular  systems  of  interest  in  this  paper.  We  were  not  able  to  conclude 
that  the  static  criterion  62E<0,  is  equivalent  to  the  dynamical  criterion 
for  instability,  w2<0,  for  these  systems,  although  in  certain  cases  the 
two  tests  can  be  shown  to  agree.  Even  if  one  could  decide  on  the  nature 
of  the  stability  of  a  solution  to  a  classical  field  theory,  one  still 
has  to  answer  the  question  of  whether  or  not  stability  should  even  be 
required  at  the  classical  level.  There  is  much  room  for  investigation 
in  the  whole  area  of  stability. 
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APPENDIX  1 


FOUR-MOMENTUM  IN  THE  TIME-INDEPENDENT  CASE 


Al.l  Definition  of  the  Four -momentum 


The  invariance  of  the  Lagrangian,  for  a  closed  system,  under  trans¬ 
lations  in  space  and  time  leads  directly  to  the  conservation  of  the  four- 
•f* 

momentum,  P  . 

y 


(Al.l) 


y 


a 

is  the  integral  of  the  normal  component  of  the  canonical  energy-momentum 

tensor,  T  ,  over  a  space-like  surface,  a.  The  tensor,  T  ,  is  defined 
yv  yv 

in  terms  of  the  Lagrangian  density,  •£=#£( A  ,  A  ,f,f  )by: 

y  y,v  >v' 


(A1.2) 


T  =  6  eC.  -  k 
yv  yv  a,y  3A 


3*C  f 


a,v 


»*  3£,v  ' 


If  does  not  depend  explicitly  on  the  coordinates  x  ,  as  is  assumed 

y 

here,  then  T  =  0  ,  when  evaluated  for  fields  satisfying  the  field 
yv,  v 

equations.  This  condition  implies  that  P  is  independent  of  the  surface, 

y 

a.  This  is  an  expression  of  the  conservation  of  the  energy-momentum 


four-vector  P  . 

y 


If  we  take  a  to  be  a  surface  of  constant  time,  then 


(A1.3) 


P  =  -i 

y 


T  .  d3x  . 
y4 


The  space  components,  P  ,  define  the  three-momentum,  ?,  and  the  energy, 


E,  is  given  by 
(A1.4) 


E  =  -i  P,  =  - 
4 


T. ,  d3x 
44 


t 


See,  for  example,  Aharoni  (1965),  Barut  (1964),  or  Rzewuski  (1967). 
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A1.2  Three -momentum  in  the  Static  Case 


In  a  reference  frame  in  which  the  fields  are  independent  of  time, 

one  can  show  that  the  components,  P  ,  of  the  three-momentum,  vanish.  We 

will  call  this  a  rest  frame  of  the  system. 

To  show  that  P^=0  for  time-independent  systems,  we  introduce  a 

symmetrical  energy-momentum  tensor  0  ,  which  is  known  to  exist  for 

yv 

Lorentz  invariant  Lagrangians.  The  tensor  0  =  0  ,  is  defined  such 

yv  vy 

that  0  =0  and  it  gives  the  same  momentum  when  integrated,  as  does 

yv,v 

the  canonical  tensor,  T  : 

yv 


(A1.5) 


pk  =  -1 


Tk4  d3x 


-1 


0k4  d3x 


=  “I 


®4k  d3x 


But,  this  last  integral  is  here  shown  to  vanish  under  quite  general 

circumstances.  Using  0  =0  we  have: 

y  v ,  v 


(A1.6) 


0  .  d3x  = 
yk 


(x,  0  )  d3x 

k  yv  ,v 


(x.  0  .)  .  d3x  ,  in  the  static  case, 
k  yj  '  ,j 


=  I  x^0  .  dS .  ,  where  S  is  a  surface  at  r=°° 

y  3  3 

s 

_3 

=  0  ,  if  0  .  falls  off  faster  than  r  as  r-*». 

yj 

So,  P,  =0. 
k 

Therefore  we  can  always  find  a  rest  frame  for  the  system,  if  there 


is  a  time-independent  solution  in  which  the  energy-momentum  tensor  falls 

_3 

off  faster  than  r  as  r  -*■  (A1.6)  also  shows  that  the  self  stress. 


0  ,  d3x  ,  vanishes  under  these  conditions. 

jk 


0  was  originally  introduced  by  Belinfante  (1940)  and  Rosenfeld  (1940), 
but  one  could  also  see  any  of  the  books  in  the  previous  footnote. 


■ 


. 
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A1 . 3  Rest  Energy 


We  now  look  at  the  energy,  E  =  - 


T..  d3x 
44 


In  the  case  of  a 


time-independent  system,  E  is  the  rest  energy  and  is  given  by 


E  =  - 


•£  d3x 


Look  at  the  particular  example  considered  in  Chapter  2.  In  the 
static  case,  with  the  vector  potential,  X  =  0  ,  the  Lagrangian  density, 
is  given  by  (2.16) 


(2.16) 


\  B  (Vg)2  -  \  D  (V*)2  +  G(4,,g)  |  , 


where  B  and  D  are  constants.  The  fact  that  the  field  equations  we  con¬ 
sider  are  derived  from  a  Lagrangian  through  a  variational  principle, 
allows  one  to  obtain  integral  relations  satisfied  by  solutions  of  the 


equations.  In  particular,  the  static  Lagrangian,  L  = 


jO  d3x  =  -  E, 


is  an  extremum  for  a  solution  of  the  field  equations,  compared  to  func¬ 
tions  obtained  from  this  solution  by  scale  or  amplitude  variations.  In 
the  expression  (2.16)  replace  <J>(r)  by  v4>(Ar),  and  g(r)  by  yg(Ar),  with 
A>0,  and  calculate  the  Lagrangian 


(A1.7)  L  =  L(y , v, A) 


_1_ 

4n 


(  fv2(Vg(X?))2-  fv2(?,|>(Xr))2  +  G(gg(Xr)  ,v$(Xr»)  d3x 


Assuming  that  4> (r)  and  g(r)  are  solutions  of  the  field  equations,  the 
extremum  conditions  mentioned  above  take  the  form 


(A1.8) 


8  L 

9y 


8  L 
8v 


8  L 

8A 


=  0 


where  all  the  derivatives  are  evaluated  at  y=v=A=l. 

Carrying  out  the  differentiations  in  (A1.8)  gives  the  following 
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integral  relations: 


(A1.9) 


(  -  D  (V<f>)2  +  t  ||  )  d3x  =0 
(  B  (Vg)2  +  g  ||  )  d3x  =0 
(  B  (Vg)2  -  D  (V<j>)2  +  6  G  )  d3x  =0  . 


Taking  a  general  linear  combination  of  the  integrals  in  (A1.9), 
one  finds  that  the  rest  energy,  E  =  -  L,  can  be  expressed  as 


(A1.10)  E  =  ±  j  (  ({H3+*)  B  (Vg)2  -  (|+ct+£)  D  (V4,)2 


+  “  ♦  If  +  6  8  If  +  (1+3^  G 


d3x 


Here,  a,  B,  and  y,  are  arbitrary  constant  parameters.  In  the  particular 
case  where  a  =  B  =  (1+y)  ,  we  get 


(Al.ll)  E  = 


4tt 


— \  (i+y)  (  ‘f’  +  8  )  +  (1+3y)  g 


d3: 


Choosing  y  =  -1  gives  the  particularly  simple  expression 


(A1.12) 


E  = 


2tt 


G  d3x 


In  Chapter  3  we  studied  systems  based  on  the  Lagrangian  density 


(Al. 13) 


4tt 


1 


2  (VxX)2  -  \  +  \  (^g)2  +  |  (v4) 2 


+  (<f>2-£2)  g2  +  K(g)  +  M(g) 


Using  the  technique  of  deriving  integral  relations  from  scale  and  ampli¬ 
tude  variations,  one  finds  in  this  case  that  the  energy  of  solutions  of 
the  field  equations  can  be  written  as 


■ 


. 


. 
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(A1.14)  E 


Here,  a,  3, 
3=0  and  a 

(Al. 15) 


/:  d3: 


4tt 


(  -kvxA)2  +  |(V-A)  )  (l+2y+A)  -  |(V4>)2(l+2a+A) 


+  j(Vg)2  (1+23+A)  +  K  V«£  (1+7+2 A)  +  3  g  v4 


+  <})2g2  (l+2a+23+3A)  -  £2g2  (l+23+2y+3A) 


+  M  (1+3  A)  +  3  g 


dM 

dg 


d3x 


7,  and  A  are  arbitrary  constants.  With  the  particular  choice 
=  7  =  A  =  — j  ,  we  find  from  (A1.14) 


E  = 


12it 


(Vg)2  +  2  (i2-«f)2)  g2 


d3] 
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APPENDIX  2 


DERIVATION  OF  CONDITIONS  TO  FIX  THE  CHARGE 

In  this  appendix  we  derive  some  of  the  conditions  which  may  be 
required  for  the  Lagrangian 

(2.16)  }b  (Vg)2  -  |  D  (V^.)2  +G(<|>,g) 

to  lead  to  a  fixed  charge  system.  We  assume  that  G(<j>,g)  is  given  by 


(2.20) 

G(4>  »g)  ~  I  Gmn  <f>m  gn  • 

m,n 

We  look  for  solutions  of  the  form 


(2.34) 

00  b 

<t>  =  y  — 

L.  n 
n=l  r 

Z  Sn 

8  =  \  IT 

n=  1  r 

to  the  field  equations 


(A2.1) 

1  3G  n  .m  n 

8  =  I  IE  '  L  °mn  ^  8 

m,n 

(A2.2) 

1  3G  y*  _  .m  n 

^  D  3(J>  ^  mn  ^  ® 

m,n 

Substitution  of  the  series  (2 . 34)  into  the  field  equations  (A2.1)  and 
(A2.2)  and  equating  the  coefficients  of  like  powers  of  r  on  each  side 
gives  the  following  conditions;  (  P^,  Sm»  and  Tm  are  defined  by  equa 


tions  (2.36) 

through  (2.39).  ) 

(A2.3.0) 

0 

0 

0 

ti 

0 

PU 

11 

0 

(A2.3.1) 

0  =  pi  =  “io  +  e  si 
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< 


- 

Ill 


(A2. 3. 2)  0  =  P2b2  +  S^-eb^ 


(A2.3.3) 

0  =  P3b3  +  2P2b1b2  +  S2b1(g2-eb2)  +  S1(g3-eb3) 

(A2.3.4) 

2g2  =  P4bJ  +  3P3b2b2  +  2P2b3b3  +  S3b21(g2-eb2) 

+  S2bl +  a20b2+allb282+a0282 

(A2. 3 .n) 

*sn(n-l)  gm  =  Pb"  +  nP  b"  V+  .  .  .  +  S  <g  -eb  )  +  ... 
n~Z  n  1  n-I  1  Z  Inn 

and. 

(A2.4.0) 

o 

O' 

II 

o 

(A2.4.1) 

°  =  Ql  =  6l0  +  e  T1 

(A2.4.2) 

0  =  Q2b2  +  T1(g2-eb2) 

(A2.4.3) 

0  =  Q3b3  +  2Q2b1b2  +  T2bl(g2"Eb2)  +  Tl(s3_eb3) 

(A2.4.n) 

^n(n-l)  b  =  Ob"  +  nQ  ,b"_2b  +  .  .  .  +  T  (g  -eb  )  +  ... 

n~Z  n  1  n-1  1  Z  Inn 

Here,  equations  (A2.3.n)  and  (A2.4.n)  are  obtained  from  the  r  n  term  in 
the  g  and  <f>  equations,  respectively. 

We  want  the  solutions  of  these  field  equations  to  have  finite 
energy.  From  Appendix  1,  the  energy  may  be  expressed  as 


(A1.12) 

* 

E  =  G  d3x 

Ztt 

. 


If  G  goes  to  zero  slower  than  r“4 ,  as  r  -*  °°,  the  energy  will  not  be  fin¬ 
ite.  To  examine  the  asymptotic  properties  of  G,  substitute  the  expansions, 
(2.34)  for  <j>  and  g  into  the  expression  (2.20)  for  G.  One  finds 
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G  ■  G00  +  7  (  eBVDQ0  >bl  +  2^  (  eBP1-DQ1 


+  3^t  (  eBP2-DQ2  )bj  +  OCr'4). 


If  Gqq  =  0,  then  one  finds,  using  the  field  equations,  that  all  terms  up 
to  and  including  r~3,  are  identically  zero,  without  placing  any  further 
constraints  on  the  parameters  appearing  in  the  Lagrangian.  Therefore, 
the  only  condition  we  get  from  the  requirement  of  finite  energy,  is 


=  0. 


We  now  show  that  the  field  equations  can  fix  the  charge,  if  the 
Lagrangian  is  chosen  properly.  Certain  conditions  must  be  satisfied  by 
the  parameters  in  the  Lagrangian,  if  solutions  with  non-zero  charge,  b^, 
are  to  exist,  and  in  particular,  if  this  charge  is  to  have  a  fixed  mag¬ 
nitude.  As  a  starting  point,  in  this  section,  we  assume  that  the  para¬ 
meter,  T^  =  3  ^  0.  Also,  throughout,  it  is  assumed  that  b^  ^  0. 

Equations  (A2.3.0)  and  (A2.4.0)  require  the  vanishing  of  the  coef¬ 
ficients  01qq  and  3qq>  the  terms  in  the  field  equations  which  are  indep¬ 
endent  of  <J>  and  g.  Equation  (A2.3.1)  serves  to  determine  the  value 


(A2.5) 

Consistency  of  equations 
(A2.6) 

Eliminating  g^-eb^  from 
and  (A2.4.3)  gives: 


(A2.3.1)  and  (A2.4.1) 

aioTi  =  eiosi  * 

(A2.3.2)  and  (A2.4.2); 


requires 


and  g^-eb^ 


from 


(A2.3.3) 


(A2.7) 


T1P2  =  S1Q2  >  3nd 


' 
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(A2.8) 


T  P  +  T  P  = 
13  2  2 


SlV  S2Q2 


The  requirements  listed  thus  far  are  necessary  for  the  existence  of 
solutions  with  non-zero  charge,  but  the  question  of  fixing  the  charge 
does  not  arise  until  one  considers  equations  (A2.3.A)  and  (A2.4.4). 
Eliminating  g^-eb^  from  these  equations  gives 

(A2.9)  -2  (S1-eT1)b2+Q2b| 


If  the  coefficient  of  is  non-zero,  and  if  S^-eT^  =  0,  the  above  exp¬ 
ression  involves  only  b^,  and  therefore  fixes  the  charge.  If  S^-eT^  4  0, 
the  equations  (A2.3.4)  and  (A2.4.4)  do  not  fix  the  charge  and  the  higher 
order  equations  (A2.3.5)  and  (A2.4.5)  must  be  considered.  However,  it 
was  found  that  the  charge  did  not  appear  to  be  fixed  by  these  higher  equ¬ 
ations  unless  S^-eT^  =  0. 

In  addition  to  the  constraints  (A2.6),  (A2.7),  and  (A2.8)  already 
mentioned,  the  parameters  are  subject  to  the  condition 


(2.40)  mBP  =  -DT  +  eBS 

m  m  m 

that  the  field  equations  be  derivable  from  a  Lagrangian.  When  m  =  1, 
this  relation,  combined  with  equations  (A2.3.1)  and  (A2.4.1),  and  the 
condition  S^-eTj  =  0,  gives  the  result 

(A2.10)  D  =  e2  B  . 

To  summarize,  if  the  conditions  (A2.6),  (A2.7),  (A2.8),  and  (2.40) 
are  satisfied,  with  =  e^,  and  e  determined  by  (A2.5),  then  all  charged 


. 
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solutions,  if  any  exist,  with  a  power  series  expansion  of  the  assumed 
form,  have  the  same  magnitude  of  charge,  q  =  |b^|,  where 

-2Q? 

(A2.ll)  q2  - - i - 

(TlP4-slQ4)+(T2P3-S2S)+(T3P2-S3Q2)  +  (Tla02-Slg02)4 

To  fix  the  charge  at  a  non-zero  value,  using  equations  (A2.3.4) 
and  (A2.4.4),  the  parameter  must  be  non-zero.  If  =  0,  one  finds 
that  it  is  still  possible  to  fix  the  charge  by  using  equations  (A2.3,5) 
and  (A2.4.5).  In  this  case,  the  necessary  conditions  are: 

(A2.12)  p0=Pi=P2=0=Qo=Q1=Q2  ;  P3=eQ3  ;  P4=eQ4  ;  Sl=eTl  ;  and  S2=eT2  * 

The  value  of  the  fixed  charge  is  q,  where 

~6Q 

(A2.13)  q2  =  - - - T  • 

T1(P5-£Q5)  -  3Q32 

Investigation  of  the  case  =  0,  did  not  reveal  any  further 
systems  in  which  all  the  spherically  symmetric  solutions,  of  the  assumed 
form,  have  the  same  charge. 


. 


. 


APPENDIX  3 


IMPOSING  THE  CONDITION  V-t  =  0 


Given  the  expression  (3.29),  reproduced  below,  for  the  Cartesian 

components,  A^,  of  the  vector  potential,  we  want  to  impose  the  supple- 

-+ 

mentary  condition  (3.28):  V*A  =  A.  .  =  0.  Here  A.  is  given  by 

1  y  1  X 

+  n-l  1 

(3.29)  A±(r)  =  l  l  (k.)lra  Ylm(6,t?)  ^(r)  ;  n-l, 2,3...  . 

1=0  m=-l 

For  small  values  of  n,  when  there  are  few  terms  in  the  sum  (3.29),  it  is 
relatively  easy  to  calculate  A^  but  for  larger  n  the  problem  is  much 
more  difficult.  An  alternative  way  to  proceed  is  based  on  the  equation 
satisfied  by  the  function  (3.29),  namely 

(3.27)  V2it  =  ~2g2l  . 

Taking  the  divergence  of  (3.27),  using  the  vector  identity 

-  V(V*it)  =  - Vx(Vx^) 


and  the  fact  that  the  divergence  of  a  curl  vanishes,  one  finds 


(A3. 1) 


0  =  V2(V«£)  +  2g2(v-£)  +  4  git* 


4- 

vg 


This  shows,  that  for  solutions  in  which  V«&=0,  and  g=g(r)^0,  the  radial 
component,  A  ,  of  the  vector  potential,  vanishes. 

When  n=l,  so  only  the  1=0  term  contributes  to  (3.29),  one  finds 
the  condition,  Ar=0,  can  be  satisfied  at  all  points  in  space  only  when 
it=0.  When  n=2,  a  non-zero  result  is  obtained,  of  the  form 


(A3. 2) 


t  = 


txr 

R3 


txr 


(r2+a2) 
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where  k  is  an  arbitrary  real  constant  vector.  The  form  (A3. 2)  is  of  part 
icular  interest  because  in  the  limit  a -*■(),  or  r-*o,  this  is  exactly  the 
form  of  the  vector  potential  of  a  particle  with  a  magnetic  dipole  moment, 
lc.  For  all  n-2,  solutions  exist  with  the  dipole  structure,  character- 
ized  by  the  icxr  factor  and  the  ^  asymtotic  (r-*»)  behaviour.  Specifical¬ 


ly, 


A  =  4*4  F(ti+2,-n+2;|;ff) 


R3 


is  a  solution.  Note  this  includes  contributions  only  from  the  1=1  terms 
in  the  expansion  (3.29),  so  we  expect  that  there  exist  more  general  solu¬ 
tions.  This  is  found  to  be  the  case. 

When  n=3,  the  condition  A  =0  no  longer  is  sufficient  to  ensure 

—y  -y 

that  the  constraint  V*A=0,  is  satisfied.  The  form  of  the  solution  for 

n=2  suggests  a  different  way  to  proceed.  Choosing  a  coordinate  system 

in  which  k  is  along  the  z  axis,  the  expression  (A3. 2)  is  of  the  form 
->•  ->  ~y 

A(r)  =  A(r)  1^,  where  is  a  unit  vector  in  the  direction  of  increasing 
fe.  To  generalize  this  result,  we  look  for  a  solution  of  the  form 


(A3.  3) 


A(r)  =  A(r ,0)  T 


for  which  the  condition  V*A=0  is  identically  satisfied, 
Separation  of  variables  in  the  equation 


(3.7.1) 


Vx(VxA)  =  2g2A 


can  be  used  to  find  a  solution  of  the  form  (A3. 3).  Let  A(r,0)  —  X(r)Y(0) 
Writing  the  separation  constant  as  1(1+1)  for  convenience,  one  finds 
that  X(r)  and  Y(0)  satisfy  the  equations 


(A3. 4) 


1  d^(rX) 

7P 


(  2c2 


R2 


n  _  1(1+1) 


X  =  0 


1 


■ 


- 


■ 
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and 
(A3. 5) 


1  d 


sin0  d0 


(sine  ~  )  +  (  1(1+1)  -  dbr )  Y  ■  0  . 


sin20 


respectively.  Here,  we  have  used  the  expression  g  = 


n 


f4n2-l] 

h 

a 

2 

a  +r 


The  equation  (A3. 4)  for  X  has  already  been  solved  because  it  is  identical 
in  structure  to  (3.11).  It  has  the  regular  solutions 

(A3. 6)  X(r)  =  X.  (r)  =  <f>  _  ;  1=0, 1,2,  ...  ,  n-1  ;  n=0,l,2,  ...  , 

i  nl 

with  <f>  _  given  by  (3.20). 
nl 

In  treating  the  Y  equation,  (A3. 5),  it  was  found  that  the  substi¬ 
tution  Y(0)  =  sinQ  W(6)  reduces  the  equation  to  a  hypergeometric  equation 
of  the  form  (3.14)  for  W,  in  the  independent  variable  sin  9.  This  leads 
to  the  solution 

(A3. 7)  Y(0)  =  sin0  F(  14~, | 2 ; sin20)  ;  1=1, 2, 3,  ...  . 

For  1=0,  a  solution  exists  of  the  form 

sin0 


Y(0)  = 


1±COS0 

but  we  exclude  this  case  because  this  solution  is  not  continuous  in  the 

closed  interval  (0 , tt)  .  Y(0)  in  (A3. 7)  can  also  be  expressed  in  terms  of 

+ 

Legendre  polynomials 


(A3. 8) 


Y(0)  = 


2  sine  d  pi(cos6> 
1(1+1)  d  (cos0) 


•  1  =  1  2 

)  1  j  i.  j  •  •  •  • 


one 


Combining  the  results  (  A3. 6)  and  (A3. 7)  for  X  and  Y  respectively, 
finds  that  t.  =  A(r,0)  t .  is  a  solution  to  the  equation  (3.7.1),  with 


n 


g  =  ^2  ,  if  A(r , 0)  is  given  by 


t 


Erdelyi  et.al.  (1953),  p  147 


' 


*'  . 
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(A3. 9) 


A(r,e) 


l  K1  ~2Th  F(l^+l.l-n+l;l+|;f|) 

1  K 

X  Sine  F(l+|,-|--^-;2;sin2e) 


y 


where  are  constants  that  may  depend  on  n.  n  can  take  on  the  values 
n- 2^3)4}  ...  . 

Alternatively,  in  terms  of  orthogonal  polynomials,  from  (3.22) 
and  (A3. 8),  we  find 


n-1 

(A3. 10)  A(r ,0)  =  l 

1=1 


1  r21+1 


,1+1 

n-1-1 


'r2-a2l 

r2+az 


sin9 


d  P1(cos0) 
d  (cos0) 


where 

"  _  2  1 
Knl  “  1(1+1)  ,  n+1  .  Knl  * 

4i+i; 

It  should  be  noted  that  the  preceding  expressions  may  not  be  the 
general  solution  to  the  equation  (3.7.1)  for  A.  For  n>2  it  may  even  be 
possible  to  find  a  solution  satisfying  V»A=0  in  which  A  is  not  of  the 
form  (A3. 3). 


■ 


. 


APPENDIX  4 


RELATION  BETWEEN  ORBITAL  AND  SPIN  ANGULAR  MOMENTA 


We  want  to  calculate  the  spin  and  orbital  angular  momenta  for  sys¬ 
tems  described  by  Lagrangians  of  the  form 


(5.39) 


=  Uy  +  Vv  +X(“»X»C,C,g), 


where  U  and  V  are  constants  and  y,v,to,x,£,  and  C  are  the  scalar  functions 
of  the  fields  A^  and  g  defined  by  (5.38).  From  (5.38)  and  (5.39)  one  can 
calculate 


(A4.1) 

(A4.2) 


aX 

3A 


=  U  A, 


a  ,8 


8  ,a 


+  V  A 


3.C 


a, 8  3u)  a8 


a<£  _  a#C  .  a«C  A 

3g  „  3x  8,a  Aa 


Using  (A4.1)  and  (A4.2),  one  finds  the  field  equations  derived  from  the 
Lagrangian  (5.39)  are 


U  Ac  0  4-  V  A  OD  + 
8,a8  a, 88 


faX 

3o) 


(A4.3) 


•  a 


fl«C  +  »£.  A 

Ux  ,a  35  a 


„  +  2^-A 

3C  8,a  3C  a 

_  a£ 

8g 

»ot 


The  other  quantity  we  will  need  in  the  calculation  of  the  angular  momenta 
is  the  canonical  energy-momentum  tensor 


(5.10) 


3<r 


T  =  6  -  A  t\ , 

yv  yv  o,y  3A 


-  g 


a  «C 


a,v 


,v»  9g 


v 


In  particular,  for  a  time-independent  case,  we  have 

9«£ 


(A4.4) 


T  ,  =  -A 
m4 


=  -i 


a ,  m  3  A  . 

-  g 

a,4 

U  A. 

+ 

,m  ag 


,4 


1  r  .  3  <J>  +  | jr-  <f>g 

n,m  ,n  3u)  ,m  ot, 


,m 
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Substituting  (A4.4)  and  (A4.1)  into  the  definitions  (5.32)  and 
(5.33)  of  the  orbital  and  spin  angular  momenta  and  S^,  respectively, 
one  finds,  in  the  static  case: 


(A4.5) 

(A4.6) 


Lk  = 


klm 


X. 


UA  <j>  +  (j, 


n,mT  ,n  9cu  ,m  9£ 


,  9«C  , 

+  4»g 


m 


d3x  and 


S  =  U  e.  _ 
k  klm 


d3: 

1T  ,m 


,  or. 


using  the  form  of  the  electric  field  in  the  static  case:  E  =  -*V<j), 
we  get  the  more  familiar  expression 


(A4.7) 


t  =  U 


d3x 

4 


Using  the  field  equations  one  can  manipulate  the  expression  for 
into  a  form  involving  S^.  The  field  equations  (A4.3)  in  the  static 
case  are: 


(A4.8) 

(A4.9) 


UA  +  VA  + 
n,mn  m,nn 


'a*' 


9w 


-  g  -  2  A  =  0 
,m  9£  ,m  9c  m 


Vd>  -  2  <t>  =  0 

Y,nn  9?  * 


Subtracting  A^  times  (A4.9)  from  <j>  times  (A4.8)  and  integrating  gives 


(A4.10)  0  =  e. 


klm 


<b  fUA  +VA  + 
v  n,mn  m,nn 


'zjC  fcC  ,  TTA  . 

9o)  __  9£  *%m  m  ,nn 


,m 


d3x. 


Using  the  expression  (A4.5)  for  L  ,  the  relation  (A4.10)  becomes 


(A4.U)  Lk=-eklm 


X, 


U(<f>A  )  + 

n,m  ,n 


a jC 

4>=- 

dO) 


,m 


m ,  nn  m  ,  nn  j 


ddx 


Now  look  at  each  of  the  terms  in  (A4.ll)  separately.  The  term 
multiplied  by  U  can  be  written  as 


-  e 


klm 


x_  (<f>A  )  d  3x  =  -e 

1  n,m  ,n 


klm 


(Xi  <f*Ar.  _.)  _  "  (MJ  +4>  A 
1  n,m  ,n  1  ,m  ,m  1 


d3x 


. 


. 


jj  :  J 

. 


■ 
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(A4.12)  (continued) 


=  -e 


klm 


(x  4>A  )  -  (<j>A. ) 

1  n,m  ,n  Y  r  ,m 


d  ^3 


U  Sk  * 


using  the  expression  (A4.6)  for  S  .  The  middle  term  in  (A4.ll)  is 


(A4.13)  -e 


klm 


X, 


v  3u> 


d3x  =  -e 


.m 


klm 


» 

\ _ 

J 

d3x  , 


,m 


because  .  Finally,  the  term  in  (A4.ll)  multiplied  by  V  is 


-e 


klm 


x-  (4>A  -A(b  )  d3x 
1  m,nn  mT,nn 

(x  (cJ)A  -A  d>  ))  -(<f»A  )  -+2(f>  , A 

v  1  m,n  m  ,n  '  ,n  m  ,1  ,1  i 


-  -e 


klm 


d3x 


=  — £ 


klm  . 


(X1  (M  -A<J>  ))  ~(4>A  )  ,  d3x  +“S  , 

JL  in  y  n  m  y  n  y  n  m  y  JL^  u  R 


using  (A4.6)  for  S^_  once  again. 


Substituting  (A4.12),  (A4.13),  and  (A4.14)  into  the  expression 


(A4.ll)  for  gives 


(Ai.iS)  Lk=(f-l)sk 


-  £ 


klm  . 


klm 


,('(xilr+<v-u>Ai> 


,m 


d3x 


x,  (U<f>A  +V(J>A  -VA  <j>  ) 

x  n )  in  hi  y  n  m  ,n  ^ 


»n 


d^x 


The  integrals  appearing  explicitly  in  (A4.15)  can  all  be  converted 
to  surface  integrals  at  spatial  infinity  and  will  therefore  vanish  if  the 


fields  A^  and  g  go  to  zero  fast  enough  as  r-*°°.  For  example,  if,  as  i-*”, 

cj)  goes  as  X  goes  to  zero  faster  than  ~  and  goes  to  zero  faster 

than  -V,  then  the  integrals  vanish.  Therefore,  for  time  -  independent 
r*- 

sufficiently  localized  solutions,  for  all  Lagrangians  of  the  form  (5.39), 
that  is, 

=  \  VpAp,o  +l  A0,pA0,p  +<4“(Aa,a  :  li(g,a)2;  «,aAo  :  Ao  ;  «>• 


the  following  relation  exists  between  the  spin  and  orbital  angular  momenta: 


. 


o 


B90}{ 


■ 
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(A4. 16)  Lk=(T-l)Sk  . 

Using  the  explicit  expression  (A4.7)  for  we  have 


S  =  U 

4 

ExA 

d3x 

(A4.17) 

L  =  (2V-U) 

« 

► 

ExX 

d3x 

J  =  L  +  S  =  2V 

• 

ExA 

d3x 

where  U  and  V  are  constants,  and  E  -  -V<j>,  is  the  electric  field. 


. 


